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DEVELOPMENT  OF  A  RUN  TIME  MATH  LIBRARY 
FOR  THE  I750A  AIRBORNE  MICROCOMPUTER 


The  purpose  of  this  thesis  was  to  develop  and  evaluate  the 
performance  of  a  run  time  math  library  for  those  architectures  conforming 
to  HIL-STD-1750A,  and  is  intended  to  serve  as  a  benchmark  for  future 
contractor  development.  The  routines  implemented  include  all  circular 
trigonometric  functions  and  their  inverses.  Appendix  A  contains  the 
descriptions  of  pseudo-operations  I've  used  in  explaining  the  design  of 
these  functions,  and  will  prove  useful  to  you  as  you  try  to  follow  my 
logic. 
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the  sponsoring  organization,  for  all  the  help  that  they  gave  me  in  getting 
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Abstract 

This  project  produced  a  run-time  math  library  for  the  MIL-STD-1750A 
embedded  computer  architectures.  The  math  library  consists  of  the 
circular  trigonometric  functions  and  their  inverses.  In  addition,  the  steps 
required  for  the  performance  analysis  of  the  math  library  have  been 
outlined. 

Several  approximation  methods  were  investigated,  but  the  Chebyshef 
Economization  of  the  Maclaurin  series  polynomials,  and  rational 
approximations  derived  from  the  second  algorithm  of  Remes  were 
determined  to  be  the  best  available.  Each  functions  implementation  was 
designed  to  take  advantage  of  features  of  MIL-STD-1750A  architectures. 
The  recommended  test  procedures  will  provide  measures  of  the  average 
and  worst  case  generated  errors  within  each  approximation. 


of  its  avionics  weapon  systems.  Standardization  of  high  order  languages 
and  an  Instruction  Set  Architecture  (ISA)  is  one  way  the  Air  Force  feels 
it  can  reduce  these  costs.  In  the  past,  a  major  cost  contributor  was  the 
proliferation  of  unique  avionics  systems  and  subsystems.  Costs  increased 
with  respect  to-  purchasing  and  inventorying  small-lot  spares  at  many 
bases;  training  technicians  to  maintain  complex  and/or  unique  flight  and 
test  equipment;  developing  and  maintaining  software  development 
facilities;  training  programmers  to  write  application  programs  in  seldom 
used  high  order  languages;  and  in  training  programmers  to  maintain 
software  (especially  operating  systems)  in  seldom  used  machine 
languages.  (1;  8.1) 

MIL-STD-1750A  defines  a  standard  16-bit  instruction  set 
architecture  intended  primarily  for  avionics  weapon  systems.  The  major 
cost  advantage  of  this  standard  ISA  will  come  in  the  form  of  common 
support  software  tools.  An  extensive  set  of  support  software  tools  has 
already  been  developed  and  includes:  a  1750A  assembler/crossassembler; 
a  J73  compiler  with  1750A  ISA  code  generator;  linker/loader  programs; 
and  a  I750A  acceptance  test  program  (1 :  8 A).  Other  cost  benefits  will  be 
realized  through  the  independent  development  of  software  and  hardware, 
(2:D  and  common  maintenance  and  test  equipment.  (3: 1 68) 


reduction,  "in  1978  the  Department  of  Defense  had  in  its  inventory, 
software  written  in  about  150  different  programming  languages.  This 
linguistic  proliferation  increased  maintenance  problems  due  to 
programmer  training  requirements  and  lack  of  support  tools  for  many  of 
the  languages".  (1=6.1)  The  D.O.D.  and  Air  Force  recognized  this  as  a 
problem,  and  they  took  steps  to  correct  it.  D.O.D.  Instruction  5000.31, 
“Interim  List  of  D.O.D.  Approved  High  Order  Programming  Languages", 
states  that  only  approved  languages  may  be  used  for  new  defense  system 
software.  JOVIAL  is  one  of  the  languages  approved  by  this  instruction. 

As  previously  mentioned,  the  development  of  a  standard  ISA  such  as 
MIL-STD-J750A  will  help  reduce  total  life-cycle  costs  of  Air  Force 
avionics  weapon  systems.  This  reduction  will  come  partially  through  the 
use  of  common  support  software  tools,  many  which  have  already  been 
developed.  It  was  also  mentioned  that  one  of  the  support  software  tools 
already  developed  includes  a  JOVIAL  compiler  that  generates  1750A  ISA 
code;  however,  a  math  library  containing  all  the  algebraic  and 
trigonometric  functions  required  by  these  languages  has  not  been 
developed.  The  sponsor  for  this  thesis  is  the  Aeronautics  System  Division, 
Language  Control  Branch.  They  are  the  D.O.D.  JOVIAL  and  ADA  compiler 
validation  site,  and  are  responsible  for  the  development  of  such  libraries. 
Completion  of  this  thesis  will  help  the  Air  Force  reduce  avionics  weapon 
systems  cost  through  the  development  of  a  math  library  for  software 
support  of  all  1750A  systems. 


Problem 


At  the  time  of  this  writing  there  are  no  math  libraries  written  that 
take  advantage  of  the  1750A  instruction  set.  In  keeping  with  the  intent  of 
recent  standardization  policies  of  both  the  D.O.D  and  Air  Force,  the  library 
provided  is  written  in  the  D.O.D  approved  language  JOVIAL.  The  coding  of 
the  library  was  only  a  small  effort  of  this  thesis.  Host  of  the  detail  has 
gone  into  verification  .validation,  and  performance  evaluation  of  the 
product.  Inasmuch,  the  focus  of  this  report  is  divided  into  two  primary 
categories;  software  development  and  software  testing. 

Hath  libraries  are  important  because  they  provide  the  programmer 
several  tools  that  serve  as  building  blocks  for  applications.  Math  libraries 
prevent  programmers  from  having  to  reinvent  the  wheel  each  time  a 
function  is  needed.  Libraries  also  provide  a  means  for  using  functions 
that  take  full  advantage  of  a  particular  computer  architecture  computer 
architecture. 

The  design  of  a  procedure  for  computing  the  value  of  functions  is  not 
mathematically  complete  in  itself.  An  understanding  of  a  computer 
architecture’s  operation  is  necessary  to  insure  that  the  computation  of 
any  given  function  is  as  efficient  as  possible,  while  also  providing  the 
highest  degree  of  accuracy.  Such  architectural  considerations  include: 
word  size;  number  of  bits  in  both  the  exponent  and  coefficient  fields  of  a 
floating  point  number,  the  number  of  integer  and  fraction  bits  in 
fixed-point  numbers,  the  way  mathematical  operations  are  performed  by 
the  architecture,  memory  size  of  the  architecture,  and  execution  time. 
Other  considerations  include  overflow,  underflow,  and  precision.  These 


considerations  for  the  functions  define  the  problem  addressed  by  this 
thesis  effort. 

Scope 

This  effort  was  limited  to  the  design,  code,  and  performance 
evaluation  of  the  circular  and  inverse  circular  trigonometric  functions. 
The  functions  were  included  in  a  math  library  targeted  for  MIL-STD-1750A 
computer  architectures,  and  are  the  ones  typically  found  in  most  FORTRAN 
libraries.  Specifically,  these  functions  include:  sine  (sin),  cosine  (cos), 
tangent  (tan),  cotangent  (cot),  arcsine  (as in),  arccosine  (acos).  and 
arctangent  (atan  and  atan2). 

All  functions  have  been  written  to  either  accept  and  return  double 
precision  fixed-point  values,  or  to  accept  and  return  extended  precision 
floating-point  values.  Floating-point  functions  are  distinguished  from  the 
fixed-point  functions  by  the  name  used  to  invoke  them.  All  floating-point 
functions  have  an  "P  concatenated  to  the  end  of  them;  otherwise,  they 
have  the  same  names  as  those  used  by  the  fixed-point  functions.  For 
example:  the  fixed-point  sine  function  is  invoked  by  using  the  name  “sin“. 
and  the  floating-point  sine  function  is  invoked  by  using  the  name  "sinf“. 

Also  provided  are  performance  summaries  for  each  of  the  functions, 
and  algorithms  that  may  be  used  to  determine  the  polynomial  coefficients 
for  computing  any  of  functions  addressed  by  this  paper.  These  algorithms 
can  be  found  in  Appendix  A,  and  produce  coefficients  that  are  valid  for  any 
nonvector  architecture. 


During  a  design  review  held  in  nay  of  1985,  it  was  made  clear  that 
certain  events  could  cause  overflow  errors  and  underflow  errors,  and 
division  by  zero.  Since  the  functions  are  to  be  used  within  an  embedded 
avionics  weapon  system,  it  is  necessary  that  such  conditions  are  detected 
and  handled  gracefully.  The  consensus  of  opinion  from  all  participants  of 
the  design  review  was,  that  the  functions  should  not  be  aborted,  and  that 
default  values  should  be  returned.  The  error  conditions  and  values  returned 
are  discussed  in  the  individual  design  sections  of  this  thesis.  This 
constitutes  an  important  assumption  on  how  to  handle  such  error 
conditions,  and  bears  further  investigation  before  implementing  on  a 
real-time  system. 

Another  factor  discussed  during  the  design  review  was,  that  there  is 
a  need  for  both  fixed-point  and  floating-point  functions.  Floating-point 
functions  are  more  precise  than  fixed-point  algorithms,  but  take  longer  to 
execute.  Conversely,  fixed-point  functions  don’t  have  the  precision,  but  are 
much  faster.  In  addition,  it  was  mentioned  that  many  avionics  applications 
also  use  what  is  termed  as  pi-radians.  Pi-radians  are  angular  units  of 
measure  expressed  in  terms  of  multiples  of  pi,  and  are  equal  to  radian 
measures  divided  by  pi.  For  example:  180°  is  equivalent  to  3. Ml 596 
radians,  or  1.0  pi-radian.  If  the  algorithms  use  pi-radian  measure  rather 
than  radian  measure,  there  is  no  need  for  overflow  condit  on  checks,  and 
therefore,  a  significant  amount  of  work  is  eliminated  from  the  domain 
reduction  computations  . 


As  stated  earlier,  algorithms  for  both  fixed-point  and  floating-point 
type  functions  have  been  written.  Because  of  time  limitations,  both 
methods  of  unit  measure  could  not  be  implemented,  and  only  the 
fixed-point  algorithms  use  the  pi-radian  metric.  The  assumption  is  that  if 
the  need  arises,  this  thesis  can  serve  as  a  guide  for  implementing  both 
function  types  in  either  unit  of  measure. 

General  Approach 

The  approach  used  during  this  thesis  effort,  is  termed  the 
"  logical  ized"  model  of  a  software  system  development  cycle.  This  approach 
was  considered  a  better  alternative  to  the  more  commonly  used 
"waterfall"  method  of  software  system  development.  The  waterfall 
consists  of  the  apparently  neat,  concise  and  logical  ordering  of  a  series  of 
steps  that  must  be  accomplished  to  obtain  a  final  software  product.  These 
steps  are  performed  in  order  and  include;  systems  analysis,  requirements 
definition,  preliminary  design,  detailed  design,  coding,  testing,  and 
implementation. 

The  logicalized  model  is  similar  to  the  waterfall  model  just 
described,  but  is  more  concerned  with  the  problem  definition  side  of  the 
cycle  (see  figure  1).  This  approach  makes  it  more  useful  for  eliminating 
errors  that  are  typically  occur  during  the  waterfall’s  requirements 
definition  and  design  phases.  Errors  generated  during  these  phases  of  the 
waterfall  model  typically  occur  because  designers  have  a  tendency  to  shift 
between  (abstract)  high-level  design  issues  and  (physical)  implementation 
considerations  .  Thayer  (5:  335-41)  and  Boehm  et  al. 
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Figure  1  Chart  of;  (a)  Waterfall,  and  i.b)  Logicalized 
Model  of  a  Software  Development  Cycle 
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Table  1  Information  Flow  of  the  Logi cal i zed  Software 
Development  Cycle 


(6:  125-33)  made  it  clear  that  these  problems  existed  ,  and  made  the 
point  that  design  errors  not  only  outnumbered  other  errors,  but  that  they 
were  also  more  persistent.  For  that  reason,  more  attention  was  given  to 
the  top-down  decomposition  and  abstract  (logical)  modeling  of  this 
particular  software  system.  Such  a  structured  approach  recommends  a 
dichotomy  between  the  logical  design  issues,  and  implementation  issues. 

The  information  flow  of  a  logicalized  model  is  summarized  in  Table 
l,  and  is  "analogous  to  an  artist’s  conception  of  a  building:  There  is 
enough  information  to  allow  the  customer  and  designer  to  communicate 
and  to  establish  the  buildings  pluses  and  minuses,  but  not  enough  detail  to 
begin  construction.  A  series  of  reviews,  refinements,  and  the  imposition 


of  local  building  ordinances  for  example,  are  necessary  before  that 
construction  can  start.”  (7:14) 

Therefore,  the  approach  taken  for  this  project  was  similar  to  that 
just  described.  The  ASD/Language  Control  Branch  established  the 
requirements  for  a  MIL-STD-1750A  run-time  math  library  written  in  the 
D.O.D  approved  high-order  language  JOVIAL.  During  a  design  review  and 
several  other  meetings,  certain  design  considerations  were  refined.  From 
there,  a  "logical”  model  was  established  as  a  baseline.  This  was 
accomplished  by  using  the  refined  problem  statement,  and  researching  the 
different  methods  for  approximating  the  different  trigonometric  functions. 

The  baseline  model  served  as  a  reference  from  which  all  decisions 
regarding  actual  implementation  could  be  made.  Before  preceding  to  the 
next  phase  of  development  two  such  decisions  had  to  be  made.  These 
decisions  were  to  determine  which  testing  methods  and  which 
performance  evaluation  techniques  would  be  used  after  coding  was 
complete.  These  decisions  determined  what  sort  of  tests  would  catch  all 
possible  errors,  and  determined  what  techniques  could  be  used  to  establish 
a  confidence  level  for  the  final  product. 

Up  until  this  point,  the  abstract  model  has  been  devoid  of  any 
implementation  considerations.  However,  after  it  was  clear  that  the 
abstract  design  was  complete  and  consistent  with  the  requirements,  it 
became  necessary  to  consider  changes  to  fit  the  problem  into  the 
N1L-STD-1750A  environment.  Before  any  changes  could  be  made,  it  was 
necessary  to  complete  the  following  steps;  study  the  architecture  and  ISA 
defined  by  HIL-STD-1750A;  determine  what  resources  were  available, 
such  as  software  support  tools  and  hardware;  and  then  to  learn  how  to  use 


the  available  resources.  From  there,  it  was  possible  to  develop  an  abstract 
model  of  an  implementable  solution.  This  model  took  advantage  of  those 
environmental  factors  that  affected  the  speed  and  accuracy  of  computation 
for  each  function  approximation. 

The  major  subset  of  the  logicalized  software  engineering 
methodology  just  described  is  called  structured  programming.  Structured 
programming  can  be  understood  as  the  decomposition  of  a  problem  in  order 
to  establish  a  manageable  problem  structure.  The  highest  conceptual  level 
represents  a  general  description  of  the  problem,  and  each  level  of 
decomposition  provides  more  detail  into  the  problem.  This  decomposition 
is  carried  out  until  the  problem  is  almost  in  coded  form,  and  is  often 
called  a  stepwise  refinement  of  the  problem.  All  implementation 
considerations  are  left  until  the  lowest  levels  of  refinement. 

The  goals  of  structured  programming  must  be;  to  minimize  the 
number  of  errors  that  occur  during  the  development  process;  to  minimize 
the  effort  required  to  correct  errors  in  sections  of  code  found  to  be 
deficient;  upgrading  sections  when  more  reliable,  functional,  or  efficient 
techniques  are  discovered;  and  to  minimize  the  life-cycle  costs  of  the 
software.  (8=32)  In  addition  it  must  reduce  the  complexity  of  the  problem. 

Structured  flowcharting  is  a  technique  used  to  support  these 
structured  programming  concepts  and  goals,  and  is  “designed  to  reduce 
labels  and  unstructured  branching,  encourage  a  single  entry/single  exit 
approach,  aid  in  the  use  of  top-down  design  techniques,  and  enhance 
modularization.  The  approach  encourages  the  designer  to  conceive  of  the 
system  in  high-level  constructs  and  not  in  terms  of  individual  detailed 
statements.-  (7-116)  The  structured  flowcharting  technique  was  used 


throughout  the  development  of  this  project,  not  only  because  of  the 
aforementioned  reasons,  but  also  for  its  simplicity  and  understandability 
from  a  reviewer  standpoint. 


Sequence _of  Presentation 

This  thesis  addresses  the  the  design  and  performance  evaluation  of  a 
run-time  math  library  that  is  targeted  to  M1L-STD- 1 750A  architectures. 
The  requirements  definition  for  this  problem  has  already  been  discussed 
(Chapter  1  -  Problem/Scope).  The  next  topic  discussed  is  the  theoretical 
development  of  this  thesis  effort  (Chapter  2).  In  particular,  abstract 
design  considerations  for  each  of  the  implemented  functions  is  discussed. 
The  next  area  for  discussion  is  the  detailed  design  considerations  that 
were  made  during  implementation  of  the  library  functions  (Chapter  3).  The 
last  aspects  covered  in  this  report  are  the  test  and  performance 
evaluation  methods  used. 

Appendices  include  algorithms  useful  for  determining  the 
coefficients  of  each  of  the  functions  (appendix  E),  pseudo-code  operations 
used  in  the  structured  flowcharts  (Appendix  A),  source  listings  for  the 
implemented  functions  (Appendix  B),  support  software  developed  in 
conjunction  with  this  thesis  (Appendix  C),  and  the  VAX  VMS  command  files 
required  to  compile  link  and  run  the  developed  product  (Appendix  D). 


The  purpose  of  this  thesis  was  to  create  and  analyze  trigonometric 
functions  developed  for  I750A  architectures.  This  chapter  is  concerned 
with  the  design  theory  of  the  algorithms  used  to  approximate  those 
functions.  Within  the  given  constraints,  the  emphasis  for  each  of  the 
designs  is  to  compute  results  as  quickly  and  as  accurately  as  possible. 

One  way  of  computing  a  value  quickly  is  to  select  an  approximation 
that  converges  rapidly  towards  the  value  of  the  true  function,  f(%).  There 
were  several  methods  of  approximation  that  were  considered;  however,  the 
polynomial  and  relational  approximations  described  by  Cody  and  Waite 
(4-  17-84)  were  found  to  be  the  best.  The  coefficients  given  by  Cody  and 
Waite  were  derived  by  using  Chebyshev  Economization  of  the  Taylor  series 
for  each  function  for  the  approximation  itself,  or  as  a  starting  point  for 
computing  a  rational  approximation  via  the  second  algorithm  of  Remes. 
An  excellent  reference  for  Cnebyshev  Economization  is  Conte  and  de  Boor 
(9;  265-273),  and  an  excellent  reference  for  the  second  algorithm  of 
Remes  is  Ralston.  (10:  301-306) 

Another  means  of  reducing  the  amount  of  processing  time  required  to 
compute  a  result  is  to  take  advantage  of  certain  aspects  of  the  computer’s 
architecture,  as  well  as  the  different  execution  times  for  different 
instructions  within  the  ISA.  For  example,  incrementing  the  exponent  field 
of  a  floating-point  value  is  not  only  faster,  but  more  accurate  than  the 


equivalent  operation  of  multiplying  by  two,  or  examining  the  sign  bit  of  a 
variable  is  faster  than  comparing  the  entire  value  to  zero.  These 
techniques  have  been  used,  and  are  referenced  in  the  design  descriptions 
as  pseudo-operations.  These  operations  are  equivalent  to  those  described 
by  Cody  and  Waite  (A-  9),  and  are  listed  in  Appendix  A. 

The  accuracy  of  an  approximation  may  be  dependent  upon  the  domain 
over  which  the  function  is  approximated.  For  example,  if  the  domain  of  an 
approximation  is  halved,  the  error  may  be  reduced  by  a  factor  of  about 
2*(tt+i)  f0r  an  polynomials  of  degree  n.  (IT-  59)  This  can  be  shown  to  be 
true  for  most  functions,  but  not  all  of  them.  Domain  reduction  has  no 
effect  on  accuracy  in  approximations  of  certain  functions;  however,  it 
still  serves  as  an  excellent  guide  when  designing  an  application.  This  is 
due  to  the  way  computer  architectures  perform  operations  and  store 
mathematical  values  for  floating-point  numbers.  The  most  significant 
bits  of  a  number  are  always  maintained,  and  since’ only  a  finite  number  of 
bits  are  available  to  represent  the  value,  it  is  possible  that  bits  from  a 
fractional  representation  may  be  lost  during  operations  on  targe  numbers. 

Fortunately,  the  trigonometric  functions  lend  themselves  to  domain 
reduction  through  the  properties  of  periodicity,  symmetry,  and 
antisymmetry.  This  allows  function  arguments  to  be  reduced  so  that  a 
more  accurate  approximation  may  be  calculated  than  what  is  possible 
without  argument  reduction.  How  these  properties  are  actually  used  in 
domain  reduction  depends  on  the  function,  and  are  described  in  the 
following  subsections. 


Approximation  Techniques 

The  filL-STD- 1 750A  ISA  doesn’t  call  for  the  implementation  of  the 
elementary  functions  as  standard  instruction  operators,  so  it  is  necessary 
to  design  software  routines  of  optimum  efficiency  to  replace  them.  The 
word  “optimum*  could  be  given  a  variety  of  precise  definitions,  but 
presumably  it  refers  to  an  average  execution  time  and  storage  space. 
Unfortunately,  there  is  no  known  way  to  derive  or  prove  such  an  "optimal" 
design.  For  these  reasons,  the  search  for  the  appropriate  approximation 
technique  was  limited  to  polynomial  and  rational  approximations. 

Some  of  the  most  popular  methods  of  approximation  used  are  called 
Chebyshev  approximations.  Chebyshev  approximations  are  often  referred 
to  as  “minimax"  approximations  because  they  are  used  to  minimize  the 
maximum  “error"  between  the  true  function  f(x),  and  the  approximation 
of  f(%).  However,  these  methods  of  approximation  are  not  without  their 
problems,  and  there  is  a  price,  even  though  it  is  small  one,  to  be  paid  for 
using  them.  For  example,  the  sum-of-squares  of  the  errors  in  a  Chebyshev 
approximation  will  be  higher  than  if  a  least-squares  method  of 
approximation  is  used.  However,  since  Chebyshev  approximations  assure 
that  an  error  is  never  greater  than  a  given  amount,  they  were  selected  by 
this  study. 

Polynomials.  The  first  class  of  approximations  discussed  are 


polynomials,  and  are  the  simplest  of  all  the  classes  of  approximations 
considered.  The  most  important  subclass  of  the  polynomials  is  the  class 


Tfl  (Chebyshev),  and  are  polynomials  not  exceeding  degree  n.  The 
Chebyshev  polynomials  are  especially  important,  and  gave  rise  to  the 
general  concept  of  Chebyshev  "approximations"  discussed  in  the  preceding 
paragraph. 

The  motivation  for  using  Chebyshev  polynomials  over  all  other 
polynomials  is  their  property  of  least  maximum  error,  and  their  error 
behavior  over  the  entire  interval  of  the  approximated  function.  Through 
the  use  of  Theorem  1,  the  Alternation  Theorem  given  below,  Chebyshev 
was  able  to  prove  for  all  the  polynomials  of  degree  n  with  a  leading 
coefficient  of  I,  that  the  Chebyshev  polynomial  divided  by  2  has  the 
least  maximum  error  in  the  interval  [-1,1].  In  other  words,  no  other 
polynomial  of  the  type  mentioned  will  have  a  smaller  error  than 
xft(x)/2ft''.  In  order  for  a  polynomial  Ptt(x)  to  be  considered  a  Chebyshev 
approximation  of  the  function  f(%),  the  theorem  requires  that  the 
maximum  discrepancy  between  f(%)  and  Pa(%)  occur  with  alternating 
signs  at  n+2  points  over  the  interval  [-1,1]. 

Alternation  Theorem:  The  polynomial  Pa  of  degree  in  that  (1) 

best  approximates  f  is  characterized  by  the  existence  of  at 
least  n+2  "points  of  alternation" 

The  other  motivation  for  the  use  of  Chebyshev  polynomials  is  that 
its  generated  errors  are  more  well  behaved  than  the  errors  generated  by 
other  polynomials.  For  example,  approximations,  based  on  the  Maclaurin 
series  whose  ir.;erval  includes  zero,  have  errors  that  are  very  nonuniform 


—  small  near  the  middle,  but  very  large  at  the  end  points.  It  is  more 
desirable  to  use  an  approximation  whose  behavior  is  more  uniform  instead 
of  powers  of  x.  Since,  as  stated  in  Theorem  1.  the  Chebyshev 
polynomials  spread  the  error  over  the  entire  interval,  they  provide  this 
more  desirable  behavior. 

Definition  of.  Ids.  Chebyshev  Polynomials.  The  Chebyshev  polynomials 
form  an  orthogonal  set,  and  are  defined  by  the  following  equation. 

x_  (x)  ■  cos  (ne)  0  ■  arccos  (x)  ( I ) 

n  -  1.  2.  .  .  . 

From  elementary  trigonometry,  cos(n0)  is  a  polynomial  of  degree 
n  in  cos(0),  and  cos(arccos(x))  =  x;  therefore,  it  follows  that  the 
Chebyshev  polynomials  defined  by  Ta(x)  =  cos(n  arccos(x))  are 
polynomials  of  degree  n. 

By  substituting  arccos(x)  for  0  and  xft(x)  for  cos(n  arccos(x)) 

in  the  identity  function  shown  in  equation  (2),  the  recurrence  relation 
defined  in  (3)  is  formed. 

cos  ((cl  ♦  1)0)  ♦  cos  ((n  -  1)0)  =  2  cos  (0)  cos  (n0)  (2) 

Tft+1(x)  =  2x  xft(x)  -  Ta_,  (x)  (3) 


Let  x0  =  1  and  x,  =  x  then  from  the  recurrence  relation 
defined  in  (3),  successive  polynomials  of  greater  degree  can  be  generated 
as  in  column  A  of  Table  2. 


A 

B 

Tn  -  ' 

Tl  =  X 

x2  =  2xT1  -  X0  =  2x2  -  1 

T,  =  2x  x0  -  x,  =  4x3  -  3x 

x°  =  1 

X1  =  X 

II  II 

- 

X**  =  1/2 
x3  =  y4 

(  2x2  -  1  +  1  ) 
(4x3  -  3x  +  Zy 

=  V2(t2  .  t0) 
0-  V4(x3  +  3x,) 

Table  2  (A)  Chebyshev  Polynomials;  (B)  Powers  of  Chebyshev  Polynomials 


By  using  the  results  in  column  A  of  Table  2,  the  powers  of  the 
Chebyshev  polynomials  can  be  found.  That  is,  it  is  possible  to  express  the 
powers  of  %  in  terms  of  xft .  An  example  of  the  powers  of  x  are 
shown  in  Table  2  column  ft  Appendix  E  contains  an  algorithm  that 
generates  both  the  Chebyshev  polynomials,  and  their  powers. 

Chebyshev  Economization.  As  already  mentioned,  the  Maclaurin 
series  can  be  used  to  approximate  many  functions.  In  addition  to  the 
disadvantages  that  have  already  been  mentioned  for  using  this  series  as  an 
approximation,  the  Maclaurin  series  also  converges  very  slowly.  That  is, 
it  takes  several  multiplications  and  additions  to  obtain  a  desired 
accuracy.  One  way  of  obtaining  a  lower  degree  polynomial,  and  still 
maintain  the  desired  accuracy,  is  to  use  a  technique  that  is  called 
“telescoping"  or  "Chebyshev  Economization".  In  other  words,  the 
polynomial  can  be  expressed  in  a  manner  similar  to  that  shown  in  (4). 


To  compute  the  economized  polynomial  approximation  to  the  function 
j(%)  of  absolute  accuracy  e  on  the  interval  (-1,  l],  use  the  following 
procedure  as  outlined  by  Conte  et.al.  (9=  271-272) 


Step  1.  Get  a  power  series  expansion  for  f(%)  valid  on  [-1,1]; 
typically,  calculate  the  Maclaurin  or  Taylor  series  expansion  for  f(x) 
around  %=0. 


Step  2.  Truncate  the  power  series  to  obtain  a  polynomial  as  in  (5), 
which  approximates  f(%)  on  [-1,  I]  within  an  error  where  cn  is 

smaller  than  e,  and  em  is  defined  as  in  (6).  The  result  of  is  the 

SI  A 

!•  maximum  absolute  value,  within  the  interval  [-1,1],  of  the  product  of 

the  first  truncated  coefficient,  %  to  the  power  of  n  ♦  1,  and  the  n  +  I 
derivative  of  the  function  f(x). 


. 


Step  3.  By  making  use  of  a  table  similar  to  that  shown  in  Table  2 
column  B,  expand  the  polynomial  Pfl(%)  into  a  Chebyshev  series  as 
defined  in  (4).  In  other  words,  substitute  the  far  right-hand-side  of  the 
equations  in  Table  2  column  B,  with  the  appropriate  powers  of  % 


r 


•  .  i 

r 


r 
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contained  in  the  polynomial  formed  by  Step  2  of  this  algorithm.  The 
result  is  similar  to  that  shown  in  (7),  but  of  a  greater  degree. 


Step  4.  Retain  the  first  k  +  1  terms  in  this  series,  i.e.  find 
equation  (7).  choosing  k  as  the  smallest  possible  integer  such  that 
equation  (8)  holds  true. 


P  k(%)  =  d0T0(%)  +  .  . .  +  (7) 

E.  *  *  ■  ■  •  *  <*,  *  E  (8) 

Step  5.  Convert  the  result  of  Step  4  into  a  power  series  polynomial 
similar  to  (5).  by  making  use  of  a  table  similar  to  that  in  Table  2  column 
A.  In  other  words,  substitute  the  right-hand-side  values  of  Table  2 
column  A,  into  the  equation  formed  by  Step  4.  Simplify  the  result. 

Rational  Approximations.  In  most  instances,  rational  approximations 
will  generate  a  least  maximum  error  that  is  as  small  or  smaller  than  a 
Chebyshev  polynomial,  and  will  also  cost  less  in  terms  of  the  number  of 
multiplications  and  additions  required  to  compute  them.  Therefore,  they 
deserved  attention  in  this  study. 

As  stated  earlier,  the  approximation  techniques  considered  by  this 
thesis  are  classified  as  Chebyshev  approximations.  These  methods, 
through  their  exploitation  of  Theorem  1,  provide  approximations  whose 
maximum  error  is  less  than  those  generated  by  other  techniques.  There 


are  several  algorithms  that  generate  rational  approximations  that  can  be 
considered  Chebyshev  approximations;  however,  the  ones  that  generate  the 
most  uniform  approximations  are  those  generated  by  the  second  algorithm 
of  Remes.  This  algorithm  is  easily  automated,  and  is  described  in  detail 
by  the  following  subsection. 


Ih£  SfigQQfl 


qL  Remes.  The  method  used  in  this 


description  is  similar  to  that  outlined  by  Ralston  (10:  301-305),  and  is 
summarized  in  Figure  2. 

Let  f(%)  be  a  continuous  function  that  is  to  be  approximated  over  the 
the  interval  [a.  b],  and  let  the  interval  include  the  point  0.0  . 
Furthermore,  let  (9)  equal  the  error  of  any  rational  approximation  of  the 
form  shown  in  (10). 


r  »(x)  -  max  |  f(%)  -  Rmjt(x) 


o 


Step  1  of  the  algorithm  names  the  input  required  for  this  algorithm. 
The  input  value  f(%)  is  the  function  being  approximated.  If  the  algorithm 
is  being  run  on  a  machine  with  higher  precision  than  the  error  for  which 
the  function  is  being  approximated,  then  the  built-in  functions  of  the 
machine  can  be  used  for  f(%).  (f  the  machine  that  the  algorithm  is  to  run 
on  is  of  the  same  precision  for  which  the  approximation  is  to  be  made, 
then  a  reasonable  substitute,  such  as  a  truncated  power  series  that  is  of 
equal  or  greater  precision  than  what  is  being  approximated,  can  be  used. 

The  other  inputs  include;  m,  k.  [a,  b],  and  C0  .  .  .  (^ .  The  values 
m.  and  k  represent  the  degree  of  the  polynomials  found  in  the 
numerator  and  denominator,  respectively.  The  interval  [  a,  b  ]  is  the 
interval  for  which  the  approximation  is  valid,  and  should  include  the  point 
zero,  as  it  will  allow  the  coefficient  b0,  of  the  denominator,  to  always 

be  one.  The  values  C0  .  .  .  Cg  represent  the  first  N  ♦  1  coefficients  of 
the  power  series  polynomial  that  is  being  converted  to  a  rational 
approximation.  The  value  N  represents  the  sum  of  the  degree  of  the 
polynomials  used  in  the  numerator  and  the  denominator  (m  ♦  k). 

The  second  step  of  the  algorithm  is  to  compute  a  series  of  Pade 
approximations  and  their  error  coefficients.  The  Pade  approximations  are 
of  the  form  depicted  in  (1 1 ),  with  the  restrictions  that  OHsm  and 
0  i  j  -  i  i  k.  For  example,  the  sequence  of  Pade  approximations  computed 
for  an  R^2  approximation  would  only  include  R<0>00(%),  R^i.oM- 
R<»u(x),  R<\,(  x.).  and  The  error  of  the  approximations  is 

equal  to  the  first  power  of  %  truncated  from  the  power  series,  multiplied 


by  the  error  coefficients  shown  in  (12).  The  error  calculations  used  would 
only  include:  d*0'0\,  d*1#®2.  d<u>5.  and  d<2,2)5  . 


M  pVl>(x‘i 

R  >  _,(x)  =  -L-i- 

u  <$(*> 


i  =  o _ N  -  I 


=  ,4 


The  coefficients  for  each  of  the  sequence  of  Pade  approximations  are 
computed  using  (13)  and  (M).  Equation  (13)  forms  a  set  of  m  linear 
equations,  which  when  solved,  determines  the  value  of  the  coefficients 
used  in  the  denominator.  Those  values  can  then  be  directly  substituted 
into  the  set  of  equations  formed  by  (14),  and  will  determine  the  value  of 
the  coefficients  for  the  numerator. 


2  CM  .  ,b.  =  0 

>*o  N-s*j  J 


s  =  0.,  1 , . . . N  -  rn  -  1 
(C,  =  o  If  1  <  0,  t’0  =  1 ) 


r  =  0,  1 , .  . ,  m 
(bj-o  if  j  >k  ) 


The  third  step  of  the  Remes  algorithm  is  to  compute  the  economized 
approximation  To  complete  this  step,  it  is  necessary  to  compute 


This  polynomial  can  be  determined  by 


the  Chebyshev  polynomial  x^. 
using  equation  (3)  of  the  previous  subsection.  Once  the  coefficients  of 
thm  are  f°unc!-  then  values  y  from  (15)  can  be  directly  substituted 
into  ( 1 6),  and  thus  solve  w  x).  The  value  tj  in  (15)  is  the  coefficient 
for  u*  in  x^u).  The  rational  approximation  must  also  be  normalized, 
that  is.  the  numerator  and  denominator  must  be  divided  by  b0.  such  that 
b0  will  remain  equal  to  1. 


>vk)  f  ,-,H 

'^N+l  “ 


j  =  o _ N  -  1 


(15) 
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Qk(x)  *  ft  Vi  ^'i(x) 


( 16 ) 


The  final  step  of  the  Remes  algorithm  is  an  iterative  one.  Now  that 
the  initial  approximation  to  the  function  has  been  found,  it  becomes 
necessary  to  find  the  N  +  2  points  of  alternation.  This  can  be  done 
through  interpolation,  or  by  dividing  the  interval  into  several  small  pieces 
and  solving  for  each  point  on  a  division.  This  method  works,  and  all  that 
is  necessary  is  a  little  bookkeeping  to  maintain  a  list  of  the  N  +  2 


. ' 


-&■ 

I 
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points  of  alternation.  This  step  consists  of  the  following  three 
procedures. 


Procedure  1.  Solve  the  system  of  N  +  2  equations  for  the  N  +  2 

unknowns  a0  . ^ ,  bj  . bk  ,  and  rf®  as  shown  in 

expression  (17).  Note  that  E*0>  is  the  magnitude  of  error  in  the 
approximation  at  each  of  the  points  xj<0>  ,  and  for  the  first  iteration 
can  be  assumed  to  be  0. 


2  *(  I  xi®  1 1 

f  t  x|°> )  -  —  —  1-  ,  H)‘E 


Procedure  2.  Find  h0(x)  as  shown  in  (18).  The  function  h0(x)  then 

has  a  magnitude  of  |tf0>|  with  alternating  signs  at  %l.  i=0 . N+  1. 

In  the  neighborhood  of  each  x£<0>,  there  is  a  point  x^15  at  which  h0(x) 
has  an  extremum  of  the  same  sign  as  that  of  f(x)  -  R<(V*(%) 
at  Xj<0>  .  Replace  each  x^®  by  the  corresponding  x£<l>  .  If*  ,  the 
point  at  which  h0(x)  has  its  maximum  magnitude,  is  one  of  the  points 
x^ ,  do  not  perform  procedure  3.  If  not,  replace  one  of  the  points  x^1* 
by  x  in  such  a  way  that  h0(x)  still  alternates  in  sign  on  the  points 

X«> 

Xj 


Procedure  3.  Repeat  procedures  1  and  2  using  %0a' ..... 
from  (17).  This  process  generates  a  sequence  of  rational  approximations 
which  will  converge  to  an  optimum  if  the  initial  extrema  were 
sufficiently  close. 


Design  Considerations 


Sin  and  Cos.  Close  examination  of  figures  3a  and  3b  reveals  a 
relationship  that  can  be  used  to  reduce  the  amount  of  storage  space 
required  by  the  routines  that  compute  the  values  of  these  two  functions. 
The  relationship  is  expressed  in  (19).  Therefore,  by  adjusting  the  argument 
that  is  passed  to  the  Cos  function,  both  it  and  the  Sin  function  can  be 
written  as  short  linkages  to  another  procedure  that  performs  a  majority 
of  the  calculations. 


cos  (a)  =  sin  (  |  oc  |  +  tt/2) 


09) 


The  three  properties  discussed  in  the  preceding  section,  periodicity 
symmetry  and  antisymmetry,  are  illustrated  in  figures  3a  and  3b  .  These 


Figure  3  Periodicity,  Symmetry.  Antisymmetry,  identity 
Properties  for;  (a)  sin;  (b)  cos 


properties  are  exploited  to  reduce  an  argument  to  within  the 
smallest  possible  domain  of  values  accurately  representing  sin(a)  or 
cos(a)  for  all  a  .  By  the  periodicity  property  shown  in  (20)  and  displayed 
in  Figure  3a,  it  is  easy  to  see  that  the  argument  a  can  be  reduced 
to  within  the  interval  [-271.  2tt].  To  compute  sin(cc) .  determine 
p  such  that  oc  =  p  ♦  2nk  .  where  k  is  an  integer.  Then  sin(a)  equals 
sin(p). 


sin  (a)  =  sin  (a  ♦  2n) 


(20) 


Examination  of  Figure  3a  also  illustrates  that  sin  (a)  is 
antisymmetric  about  the  point  zero-pi.  That  is,  sin(a)  equals  -sin(-oc)  . 
By  noting  the  sign  change  of  the  result,  this  property  can  be  used  to 
reduce  an  argument  to  within  the  interval  [0,  2tt]  .  The  interval  [tt. 
2tt]  in  sin(<x)  is  also  antisymmetric  with  respect  to  the  interval  [0,  Tt]. 
Therefore,  by  noting  the  sign  change  of  this  result  and  the  relationship 
shown  in  (21),  a  can  be  set  to  oc  -Tt  .  The  argument  is  then  reduced  to 
within  the  interval  [O.tt]. 

sin  (a  +  tt)  =  -sin  (a)  (21) 

Lastly,  examination  of  Figure  3a  also  illustrates  that  sin(cc)  is 
symmetric  about  the  point  tt/2.  As  shown  in  (22).  the  argument  a  can  be 
reduced  to  within  the  interval  [0.  tt/2]  ,  by  setting  a  to  tt-cc  . 


sin  (a  ♦  Tt/2)  =  sin  (tt/2  -  a)  (22) 


Tan  and  Cot.  As  with  the  sin  and  cos  functions,  close  examination  of 
the  graphs  for  Tan  and  Cot  (figures  4a  and  4b)  reveal  a  relationship  that 
can  be  used  to  reduce  the  amount  of  storage  space  required  to  implement 
both  functions.  The  relationship  is  expressed  in  equation  (23).  Therefore, 
by  flagging  which  routine  was  originally  called,  both  routines  can  serve  as 


short  linkages  to  another  procedure  that  performs  a  majority  of  the 
calculations. 

cot  (a)  =  1  /  tan  (a)  (23) 

The  properties  of  periodicity  and  antisymmetry  are  also  illustrated 
in  figures  4a  and  40.  These  properties  are  exploited  to  reduce  an  argument 
to  within  the  smallest  possible  domain  of  values  accurately  repre¬ 
senting  tan(cx) .  By  the  property  of  periodicity  shown  in  (24)  and  displayed 
in  Figure  4a,  it  is  easy  to  see  that  the  argument  oc  can  be  reduced  to 
within  the  interval  [-tt/2,  tt/2].  To  compute  tan(oc),  simply  determine 
p  such  that  a  =  p  +  Ttk  ,  where  k  is  an  integer.  Then  tan(a)  equals 
tan(p). 

tan  (<x)  =  tan  (a  +  tt)  (24) 

Examination  of  Figure  4a  also  illustrates  that  tan(a)  is 
antisymmetric  about  the  point  zero-pi.  That  is,  tan(a)  equals  tan(-a) .  By 
noting  the  sign  change  of  the  result,  this  property  can  be  used  to  reduce 
the  argument  to  within  the  interval  [0,  tt/2]. 

Knowledge  of  plane  trigonometry  will  allow  reduction  of  the 
argument  to  within  an  even  smaller  interval.  Let  a  be  any  angle  whose 
vertex  is  at  the  origin  of  a  rectangular  Cartesian  coordinate  system.  Also, 
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let  P(x,y )  be  any  point  on  that  angle.  The  variable  %  is  measured 
along  the  horizontal  axis,  and  y  is  measured  along  the  vertical  axis.  Then 
the  function  tan(<x)  is  defined  as  y/x, .  From  this,  it  can  be  shown 
that  for  any  angle  within  the  interval  [tt/4,  n/21,  equation  25  holds 
true. 

Therefore,  by  setting  cc  to  rt/2  -  oc  and  returning  the  reciprocal  of 
tan(a),  the  angle  <x  can  effectively  be  reduced  to  lie  within  the  interval 
[0.  tc/4]  . 

tan  (a)  =  1/tan  (tc/2  -  oc)  (25) 

As  in  and  Acos.  The  equation  p  =  sin(oc)  determines  an  infinite 
number  of  real  values  for  p  within  the  interval  [-1.  11.  Therefore,  the 
inverse,  a  =  sin_,(p) ,  has  many  possible  solutions.  However,  for  sin_1(p) 
to  be  a  true  function,  there  there  is  the  restriction  of  one  a  for  every 
value  p  .  For  this  reason,  it  is  necessary  to  pick  a  range  of  principal 
values  that  will  satisfy  this  restriction.  Figure  5  depicts  the  range  of 
principle  values  defined  for  this  function,  and  is  the  interval  Ht/2,  tt/2]. 

Careful  examination  of  Figure  5  reveals  a  relationship  that  can  be 
used  to  reduce  the  amount  of  storage  space  required  by  the  routines  that 
compute  the  values  of  these  functions.  The  relationship  is  expressed  in 
(26).  From  the  range  of  values  defined  for  sin_1(oO  (see  Figure  5),  and 
the  relationship  defined  in  (26);  the  the  range  of  values  over  which 
cos' ’(a)  is  defined  includes  [0,  ttI.  By  adjusting  the  results  that  are 
computed  in 


Figure  5  Graph  of  the  Inverses  of  bine  and  Cosine 


a  separate  procedure,  both  Asin  and  Acos  can  serve  as  short  linkages  to 
that  procedure. 


cos"1  (a)  =  xr/2  -  sin'1  (a)  (26) 

By  restricting  the  argument  of  the  functions  to  the  interval  [0,  .5), 
the  identity  function  shown  in  (27)  can  also  be  used.  Note  that  by  solving 
for  sin'KcO  in  equation  (26)  and  substituting  the  result  into  the 
left -hand-side  of  (27)  ,  we  come  up  with  the  definition  for  cos'^cc) 
shown  in  (28). 
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Figure  6  Graph  of  the  Inverse  of  Tangent 


sin’1  (<x)  =  tc/2  -  2  sin-1  VO -a)/ 2 
cos-1  (a)  =  2  sin"1  V(l-a)/2 


Atan.  The  equation  (3  =  tan(a)  determines  an  infinite  number  of 
real  values  for  (3  within  the  interval  [-1,  I]  .  Therefore,  the  inverse, 
a  =  tan_1(p)  ,  has  many  possible  solutions.  However  for  tan_,(p)  to  be 
a  true  function,  there  can  only  be  one  a  for  every  value  (3  .  For  this 
reason,  it  is  necessary  to  pick  a  range  of  principal  values  that  will 
satisfy  this  restriction.  Figure  6  depicts  the  range  of  principal  values 
defined  for  this  function,  and  is  the  interval  [-tt/4,  n/4]  . 

Examination  of  Figure  6  illustrates  that  tan_,(oO  is  antisymmetric 
about  the  point  zero.  That  is,  tan"’(<x)  =  -  tan'’(-oc)  .  By  noting  the  sign 
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change  of  the  result,  this  property  can  be  used  to  reduce  an  argument  to 
within  the  interval  [0,  «>]. 

In  a  manner  similar  to  the  reduction  method  described  in  equation 
(24)  of  the  tangent  subsection,  the  argument  can  be  reduced  to  within  the 
interval  [0,1]  .  The  identity  function  depicted  in  (29)  is  valid  for  those 
arguments  greater  than  or  equal  to  one.  This  same  identity  function  can 
be  used  to  reduce  the  argument  to  within  the  interval  [0,  tt/4]  .  This  is 
performed  by  subtracting  from  tt/2  ,  the  arctangent  of  the  inverse  for  the 
reduced  argument. 


Atan(a)  =  xr/2  -  Atan(1/a) 


(29) 


The  final  range  reduction  reduces  the  argument  to  within  the  interval 
[0,  2  -  VT  ]  .  This  is  done  by  making  use  of  the  identity  functions  shown 
in  equations  (30)  and  (31). 


p  =  (a  «/3  -  I)  /  (7T  ♦  a) 


(30) 


Atan(a)  =  tt/6  ♦  Atan(p) 


(31) 


This  chapter  deals  with  the  detailed  design  of  each  of  the  specific 
functions.  As  was  stated  earlier,  two  type  representations  have  been 
designed  and  implemented.  Therefore,  each  of  the  following  sections 
contain  a  discussion  on  the  design  of  each  function's  fixed-point  and 
floating-point  implementation. 

Each  design  has  an  associated  structured  flowchart,  and  each  box 
within  the  flowchart  has  been  numbered  for  ease  of  reference. 
Pseudo-operations  are  used  throughout  each  of  the  flowcharts,  and 
includes  those  defined  by  Cody  and  Waite  (4=  9-10).  Furthermore,  a  few 
additional  pseudo-operations  have  been  introduced,  (see  Appendix  A) 

Although  the  approximation  methods  used  are  those  suggested  by  Cody 
and  Waite,  the  actual  design  implementations  are  significantly  different. 
The  designs  proposed  by  Cody  and  Waite  are  guidelines  for  a  broad  class  of 
computer,  and  weren't  specifically  targeted  towards  a  1750A  architecture. 
Therefore,  the  designs  have  been  tailored  somewhat.  In  addition, 
fixed-point  algorithms  were  designed  to  be  invoked  with  arguments 
expressed  in  pi-radian  measure.  This  metric  was  discussed  in  the 
"Assumptions''  section  of  chapter  one. 

The  coefficients  for  each  of  the  functions  were  either  taken  from 
Cody  and  Waite,  or  are  modifications  of  those  provided  by  Cody  and  Waite. 
These  modifications  are  discussed  in  their  appropriate  subsection. 


Fixed-Point.  In  Chapter  two,  it  was  noted  that  the  identity  function 
for  "Cos",  equation  (19),  could  be  used  to  design  a  procedure  that  will 
compute  the  values  for  both  the  sine  and  cosine  functions.  Indeed,  if  one 
looks  at  Figure  3,  this  concept  can  be  verified.  For  example:  examine  the 
graph  of  cosine  and  note  that  Cos(-tt)  is  equivalent  to  Sin(|  -Tt|  *n/2). 
Therefore,  by  adjusting  the  cosine’s  argument  a  ,  to  its  absolute  value, 
plus  n/2  ,  one  routine  can  be  developed  that  will  compute  the  values  of 
both  functions. 

JOVIAL,  the  high  order  language  used  to  implement  these  functions, 
doesn't  have  the  facility  for  supporting  multiple  entry  points.  As  a  result, 
both  ’’Sin"  and  “Cos”  have  been  implemented  as  shown  in  figures  7a,  7b,  and 
8.  Both  functions  have  their  own  unique  entry  points,  but  both  call  the 
function  “5inCos“  to  compute  the  desired  results. 

Note  that  steps  2,  3,  and  4  of  Figure  7a.  represent  the  identity 
function  for  “Cos”.  The  fixed-point  algorithms  use  pi-radian  measures 
(pi-radians  =  radians/7t);  therefore,  one  half  is  used  in  lieu  of  n/2. 

Since  an  argument's  legal  range  of  values  for  this  fixed-point 
implementation  is  [-1,11  ;  these  functions  have  been  defined  as  having 
one  sign  bit,  one  integer  bit,  and  30  fractional  bits  (-2<%<2).  The  JOVIAL 
compiler  will  not  let  programmers  invoke  this  routine  with  an 
incompatible  type  ;  however,  it  is  still  possible  for  "Cos”  to  generate  a 
value  outside  the  function's  defined  range.  This  could  potentially  cause 
a  problem  on  another  machine,  but  not  on  a  1750A.  Adding  .5 


pi-radians  to  the  argument  whose  absolute  value  (through  programmer 
error)  is  already  within  the  interval  [1.5,  2).  will  cause  an  overflow. 
However,  the  MIL-STD- 1 750A  specifies  that  fixed-point  overflows  will 
result  in  the  most-significant-bits  being  lost,  and  the  least  significant 
retained.  This  feature  of  the  1750A  architecture  is  actually  a  benefit  to 
this  design  implementation.  The  cosine  function  has  a  period  of  2tt 
radians  (2  in  pi-radian  measure),  and  since  an  overflow  causes  the  result 
to  be  a  modulo  of  the  period  length,  the  final  result  is  in  its  reduced  form. 

As  stated  earlier,  the  legal  range  of  values  for  these  functions  is 
[-1,1]  .In  addition  to  the  aforementioned  situation,  it  is  possible  that  an 
adjustment  to  the  cosine  argument  would  result  in  its  being  outside  the 
legal  interval  yet  not  overflow.  No  further  reduction  is  performed; 
however,  as  the  same  potential  problem  exists  in  programs  that  invoke  the 
“Sin"  function.  Therefore,  any  further  reductions  are  handled  in  “SinCos". 

The  variable  “Sign"  shown  in  step  8  of  the  "SinCos"  function  (Figure 
8).  is  what  is  used  to  note  the  sign  changes  discussed  on  pages  26  and  27. 
Since  the  initial  step  of  this  algorithm  is  to  reduce  the  original  argument 
to  within  the  interval  [0,  .5]  .  the  approximating  polynomial  will  always 
compute  a  positive  result.  Therefore,  during  each  phase  of  the  domain 
reduction  process,  it  is  necessary  to  keep  track  of  those  arguments  whose 
functional  values  are,  in  reality,  antisymmetric  with  the  computed  result. 
The  value  returned  to  the  calling  procedure  will  be  the  result  of  the 
polynomial  approximation  multiplied  by  the  value  of  "Sign"  (+1). 

The  variable  "X”  is  the  argument  passed  to  the  "SinCos"  function. 
Since  domain  reduction  on  the  argument  is  required,  and  since  JOVIAL 


treats  formal  parameters  as  read  only,  the  value  of  "X"  is  assigned  to  "F. 
The  variable  “  F"  is  used  throughout  the  domain  reduction  phase,  (step  8) 

The  first  step  of  the  domain  reduction  process  is  to  restrict  "F"  to 
within  the  interval  (0.  2)  .  If  "F"  is  negative,  it  is  antisymmetric  with 
respect  to,  its  corresponding  positive  value.  Examine  Figure  3a  and  note 
that  sin(<x)  equals  -sin(-oO  ;  therefore,  the  sign  change  is  noted 
(Sign  =  -  Sign)  ,  and  “F"  is  set  to  its  absolute  value.  The  built-in 
function  for  absolute  value  was  found  to  give  inconsistent  results,  so  the 
absolute  value  is  found  by:  F  =  -F.  (steps  9  and  10  of  Figure  8) 

The  second  step  of  the  domain  reduction  phase  is  to  reduce  the 
argument  to  within  the  interval  [0,1]  .  Examine  Figure  3a,  and  note  that 
the  interval  [1,2]  (in  pi-radians)  is  antisymmetric,  with  respect  to,  the 
interval  [0.1]  .  Again,  the  sign  difference  must  be  noted  (Sign  =  -Sign), 
and  the  argument  reduced  to=  F  -  1  -  F.  (steps  1 1  and  12  of  Figure  8) 

The  final  phase  of  the  domain  reduction  step,  places  the  argument, 
“F\  within  the  interval  [0,  .5]  .  This  reduction  is  accomplished  by 

exploiting  the  sine  function’s  property  of  symmetry.  Examine  Figure  3a 
again,  and  note  that  the  interval  [0,1]  is  symmetric  about  the  point 

one-half.  That  is,  sin(O)  =  sin(l),  sin(.2)  =  sin(.8) . etc..  Stated 

another  way,  if  "F"  is  greater  than  .5  ;  then  sin(F)  can  be  computed  using 
sin(l-F)  .  Therefore,  “F"  is  reduced  by  setting  F  to  1  -  F.  (steps  13  and 
14  of  Figure  8) 

The  next  step  of  the  algorithm  is  to  compute  the  result  of  a 
polynomial  approximation  for  sin(F)  .  This  polynomial  is  a  "Chebyshev 
Economization"  of  the  Maclaurin  series  for  sine.  However,  if  the  reduced 
argument  is  so  small  that  its  use  within  the  approximation  would  cause  an 


underflow  (i.e.  T  is  less  than  some  epsilon  "eps"),  the  approximation  is 
set  to  rr*F  .  In  systems  using  a  radian  measure  —  sin(F)  approaches  "F" 
as  “F"  approaches  zero.  Therefore,  in  systems  using  pi-radian  measure  — 
sin(F)  approaches  F*tt  as  “F  approaches  zero.  The  lower  limit  is  set 
to  the  value  that  prevents  the  computation  rj’F3  from  underflowing 
(t-j  is  a  coefficient  from  the  approximating  polynomial).  As  with  an 
overflow  condition,  the  1750A  handles  underflows  without  interrupting 
software  processing,  but  sets  the  result  to  zero.  Underflow  may  not 
create  a  problem,  but  the  check  for  underflow  may  eliminate  the  overhead 
associated  with  computing  the  polynomial  approximation,  (steps  15  and 
18  of  Figure  8) 

The  coefficients  used  in  this  algorithm  are  modifications  of  those 

presented  by  Cody  and  Waite  (4=  132).  The  approximation  described  by  Cody 

and  Waite  is  similar  to  the  form  shown  in  (32).  The  "r“  values  represent 

* 

the  coefficients  of  the  approximation,  and  TC0  represents  the  angle  being 
approximated.  These  modifications  were  necessary  because  the  algorithms 
of  Cody  and  Waite  use  arguments  expressed  in  radian  measure  (i.e.  Tte), 
rather  than  pi-radians  as  in  this  implementation. 


P(7T8)  =  r5(TC0)11  ♦  r4(7te)9  ♦  r3(7te)7  ♦  r2(7t0)5  + 
r1(Tt9)3  +  (7T0) 


(32) 
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By  restructuring  (32)  into  the  form  of  that  shown  in  (33).  the  new 

coefficients  can  be  easily  computed.  From  (32),  (33).  and  the  definition 

for  pi -radians;  it  can  be  shown  that  the  new  coefficients,  r'  .can 

be  computed  by:  r'ft  -  ra*7t7  .  Where  "y",  is  the  power  of  the  angle 

expressed  by  the  radian  metric  —  rft(7T8F  in  (34). 

P(tT8)  -  TT1 '  r5  811  ♦  TF9  r4  8*  +  tt7  r3  87  +  tt3  r2  83  + 

n3rl  83  +  tt8  (33) 

A  cursory  examination  of  this  method  of  approximation  shows,  that  it 
would  be  very  inefficient  if  the  algorithm  were  coded  exactly  as  shown 
(32  or  33).  This  problem  is  eliminated  by  making  use  of  "Horner's  Rule".  If 
“G“  is  set  to  the  value  “F2“  ,  the  polynomial  can  be  computed  as  shown  in 
(34)  and  (35).  (steps  16  and  17  of  Figure  8) 

result  -  ((((((r5  G  +  r4)  G  +  r3)  G  +  r2)G  +  rj)G  (34) 

result  -  result  re  -*■  re  (35) 

This  last  method  of  approximation  requires  eight  multiplications  and 
five  additions;  as  compared  to  the  five  multiplications,  five  additions,  and 
five  powers  shown  in  (32).  The  significant  difference  in  the  number  of 
arithmetic  operations,  will  also  prevent  the  loss  of  precision  that 
accompanies  the  computation  of  the  individual  powers  of  “F“. 


r  -  Cody  &  Waite 


0. 1 6666_66560_883  E+0 
0.83333_30720_556  E-2 
0. 1 9840_83282_3 1 3  E-3 
0.27523_97 106_775  E-5 
0.23868_34640_601  E-7 


Power  of  rf 


1.64493-40611-400  E+0 
0.81 174-21707-751  E+0 
0.19074-76226-769  E+0 
0.0261 1-62053-162  E+0 
2.23522-40374-060  E-3 


Table  3  Coefficients  for  Polynomial  approximation  to  Sin 


The  coefficients  for  this  method  of  approximation  were  determined 
from  the  relationship,  r'  -  r  *7t*  ,  and  the  powers  of  tt  implicit  to  (34). 

The  new  values  were  computed  using  the  coefficients  given  by  Cody  and 
Waite,  the  value  of  pi  in  (36).  and  a  machine  of  higher  precision  than  that 
of  the  1750A.  The  new  coefficients  are  shown  in  Table  3. 


pi  =  3.14159_26535_90 


Floating-Point.  The  method  used  in  approximating  the  floating-point 
sine  and  cosine  functions,  is  very  similar  to  that  just  described.  So.  there 
will  be  references  to  the  discussions  of  the  preceding  subsection.  The 
major  difference  between  the  previous  algorithm  and  this  one  is.  that  this 
algorithm  accepts  arguments  expressed  in  the  conventional  radian 
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Figure  9  0)  Sinf ;  b)  Cosf  Structured  Flowchart 


Even  though  these  algorithms  accept  arguments  of  radian  measure, 
they  are  still  quite  different  from  the  algorithms  described  by  Cody  and 
Waite.  Their  algorithm  converts  the  argument  to  pi-radians, 
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m 


SinCosf 


Note  the  naming  convention  used  for  these  functions.  They  are 
distinguished  from  their  fixed-point  counterparts  by  the  letter  T  at  the 
end  of  their  name.  This  convention  is  used  for  all  functions  implemented 
as  a  result  of  this  thesis,  and  will  not  be  elaborated  elsewhere. 

Prior  to  invoking  the  function  “SinCosf",  both  “Sinf"  and  “Cosf"  must 
perform  a  limit  check  on  their  passed  argument.  This  step  is  necessary  to 
insure  a  good  approximation  for  large  arguments.  Due  to  the  way  floating¬ 
point  values  are  stored,  least  significant  bits  are  lost  as  numbers  become 
larger.  It  is  possible  for  a  number  to  be  so  large,  that  it  is  no  longer 
representable  as  a  reasonable  multiple  of  pi.  Therefore,  the  maximum 
limit  used  in  these  algorithms  is  similar  to  that  recommended  by  Cody 
and  Waite  (4:  134).  The  maximum  size  an  argument  can  take  is,  the 
integer  value  tf*2t/2  (where  "t"  is  the  number  of  non-sign  bits  of  the 
floating-point  coefficient).  If  the  argument  is  greater  than  the  maximum 
limit,  each  of  the  respective  functions  will  return  a  value  of  zero.  This  is 
an  area  in  need  of  further  analysis,  and  should  lead  to  an  acceptable 
method  for  handling  exceptions  within  this  avionics  weapon  systems, 
(steps  1  and  2  of  figures  9a  and  9b) 

Before  “Cosf"  can  invoke  the  function  “Sincosf",  it  must  make  use  of 
the  identity  function  described  in  equation  (19).  Since  the  floating-point 
argument  is  still  in  radian  measure  at  this  point,  rr/2  is  added  to  the 
absolute  value  of  the  variable  T\  (steps  3-5) 

The  first  step  of  the  "SinCosf"  function  is  to  initialize  the  variables 
that  will  be  used  during  the  domain  reduction  phases.  The  working  variable 
T",  will  contain  the  eventual  reduced  argument.  Dividing  T"  by  ff 
initializes  it  to  a  multiple  of  pi-radians.  In  the  actual  implementation,  T" 


is  multiplied  by  the  constant  in  (37)  rather  than  divided  by  n  .  This 
method  is  preferred,  because  multiplications  are  more  efficient  than 
division  in  most  architectures,  (step  7  Figure  10a) 

1  /Tf  =  0.3 1 830_9886 1  _838  (37) 

The  variable  T",  initialized  in  step  7,  will  be  used  to  find  the 
integer  portion  of  T*.  Since  JOVIAL  doesn't  include  the  built-in  function 
"INT"  for  extracting  the  integer  portion  of  a  floating-point  number,  other 
features  of  the  language  are  used  to  perform  the  equivalent  action.  This 
method  is  discussed  in  Appendix  A;  however,  at  this  point  the  variable 
must  be  initialized  to  zero. 

As  in  the  fixed-point  algorithm  described  in  the  preceding  subsection, 
it  is  necessary  to  note  any  reduction  process  that  has  an  affect  on  the 
sign  of  the  final  result.  Therefore,  the  “Sign"  flag  must  be  initialized  to 
one.  (step  8  Figure  10a) 

The  first  phase  of  range  reduction,  is  to  reduce  the  argument  to 
within  the  range  of  all  positive  multiples  pi.  If  T"  is  negative,  it  is 
antisymmetric  with  its  positive  counterpart;  therefore,  after  setting  T" 
to  its  absolute  value,  the  sign  change  must  be  noted  (Sign  =  -Sign),  (steps 
9-10  Figure  10a) 

The  next  domain  reduction  phase  reduces  the  argument  to  within  the 
interval  10,  1]  ,  and  is  represented  by  steps  11  through  14  of  the 
“SinCos"  flowchart.  If  the  argument  "F",  is  greater  than  or  equal  to  one,  it 
has  an  integer  portion  defined  within  it.  If  the  integer  portion  of  "F"  is 


Figure  11  Bit  Layout  of  1750A  Floating-Point  Number 


odd,  then  the  argument  lies  in  an  interval  that  is  antisymmetric,  with 
respect  to,  those  arguments  that  lie  within  the  interval  [0,  1]  .  If  it  is 
antisymmetric,  the  sign  difference  must  be  noted  before  proceeding. 

The  description  of  the  ODD  function  (step  12  Figure  10a)  in  Appendix 
A  makes  it  sound  as  though  the  process  is  rather  involved.  To  the 
contrary,  it  is  very  simple  and  fast  as  compared  to  the  method  described 
by  Cody  and  Waite.  This  function  can  be  implemented  by  using  JOVIAL 
specified  tables  that  take  advantage  of  the  floating-point  bit  pattern 
shown  in  Figure  1 1. 

In  a  manner  similar  to  that  described  for  the  ODD  function,  the  "INT" 
function  can  be  simulated.  If  the  exponent  points  to  the  least  significant 
bit  of  a  floating-point  number's  integer  field,  then  bits  zero  through  the 
least  significant  bit,  contain  the  actual  integer.  The  integer  bits  of  the 
argument  "F"  are  copied  into  the  equivalent  bit  positions  of  the  variable 
"Y"  ,  and  then  the  exponent  field  of  7"  is  set  to  the  value  of  the 
exponent  field  contained  in  "F".  Subtracting  "Y"  from  ”F”  will  reduce  the 
argument  "F"  to  within  the  interval  (0,1).  (steps  14  and  15,  Figure  10a) 
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Table  4  Coefficients  for  Polynomial  spproiamation  to  Mnf 


The  last  step  of  the  range  reduction  phase  is  to  place  the  argument  to 
within  the  interval  [0,  .51  .  This  step  is  identical  to  that  described  in 
the  fixed-point  algorithm,  and  involves  comparing  “F"  to  one-half.  If  “F"  is 
greater  than  one-half  the  “F"  is  set  to:  1  -  F  .  (step  17.  Figure  10b) 

The  next  two  steps  of  the  algorithm  are:  to  approximate  the  function 
sin(F);  and  then  adjust  it  by  multiplying  it  by  “F”  and  “Sign”.  These  steps 
are  identical  to  those  described  in  the  previous  subsection,  and  need  not 
be  expanded  again.  Even  though  the  coefficients  used  in  the  approximating 
polynomial  are  different  from  thost  used  by  the  fixed-point  algorithm, 
they  were  derived  in  the  same  manner  as  described  before.  These  values 
are  displayed  in  Table  4. 


Figure  12  a)  Tan;  b)  Cot  Structured  Flowchart 


Fixed-Point.  In  chapter  two,  the  Identity  function  for  "Tan"  and  "Cot", 
equation  (23),  established  that  a  procedure  could  compute  approximations 
for  both  the  tangent  and  cotangent  functions.  A  cursory  examination  of 
Figure  28  will  verify  this.  Therefore,  by  noting  which  function  Is 
invoked,  a  single  procedure  can  be  called  to  do  both  computations. 

Since  JOVIAL  doesn’t  have  the  facility  for  supporting  multiple 
entry-points,  these  functions  have  been  implemented  as  shown  in  figures 
12a,  12b,  and  13.  Both  "Tan"  and  "Cot"  have  their  own  unique  entry-points, 
but  both  call  "TanCot"  to  compute  the  desired  results. 


Note  that  the  "Cot”  function  does  not  implement  the  identity  function, 
as  do  the  “Cos"  or  "Cosf”  functions.  The  approximation  method  used,  is  the 
"Remes”  algorithm  discussed  in  the  preceding  chapter.  Since  this  method  is 
a  rational  approximation  that  uses  a  polynomial  evaluation  for  the 
numerator,  and  another  for  the  denominator;  the  identity  function  can  be 
Implemented  within  "TanCot”  .  "TanCot”  must  choose  which  polynomial 
evaluations  are  to  be  used  as  the  numerator,  and  as  the  denominator.  This 
decision  is  based  on  the  value  of  “Job",  a  flag  passed  to  it  by  "Tan"  and 
"Cot",  (step  2,  figures  12  and  12b)  which  function  has  been  invoked,  one 
routine  can  be  developed  that  will  compute  the  values  of  both  functions. 

Step  5,  of  Figure  13,  is  not  a  computer  operation.  However,  it  is  an 
efficient  means  of  Implementing  the  "Tan"  and  “Cot"  identity  function,  and 
one  of  the  range  reduction  phases.  The  rational  approximation  is  of  the  form 
shown  in  equation  (10);  where  the  polynomial  used  in  the  numerator  is 
represented  by  Pg  ,  and  the  polynomial  used  in  the  denominator  is 
represented  by  Qg.  Each  polynomial’s  evaluation  is  contained  in  a  table 
called  "Poly”:  Pg  is  in  entry  zero,  and  Qg  is  contained  in  entry  one. 
The  "alias”  pseudo-operation,  is  equivalent  to  a  JOVIAL  “Define”  .  It  is  a 
note  to  the  compiler,  that  there  are  two  methods  of  referring  to  each 
table  entry.  The  reason  for  this  may  not  be  intuitively  obvious  at  this 
point,  but  will  become  clear  later. 

As  in  every  function  discussed  previously,  it  is  necessary  to  note  any 
domain  reduction  processes  that  have  an  impact  on  the  sign  of  the  final 
result.  Since  the  final  domain  reduction  phase  of  this  algorithm  reduces  an 
argument  to  within  the  pi-radian  interval  [0,  .25]  ,  computed  results 

are  always  positive.  Therefore,  the  “Sign"  flag  is  initialized  to  indicate  a 


positive  result.  If,  at  any  point  during  the  domain  reduction  step,  it  is 
determined  that  the  true  functional  value  is  antisymmetric  with  what  is 
currently  indicated,  the  flag  is  negated  (i.e.  Sign  =  -Sign),  (step  5,  Figure 
13) 

Steps  6,  7,  and  8  of  the  TanCot  flowchart,  are  where  the  initial 
numerator  and  denominator  are  determined.  If  “TanCot"  was  called  by  "Cot" 
the  "Num"  flag  (indicates  which  table  entry  of  "Poly"  is  the  numerator)  is 
set  to  one,  and  its  counterpart,  "Den"  is  set  to  zero.  This  is  where  the  use 
of  an  alias  may  start  making  sense.  Even  though  this  function  could  have 
been  coded  without  using  aliases;  it  is  more  convenient,  and  more 
readable,  if  the  design  doesn't  have  to  be  concerned  with  "which”  entry 
contains  “what"  polynomial  evaluation.  All  that  is  necessary,  is  to  keep 
track  of  how  each  will  be  used:  either  as  a  numerator  or  as  a  denominator. 
This  is  how  the  identity  function  has  been  implemented.  In  the  event  that 
"TanCot"  is  called  by  “Tan",  “Num"  and  “Den"  are  set  to  zero  and  one, 
respectively. 

The  next  step  of  this  algorithm  is  domain  reduction,  and  the  first 
phase  of  this  step  is  to  limit  the  argument  to  all  positive  values  of 
pi-radians.  If  the  argument  "X"  is  negative,  it  is  antisymmetric  with  its 
positive  counterpart.  The  variable  "X"  is  set  to  its  absolute  value,  and  the 
sign  dichotomy  is  noted  (Sign  =  -Sign),  (steps  9  and  10,  13) 

As  was  the  case  for  "Sin”  and  "Cos",  "Tan"  and  "Cot"  accept  arguments 
expressed  in  pi-radian  measure.  The  legal  range  of  values  accepted  by 
these  routines  is  [-1,1]  ,  and  their  argument's  fixed-point  attributes  are 
the  same  as  those  defined  in  the  "Sin"  and  "Cos”  functions.  Consequently, 


WL-STD-1750A  architectures  will  allow  the  value  of  arguments  to  lie 
within  the  the  interval  (-2.  2).  After  the  aforementioned  domain  reduction 
phase,  arguments  will  be  limited  to  the  Interval  [0.  2)  .  Since  the 
architecture  will  allow  arguments  outside  the  function's  defined  interval, 
it  is  necessary  to  consider  those  exceptions  during  domain  reduction. 
Examine  Figure  4a,  and  note  that  the  period  of  tangent  is  one  pi-radian. 
The  expression,  Tan(X)  =  Tan(X-1)  .  is  a  true  equality,  and  illustrates 
that  function  calls  with  an  argument  in  the  interval  (1,2)  ,  can  still  be 
approximated  if  their  argument  is  reduced  by  one.  (steps  11  and  12. 
Figure  13) 

The  next  domain  reduction  phase  will  reduce  the  argument  to  within 
the  interval  [0,  .5]  .  Examine  Figure  4a,  and  note  that  those  arguments 
falling  in  the  interval  1.5,  I]  are  antisymmetric,  with  respect  to,  those 
lying  in  the  interval  [0,  .5]  .  That  is,  if  X>.5  ,  then  Tan(X)=-Tan(1-X). 
As  an  example,  assume  that  X  =  .75  .  Examine  Figure  4a,  and  note  that 
Tan(.75)=-1  .  Also  note  that,  Tan(1-.75)  =  Tan(.25)  =  1  =  -Tan(.75)  . 
This  example  illustrates  that,  the  argument  can  be  reduced  if  the  sign 
change  is  noted,  and  the  argument  set  to=  X  =  1  -X  .  (steps  1 3  and  14) 

The  final  phase  of  domain  reduction,  reduces  the  argument  to  within 
the  interval  10.  .251  .  Again,  examine  Figure  4a.  Note  that  all  arguments 
lying  in  the  interval  1.25,  .5] ,  are  inversely  related  to  those  arguments  in 
the  interval  [0,  .25]  .  That  is,  if  X>.25  ,  then  Tan(X)  =  1/Tan(.5-X)  . 
Further  argument  reduction  is  accomplished  by  swapping  the  values  of  the 
‘'Poly'’  subscripts,  'Num'  and  "Den',  and  setting  the  variable  'X'  to 
X=.5-X.  (steps  15  and  16,  Figure  13) 


As  stated  earlier,  fixed-point  underflow  and  overflow  will  not  create 
a  problem  for  software  packages  running  on  a  1750A  architecture. 
However,  if  the  argument  is  not  examined  for  such  conditions,  there  will 
be  an  unnecessary  amount  of  overhead  incurred  while  obtaining  an 
approximation  to  these  functions.  Step  17  of  the  "TanCot"  flowchart  is 
intended  to  eliminate  the  unwarranted  overhead.  If  the  argument,  "X",  is 
less  than  some  epsilon,  the  approximation  is  found  by  setting  the 
polynomial  “Pg“  to  X*tt  ,  and  the  polynomial  “Qg"  to  one.  Otherwise, 
each  polynomial  must  be  evaluated  separately.  The  epsilon  used  during 
implementation,  is  the  inverse  of  MAX/tt  ;  where  "MAX"  is  the  maximum 
representable  value  allowed  by  the  attributes  of  the  returned  function.  The 
attributes  of  the  returned  function  are=  one  signed  bit,  12  integer  bits,  and 
18  fractional  bits  (steps  17  and  18,  Figure  13) 

The  process  of  evaluating  the  two  polynomials,  "Pg"  and  "Qg",  is 
similar  to  that  of  P(%J,  in  the  "Sin"  and  "Cos"  implementation.  Each 
polynomial  was  derived  through  the  use  of  the  "Second  Algorithm  of 
Remes",  and  each  are  of  the  form  shown  in  (38)  and  (39),  respectively. 
Coding  the  evaluations  in  the  manner  implied  by  (38)  and  (39).  is  very 
inefficient.  This  inefficiency,  is  eliminated  through  the  use  of  "Homer's” 


Pg  =  P2(rr0)5  ♦  P.(Tre)3  ♦  p0(tt8) 


Qg  =  Q2(7T8)4  ♦  Q.(tt8)2  ♦  Q0 
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Table  5  Coefficients  for  Polynomial  approximation  to  Tan 


Application  of  "Horner’s"  rule,  will  leave  the  polynomials  in  the  form 
shown  in  equations  (40)  and  (41) . 


Pg  =  ([P2  (Tt8)2  +  P,]  (7T0)2  +  P,}  (7T8) 

(40) 

Qg  =  [  Q2  (7T8)2  +  Q,I  (Tt8)2  +  Q, 

(41) 

These  equations  use  the  radian  metric,  ire  .  The  subscripted  values  of 
"P"  and  "Q”,  represent  the  coefficients  determined  by  Cody  and  Waite 
(4-  162).  if  (40)  and  (41)  are  restructured  as  in  (42)  and  (43);  the  angle, 
8  ,  is  a  pi-radian  metric.  The  radian  coefficients,  multiplied  by  their 
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associated  power  of  pi,  determine  new  pi-radian  coefficients.  The  new 
coefficients  are  given  in  Table  5.  (step  19,  Figure  13) 


Pg  =  (  [  (P2'Jt4)82  +  PjTT2]  7T82  +  PjTC}  8 

(40) 

Qg  =  [  (Q27t^)02  +  QjTt2]  82  +  Q1 

(41) 

The  evaluation  of  Tan(X)  is  the  final  step  of  this  algorithm.  The 
result  of  "Pg"  is  stored  in  "Poly(0)“,  and  the  result  of  "Qg"  is  stored  in 
“Poly(0"  .  The  evaluation  to  be  used  as  the  numerator,  as  well  as  the  one 
to  be  used  as  the  denominator,  is  determined  by  the  values  in  "Num"  and 
“Den“.  The  value  returned  is  the  result  of  the  polynomial  division, 
multiplied  by  "Sign"  .  (step  20,  Figure  13) 

Floating-Point.  In  chapter  two,  it  was  noted  that  the  identity  function 
for  "Tan"  and  "Cot”,  equation  (23),  could  be  used  to  design  a  procedure  that 
will  compute  the  values  for  both  the  tangent  and  cotangent  functions.  A 
cursory  examination  of  Figure  28  will  verify  this.  Therefore,  by 
noting  which  function  has  been  invoked,  one  routine  can  be  developed  that 
will  compute  the  values  of  both  functions. 


Figure  14  a)Tanf;  b)  Cotf  Structured  Flowchart. 


Both  “Tanf"  and  “Cotf"  are  logically  identical  to  “Tan"  and  “Cot",  (see 
figures  12  and  14)  The  only  difference  between  the  two  type 
implementations,  is  that  the  floating-point  algorithms  are  invoked  with 
arguments  expressed  in  radian  measure. 

Step  5  of  “TanCotf"  is  similar  to  step  5  of  Figure  13.  The  only 
difference  is  that  the  argument,  "X”,  is  divided  by  tt/2  .  As  mentioned 
earlier,  multiplications  are  typically  more  efficient  than  divisions;  so  the 
division  is  implemented  as  the  product  of  the  argument  and  the  constant  in 
(42). 


tt/2  =  l.57079_63267_95 


(42) 


Steps  6  through  9,  of  Figure  15,  are  identical  to  the  same  numbered 
steps  in  Figure  13.  The  logic,  for  both  cases,  is  identical,  and  won't  be 
expanded  any  further. 

The  division  performed  in  step  5  of  this  algorithm,  expresses  the 
argument  "X"  as  a  multiple  of  tt/2  .  Examine  Figure  4a  and  note  that  every 
other  multiple  of  tt/2  ,  is  the  negative  inverse  of  the  interval  [0,  rr/2]. 
Therefore,  if  X*1  and  its  integer  portion  is  odd;  the  sign  dichotomy  is  noted, 
and  the  table  subscripts  for  "Poly"  are  swapped.  The  subscript  swap, 
changes  which  polynomials  will  be  used  as  the  numerator  and  denominator 
of  this  rational  approximation,  (step  1 1  and  12,  Figure  15) 

Step  13  of  this  algorithm,  is  used  to  express  the  argument  in  pi-radian 
measure  (i.e.  X  *  2/tt  *  1/2  =  X/tt  ■  pi-radians).  Rather  than  dividing  by 
two  or  multiplying  by  one-half,  a  faster  and  more  accurate  method  is  used. 
The  layout  of  a  1750A  floating-point  number  was  presented  in  the  "SinCosf" 
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Table  7  Coefficients  for  Polynomial  approximation  to  T anf 


section.  This  knowledge  can  be  used  to  make  multiplications,  involving 
powers-of-two,  more  efficient.  Remember,  the  exponent  field  of  a 
floating-point  number  is  expressed  as  a  power-of-two;  therefore,  a  division 
by  two  can  be  accomplished  by  subtracting  1  from  the  exponent  field.  The 
subtraction  can  be  made  through  the  use  of  JOVIAL  specified  tables, 
containing  a  signed  integer  field  that  overlays  the  exponent  of  X  . 

Table  6  compares  the  execution  times  of  the  "Exponent  Adjustment" 
method  just  described,  against  floating-point  multiplies.  It's  interesting  to 
note,  that  the  implemented  method  is  almost  three  times  faster  than  a 
straight  multiply. 

The  logic  for  the  rest  of  the  design  is  Identical  to  the  last  six  steps  of 
of  the  fixed-point  algorithm,  and  needs  no  further  explanation.  Two 
implementation  exceptions  are:  the  values  of  epsilon,  and  the  coefficients 
used  in  the  approximations.  Epsilon  is  set  to  a  value;  such  that  Pg/Qg  ,  and 
its  inverse,  will  not  cause  an  overflow,  or  an  underflow.  The  coefficients, 
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for  this  approximation,  were  found  in  the  same  manner  as  described  in  the 
“TanCot"  subsection,  and  are  summarized  in  Table  7. 


ASin  and  ACos  Implementation 


Up  to  this  point,  all  discussions  of  the  detailed  design  section  have 
been  divided  into  two  subsections;  one  concerned  with  fixed-point 
implementation,  and  a  second  concerned  with  floating-point 
implementation.  The  last  two  sections  of  this  chapter  discuss  the  inverse 
trigonometric  functions,  and  their  implementation.  The  fixed-point  and 
floating-point  implementations  have  arguments  that  are  in  the  same 
metric,  and  therefore,  do  not  differ  significantly  in  their  design.  With  but 
two  exceptions,  the  only  difference  between  the  two  type  “ASin”  and 
“ACos“  functions  are  the  coefficients  used  in  their  approximations.  For 
that  reason,  there  is  only  one  design  discussion.  The  lone  exceptions  will 
be  noted  later. 

The  identity  function  that  relates  "ACos"  to  "ASin"  is  shown  in 
equation  (26),  and  is  expanded  upon  in  (27)  and  (28)  .  These  equations 
illustrate  the  point  that  one  procedure  can  be  written  to  compute  the 
approximations  of  both  functions.  The  routine  "ASinCos"  ,  will  return  the 
appropriate  value,  depending  on  which  function  invoked  it.  As  mentioned 
several  times  prior  to  this,  JOVIAL  does  not  support  multiple 
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Figure  1b  ASin  and  ACos  Structured  Flowchart 


entry-points,  and  forces  the  two  functions  to  be  implemented  as  shown  in 
figures  16  and  17. 

Each  function  invokes  the  procedure  "ASinCos"  by  passing  the 
argument,  and  a  flag  indicating  which  function  is  making  the  call. 

The  function  approximation  of  “ASinCos"  is  performed  in  three  steps: 
reduction  of  the  argument  "Y”  to  within  the  interval  [0,  .5]  ;  the 
evaluation  of  sin"’(Y)  ;  and  reconstruction  of  sin_,(Y)  to  the 
representation  of  the  calling  function  sin_1(X) . 

The  first  step  of  argument  reduction  is  performed  by  step  3  in  the 
flowchart  of  Figure  17.  The  sign  change  is  not  noted,  because  the  final 
reduction  is  not  a  multiple  of  tl. 

The  computation  of  sin_,(Y)  is  sensitive  to  error  for  large 
arguments,  especially  for  those  that  are  close  to  1.  Therefore,  careful 
argument  reduction  is  required  to  limit  this  problem.  Steps  5  through  8 
are  designed  to  do  just  that,  and  the  reader  is  encouraged  to  reference 


Cody  and  Waite  (4--  174)  for  a  detailed  discussion  on  the  benefits  of  this 
design  consideration. 

If  the  argument,  reduced  by  step  3,  is  greater  than  .5  .then  a  table 
index  is  set  to:  one  minus  the  value  of  "Job".  That  is  .  if  "ASinCos"  was 
invoked  by  “ACos",  the  table  index  is  given  the  value  "0",  and  if  called  by 
"ASin"  it  is  given  the  value  “1“ .  The  index  value  is  used  to  indicate  which 
path  was  taken  from  step  4. 

If  the  right  path  of  the  flowchart  was  taken,  then  the  next  step  is 
insure  that  a  legal  argument  was  passed  by  a  user.  If  the  argument  is 
greater-than  one  it  is  an  undefined  argument.  The  current  design  sets  the 
return  value  to  the  maximum  representable  value  of  the  implemented 
function  type.  However,  this  is  an  area  requiring  further  analysis,  and 
should  provide  a  better  means  of  handling  such  exceptions. 

Step  8  reduces  the  argument,  such  that  step  13  will  compute  the 
arcsine  of  the  compliment  angle.  Examine  Figure  5  in  chapter  2.  If  the 
results  of  step  13  are  plotted  for  those  arguments  lying  in  the  interval 
[.5,  1]  (as  a  result  of  taking  the  right  path  of  the  flowchart);  they  would 
form  a  curve  that  is  the  the  “mirror"  of  the  arccosine  curve  (i.e.  for  the 
portion  of  the  arccosine  curve  in  the  interval  [.5.  11,  rotate  it  90  degrees 
about  the  horizontal  axis).  This  information  will  prove  useful  in 
understanding  how  sin_1(X)  is  reconstructed  from  sin-1(V). 

If,  in  step  4,  the  argument  is  less-than-or-equal-to  .5  ;  the  left 
path  of  the  flowchart  is  taken.  If  the  the  argument  is  less  than  a 
predefined  epsilon,  then  "Result"  is  set  to  “Y"  .  This  is  one  of  two  steps 
that  differ  for  the  two  implementations.  The  fixed-point  implementation 
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Table  8  ASin/ACos  Fixed  and  Floating-Point.  Coefficients 


converts  “Result"  to  pi-radian  measure  by  dividing  by  pi.  If  the  argument 
is  greater-than  the  epsilon,  then  ”G”  is  set  to  Y2  .  (steps  10-12,  Figure 
17) 

Step  13  computes  the  rational  approximation  of  arcsine,  and  is 
implemented  in  a  manner  identical  to  that  described  in  the  "Tan  and 
Cot"  section  of  this  chapter.  This  is  the  only  other  step  that  differs  in  the 
two  algorithms.  The  fixed-point  implementation  is  expressed  in  pi-radian 
measure  by  dividing  “R(g)M  by  pi.  The  coefficients  of  each  polynomial  used 
in  the  two  type  implementations  are  shown  in  Table  8. 

The  last  step  of  the  algorithm  is  necessary  to  reconstruct  sin'Hx) 
from  the  result,  sin_1(Y)  ,  generated  by  step  13.  The  two  tables  "A”  and 
"B"  are  used  in  this  process,  and  require  the  knowledge  of  which  function 
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invoked  the  "ASinCos"  function,  as  well  as,  which  path  was  taken  during 
the  argument  reduction  step. 

Examine  Figure  5,  and  also  note  the  relationships  expressed  in  (26). 
(27),  and  (28)  .  Earlier,  it  was  mentioned  that,  if  the  right  path  of  the 
flowchart  is  taken  for  argument  reduction  (i.e.  the  argument  "Y"  is  in  the 
interval  [.5,  11  );  the  arcsine  forming  the  compliment  angle  is  computed. 
A  curve  of  the  possible  results  would  “mirror"  the  arccosine  curve  shown 
in  Figure  5  .  If  the  function  “ACos“  is  the  invoking  routine,  and  the  original 
argument  is  positive,  then  the  result  lying  on  the  mirrored  curve  is 
negative,  with  respect  to,  the  arccosine  curve  shown  in  Figure  5  . 
Therefore,  the  results  from  step  13,  are  subtracted  from  zero,  (steps  14, 
15.  and  17;  Figure  17)  If  the  original  argument  is  negative,  the  results  on 
the  “mirrored"  curve  are  7t  less  than  the  arccosine  curve,  and  requires  that 
pi  (1  in  pi -radian  measure)  be  added  to  reconstruct  the  true  function, 
(steps  14-16,  Figure  17) 

If  arccosine  is  being  approximated,  and  the  left  path  of  argument 
reduction  is  taken  (i.e.  the  argument  “Y“  is  in  the  interval  [0,  .5]  ),  the 
results  of  step  13  represent  the  arcsine  curve  shown  in  Figure  5  .  If  the 
original  argument  is  positive,  then  arcsine  curve  is  reconstructed  into  an 
arccosine  curve  by  subtracting  the  results  of  step  13  from  tt/2  .  (steps 
14.  15,  and  17;  Figure  17)  If  the  original  argument  was  negative,  the  sine 
curve  in  the  interval  [0,  .5]  ,  is  tt/2  (  .5  in  pi-radian  measure)  less 
than  that  represented  by  the  arccosine  curve  shown  in  the  interval 
[-.5,  0];  therefore,  the  arccosine  is  constructed  from  the  results  of  step 
13,  by  adding  tt/2  .  (steps  14-16,  Figure  17) 


If  the  arcsine  function  is  being  approximated,  and  the  right  path  of 
argument  reduction  was  taken  (i.e.  “Y“  lies  in  the  interval  [.5,  1]  );  the 
results  of  step  13  are  represented  by  the  “mirrored"  curve.  This  curve  is 
tt/2  (.5  in  pi-radian  measure)  less  than  the  function  being  approximated, 
and  requires  that  tt/2  be  added  to  reconstruct  it  to  the  true  function.  If 
the  left  path  was  taken,  the  results  of  step  t3  are  already  represented  by 
the  arcsine  curve.  If  the  original  argument  was  negative,  then  the  results 
of  step  18  are  antisymmetric,  with  respect  to,  the  true  value  of  the 
function,  and  are  complimented. 


This  chapter  is  concerned  with  describing  the  methodology  used  for 
determining  the  correctness  and  performance  qualities  of  the  implemented 
functions.  Due  to  problems  in  the  availability  of  hardware  and  the 
associated  support  software,  the  testing  and  performance  evaluations  are 
somewhat  limited.  Hardware  became  available  towards  the  middle  of  the 
thesis  effort,  but  software  tools  used  for  development  were  incompatible 
with  those  required  by  the  available  J750A.  The  loader  used  by  the 
available  1750A  equipment,  expects  files  of  a  different  format  than  what 
is  created  by  the  software  development  tools.  Rather  than  developing  a 
new  loader,  a  routine  was  written  that  converts  load  modules  into  a 
format  required  by  the  I750A  loader.  The  reformatting  procedure  is  listed 
in  Appendix  C. 

Another  problem  that  had  to  be  overcome  before  testing  and 
evaluation  could  be  considered,  was  the  availability  of  input/output  (I/O) 
routines.  Without  I/O  routines,  further  considerations  for  testing  would  be 
fruitless.  No  I/O  packages  were  available,  and  as  a  consequence,  had  to  be 
created.  This  delayed  testing  efforts  considerably,  as  an  I/O  routine  had 
to  be  developed  with  the  use  of  the  MIL-STD-I750A  standard  ISA,  rather 


than  with  a  high-order  language.  The  I/O  package  developed  is  listed  in 
Appendix  B,  and  is  only  capable  of  writing  to  a  user  console. 

Performance  analysis  requires  the  comparison  of  1750A  results,  with 
those  generated  on  a  machine  of  higher  precision.  Unfortunately,  this 
requirement  made  the  newly  created  I/O  routine  insufficient  for  this  task. 
An  available  console  driver  has  a  routine  that  writes  user  specified  areas 
of  I750A  memory  to  magnetic  disk.  By  storing  a  function's  results  in  a 
specified  area  of  1750A  memory,  the  test  results  can  then  be  dumped  to 
disk  for  an  eventual  upload  to  a  VAX  11/780A.  The  results  are  then 
available  for  input  to  the  different  software  test  packages.  However,  the 
record  format  of  the  1750A  memory  dump  is  not  in  a  friendly  format,  and 
must  be  converted  to  a  readable  form.  At  the  time  of  this  writing,  a 
routine  for  making  the  disk  file  readable  is  not  completely  debugged. 
However,  it  is  at  a  point  where  it  could  be  completed  by  another 
programmer. 

The  aforementioned  problems  have  limited  the  amount  of  time 
available  for  designing  extensive  test  procedures.  Therefore,  validation, 
verification,  and  performance  analysis  is  confined  to;  manual  static 
analysis  methods,  critical  value  testing,  and  measurement  of  each 
algorithms  generated  error. 


To  most  people,  manual  static  analysis  is  called  "desk  checking"  . 
Static  analysis  involves  the  search  for  any  inconsistencies  between  design 
tools  (i.e.  flowcharts),  design  details  (chapter  3),  program  headers,  and 
program  comments.  This  method  is  useful  for  finding  errors  caused  by  the 
translation  of  design  into  code,  as  well  as  possible  design  errors.  An 
inconsistency  may  indicate  potential  problems.  This  methodology  was 
used,  and  all  inconsistencies  that  were  found  were  resolved. 


Critical  Value  Testing 

Critical  value  testing  is  an  attempt  to  "break"  the  software,  and 
requires  the  selection  of  specific  arguments  that  could  possibly  cause 
problems.  A  knowledge  of  each  of  the  algorithms  is  required  to  select 
proper  arguments.  Individual  test  cases  are  not  listed  here,  but  the  reader 
may  find  specific  information  by  examining  the  test  procedures  listed  in 
Appendix  B. 

it  is  possible  to  generalize  the  tests  performed  without  listing  the 
specific  test  cases.  Potential  test  arguments  are  those  whose 
intermediate  results  could  generate  an  overflow  or  underflow,  or  are 
arguments  lying  in  the  fringe  of  computational  abnormality.  These 


arguments  will  help  detect  problem  areas,  and  will  give  an  indication  as 
to  how  robust  each  function  is. 


In  addition,  arguments  that  test  each  path  of  the  algorithm  have  been 
selected.  Path  testing  is  limited  to  insuring  that  every  path  of  an 
algorithm  is  tested,  and  does  not  imply  that  every  possible  path 
combination  is  taken. 


As  was  mentioned  in  the  introduction  of  this  chapter,  screen  output 
to  the  user  console  and  hard  copies  of  computed  results  are  insufficient 
for  performance  evaluation.  Their  use  would  imply  a  visual  comparison  of 
generated  results  against  published  tables.  Such  a  technique  limits  the 
number  of  comparisons  that  could  be  made,  and  would  cause  doubt  as  to 
the  credibility  of  the  comparisons.  At  best,  it  would  provide  a  good 
feeling  for  the  quality  of  each  function's  performance.  Therefore,  it  is 
better  to  automate  the  process  completely,  and  compare  the  generated 
results  against  another  machine  generated  standard. 

The  performance  evaluation  of  the  functions  involves  the  computation 
of  two  important  statistics1  the  maximum  relative  error  (MRE),  and  the 
root-mean-square  relative  error  (RE).  Their  values  are  determined  through 
the  use  of  (43)  and  (44),  where  F(x)  is  the  test  result  and  f(x)  is  the 
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comparison  value  generated  by  the  same  extended-precision  function  call 
written  for  the  VAX  1 1/780. 
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This  method  of  error  checking  is  an  automatic  tabular  comparison, 
where  the  VAX  routines  serve  as  the  accepted  standard.  The  test  routine 
tests  densely  packed  samples  of  evenly  spaced  arguments  spread 
throughout  [-37t.  3tt]  for  floating-point  algorithms,  and  [-1,1]  for 
fixed-point  algorithms.  When  regenerating  arguments  within  the  test 
modules,  it  is  important  not  to  introduce  unnecessary  errors.  This  means 
that  arguments  in  the  VAX  should  have  its  lower  order  bits  padded  with 
zeros.  The  most-significant  bits  must  be  equivalent  to  the  number  of  bits 
in  the  1750A  argument,  and  no  extra  precision  should  be  introduced. 

The  method  of  argument  generation  just  described  is  recommended  by 
Cody  (12:  762),  and  is  the  method  used  at  the  NASA  Lewis  Research 
Center.  This  method  is  preferred  to  a  random-number  test  because  it 
measures  the  relative  error  throughout  an  entire  interval.  Using  densely 
packed  arguments  also  gives  valuable  insight  to  problems  of  different 
argument  ranges.  If  the  evenly  spaced  interval  is  set  to  a  power  of  two 


(representable  on  both  machines),  and  is  not  less  than  the 
least-significant  bit  of  the  1750A  argument,  then  an  initial  argument  can 
be  chosen,  such  that,  zero  padding  will  only  have  to  be  performed  once. 
For  example,  if  an  initial  floating-point  argument  is  -3. 1415  and  the 
chosen  interval  is  2'2  ,  the  second  argument  will  be  -3.M15  +  2"2  . 
Additional  padding  is  not  necessary,  because  "carries"  are  cascaded 
forward  and  do  not  increase  the  number  of  most-significant  bits  in  the 
next  argument.  Arguments  used  in  the  function  calls  on  both  machines 
must  be  the  same,  and  must  be  generated  in  the  same  order. 

Extra  care  is  needed  while  reading  the  1750A  results  from  disk.  Each 
of  the  1750A  results  are  stored  in  an  unformatted  file,  and  must  be  read 
into  a  binary  record.  This  record  is  moved,  bit-by-bit,  to  a  variable  of  the 
appropriate  type  (VAX  11/780  fixed-point  or  floating-point).  The 
bit-by-bit  manipulation  is  accomplished  through  the  use  of  JOVIAL 
specified  tables,  and  prevents  conversion  errors  associated  with 
formatted  input. 

Before  a  comparison  of  the  two  results  (one  from  the  1750A,  and  the 
other  from  the  VAX)  can  be  made,  the  results  generated  within  the  test 
module  must  be  reduced  to  the  same  precision  (same  number  of 
most-significant  bits)  as  those  from  the  I750A.  The  precision  reduction 
gives  a  rounded  result  that  can  be  used  to  determine  the  MREand  RE,  and 
will  give  a  meaningful  interpretation  to  the  inherited  error  of  the  1750A 
functions. 


The  purpose  of  this  thesis  was  to  develop  and  to  do  performance 
evaluation  on  a  run-time  math  library  developed  specifically  for 
MIL-STD-1750A  architectures.  The  library  consists  of  the  floating-point 
implementation  of  several  algebraic  functions.  Performance  evaluation 
was  the  major  effort  of  this  thesis,  but  not  in  the  manner  intended. 

Function  approximations  are  accomplished  through  the  use  of  either 
Chebyshef  or  rational  approximations.  The  two  different  approximation 
methods  were  discussed  in  chapter  two,  and  are  useful  in  understanding 
certain  design  considerations.  The  values  of  each  polynomial’s 

coefficients  were  derived  by  (or  were  modifications  of  those  derived  by) 
Cody  and  Waite.  (4=  17-84)  However,  the  implementation  designs  are 
significantly  different  from  those  suggested  by  Cody  and  Waite.  The 
primary  difference  between  the  implemented  designs  and  those  suggested 
by  Cody  and  Waite,  are  the  methods  of  argument  reduction  required  of  each 
function. 

Performance  evaluation  turned  out  to  be  the  major  effort,  but  not 
because  of  extensive  or  elaborate  testing  of  the  library  functions.  Host  of 
the  effort  involved  overcoming  the  following  problems: 


1. )  There  were  several  compiler  bugs  in  the  original  1750A  compiler 
used.  Assembly  listings  had  to  be  reviewed,  in  order  to  verify  each 
compilation  of  the  source  code. 

2. )  The  use  of  a  simulator  for  performance  evaluation  was  ruled  out 
because  of  the  limited  number  instructions  that  could  be  simulated,  its 
inability  to  simulate  the  use  of  floating-point  data,  and  the  relative  speed 
at  which  results  were  calculated.  The  simulator  also  lacked  a  facility  for 
writing  results  to  mass  storage.  Storage  of  results  on  an  external  device 
is  necessary  for  input  to  software  test  packages. 

3. )  A  new  compiler  and  linker  was  introduced  near  the  midpoint  of  the 
thesis  effort,  and  required  a  long  learning  curve  in  order  to  use  them. 

4. )  Once  a  1750A  machine  became  available,  it  was  determined  that 
all  its  support  software  was  intended  for  use  with  files  created  by  the 
old  compiler  and  linker. 

5. )  Rather  than  use  a  compiler  and  linker  that  had  several  deficiencies, 
or  write  a  new  loader  routine,  it  was  decided  to  write  a  support  tool  that 
would  convert  load  modules  into  a  format  expected  by  the  available  loader. 

6. )  The  reformatting  program  required  the  use  of  JOVIAL  and  its 
specified  table  features.  It  also  required  the  use  of  FORTRAN  routines  to 
perform  the  I/O  of  source  and  target  files.  The  FORTRAN  and  JOVIAL 
interfaces  did  not  operate  as  expected,  and  the  use  of  COmON/COMPOOL 
areas  wouldn't  work.  This  required  parameter  passing  between  the 
routines,  and  the  documentation  for  this  type  of  interface  was  very 
inadequate;  however,  the  problems  were  eventually  resotved. 


7. )  The  reformatting  tool  was  written  for  use  on  a  VAX  1 1/780.  It 
was  assumed  that  the  JOVIAL  compiler  was  free  of  bugs  for  a  VAX  target. 
However,  when  the  reformat  routine  was  being  debugged,  it  was 
discovered  that  JOVIAL  table  names  could  be  overlayed,  but  corresponding 
table  items  weren’t  overlayed  with  them.  This  problem  took  a  long  time 
to  discover,  and  an  additional  amount  of  time  to  design  around. 

8. )  I/O  routines  have  not  been  written  for  the  1750A,  and  had  to  be 
developed.  These  routines  are  only  capable  of  writing  to  a  console  screen. 

9. )  Screen  output  is  insufficient  for  generating  the  thousands  of 
results  that  would  be  needed  during  testing  and  evaluation,  so  another 
means  of  capturing  the  data  had  to  be  developed.  Due  to  the  lack  of  time 
and  inexperience  in  the  internal  I/O  communications  techniques  of  the 
1750A  hardware,  development  of  a  disk  I/O  routine  was  not  a  feasible 
alternative.  It  was  determined  that  results  could  be  stored  in  specific 
locations  of  memory,  and  then  an  available  console  routine  could  be  used 
to  write  the  information  to  disk.  An  additional  problem  was  encountered 
when  it  was  discovered  that  the  record  format  of  the  disk  files  is  not  in  a 
VAX  friendly  format,  and  another  routine  had  to  be  written  to  unpack  the 
stored  results. 

These  problems  limited  the  scope  of  this  thesis  effort  to  developing 
the  following;  designs;  code  that  is  free  of  syntax  errors;  the 
development  of  command  files  for  compiling,  assembling,  and  linking 
routines  written  for  the  1750A;  tools  for  formatting  load  modules  that 
are  capable  of  being  loaded  into  a  Sperry  1631  implementation  of  the 
MIL-STD-1750A;  and  tools  that  unpack  test  results  stored  on  an  RT/11 


formatted  floppy  disk.  Generic  test  algorithms  are  provided,  but  are  not 
written  in  a  high-order-language.  They  provide  the  basic  structure  for 
critical  range  testing,  and  a  means  of  evaluating  and  measuring  each 
functions  performance. 


The  products  produced  by  this  thesis  effort  are  at  point  where  design 
of  the  intended  performance  evaluation  can  begin.  All  the  groundwork  has 
been  provided,  and  should  be  adequate  for  someone  to  continue  the  effort. 
Many  of  the  aforementioned  problems  have  been  resolved,  and  support  tools 
and  command  files  are  provided  to  shorten  the  learning  curve  that 
follow-on  programmers  will  have  to  experience. 

The  following  recommendations  should  be  considered  if  this  effort  is 
continued. 

1. )  If  the  effort  is  limited  to  the  use  of  JOVIAL,  an  analysis  should  be 
made  for  determining  how  to  handle  exceptions  detected  at  run  time. 
Exceptions  include  arguments  outside  legally  defined  limits. 

2. )  Since  Ada  has  features  for  exception  handling,  all  the  library 
functions  should  also  be  developed  and  implemented  in  Ada. 

3. )  Another  point  may  be  in  favor  or  using  Ada  is  that  it  also  allows 
the  creation  of  generic  packages  and  subprograms.  The  generic 


subprograms  define  a  template,  and  generic  parameters  provide  the 
facility  for  tailoring  the  template  to  fit  a  particular  need  at  translation 
time.  In  other  words,  one  subprogram  could  provide  calculations  for  both 
fixed-point  or  floating-point  arguments,  based  on  how  it  is  used  at 
compile  time.  Because  a  generic  package  would  not  be  able  to  take 
advantage  of  the  specific  hardware  functions  unique  to  floating-point  and 
fixed  point  routines,  this  may  result  in  a  degradation  of  performance. 

4. )  Initially,  it  was  discussed  that  all  the  math  library  routines  should 
be  written  in  both  JOVIAL  and  Ada  with  the  intent  that  a  comparative 
evaluation  could  be  done  on  the  two  languages.  Unfortunately,  an  Ada 
compiler  targeted  to  the  1750A  is  not  yet  available.  When  a  compiler 
does  become  available,  it  is  recommended  that  a  new  Ada  math  library  be 
developed  and  this  comparative  evaluation  be  performed. 

5. )  The  compiler  problems,  mentioned  above,  should  be  corrected,  and 
1750A  architectures  and  associated  support  software  should  be  acquired 
before  more  time  is  allocated  to  the  effort. 


The  following  pseudo-operations  were  used  in  describing  the 
implementation  designs  of  the  different  mathematic  functions. 

ADX(X.N):  augments  the  integer  exponent  of  a  floating-point 

representation  of  X  by  N.  This  scales  the  argument  X  by  2N . 
For  example, 

ADX(1 .0,2)  =  4.0 

FIX(X):  returns  the  fixed-point  representation  of  the 

floating-point  value  X  .  This  operation  requires  explicit 
conversion  in  JOVIAL. 

FLOAT(X):  returns  the  floating-point  representation  of  the 

fixed-point  argument  X.  This  operation  requires  explicit 
conversion  in  JOVIAL. 

ODD(X):  determines  whether  the  argument  X  is  odd.  For  an 

integer,  the  least-significant  bit  is  checked  directly.  For  a 
floating-point  number,  the  integer  portion  is  checked.  A 
description  of  the  floating-point  process  for  this  determination 
is  given  below. 


Figure  18  Bit  Layout  of  1750A  Floating-Point  Number 


To  determine  whether  the  integer  portion  is  odd,  knowledge 
of  the  internal  representation  of  the  1750A  floating-point  number  is 
necessary.  The  argument  X  is  a  JOVIAL  specified  table  item  that 
makes  the  components  shown  in  Figure  18  easily  accessible.  Within 
this  table  is  an  integer  item  that  overlays  the  exponent  field  of  X. 
This  exponent  field  is  the  tool  needed  to  check  whether  the  integer 
portion  is  odd  or  even.  Since  X  has  a  value  of  one  or  greater,  and 
all  floating-point  values  are  normalized,  the  exponent  can  be  used  to 
point  to  the  least  significant  bit  of  the  integer  field.  Because  X  is 
positive,  a  one  in  the  least  significant  bit  would  indicate  the 
integer  portion  is  odd.  A  limit  on  the  maximum  value  of  the 
coefficient  has  been  imposed  by  the  functions  that  use  this  routine. 
This  limit  prevents  the  least-significant  bit  of  the  integer  portion 
from  falling  in  the  exponent  or  “LSB“  area  of  the  floating-point 
coefficient  (see  Figure  18). 


Since  the  1750A  architecture  requires  that  all  floating-point 
values  be  normalized,  the  most-significant  bit  is  in  the  first  bit 
position  following  the  sign  bit.  The  decimal-point  is  assumed  to  be 
positioned  immediately  behind  the  sign  bit,  but  immediately  in  front 
of  the  most-significant  bit.  The  exponent  represents  a  power  of 
two;  therefore,  if  e  represents  the  value  of  the  exponent  field,  the 
value  of  the  floating-point  number  is:  coefficient  *  2e  . 
Equivalently,  it  is  obvious  that  the  decimal-point  floats  t  places  to 
the  left  if  t  is  negative,  or  s  places  to  the  right  if  positive. 

Knowledge  of  how  floating-point  numbers  are  stored  can  be 
used  to  determine  whether  the  integer  portion  of  a  number  is  odd. 
The  following  example  gives  an  explanation  of  the  process. 

Given  the  following  machine  representation  of  a 
floating-point  number,  determine  whether  its  integer  portion 
is  odd.  In  the  example  below,  the  decimal-point  was  inserted 
only  for  clarity. 

0. 1 1 0000000000000000000000000000 1 0000000000000000 

Since  the  sign  bit  of  the  exponent  is  zero,  the  value  of 
the  coefficient  is  positive.  The  following  two  numbers  are 
summed  together  to  determine  the  value  represented  by  this 
coefficient: 


The  exponent  field  is  in  bold  text,  and  has  the  value 
one.  Therefore,  the  value  of  this  floating-point 
representation  is,  the  coefficient  (.75)  multiplied  by  two 
to-the-power-of  the  exponent  ( 1 ),  or  1.5. 

.75  *  21  =  1.5 

Another  way  to  compute  the  result  is  to  shift  the 
decimal-point  in  a  direction  as  indicated  by  the  exponent. 
The  exponent  in  this  case  is  +1  ,  so  the  decimal-point  is 
shifted  one  position  to  the  right.  The  number  can  then  be 
computed  in  a  similar  manner  as  described  above. 

This  last  method  demonstrates  how  to  determine 
whether  this  example  is  even  or  odd.  If  the  decimal-point 
is  shifted  1  position  to  the  right,  this  number  will  have  1 
integer  bit  and  38  fractional  bits.  The  integer  bits  always 
occupy  the  left-most  position  of  the  number.  If  the  exponent 
is  thought  of  as  a  pointer  from  the  left-most  side  of  the 
number,  the  least-significant  integer  bit  can  be  found.  The 
exponent  in  this  example  points  to  bit  position  one.  Since 
the  bit  is  set  to  1,  this  example's  integer  value  is  odd.  I 


(X):  return  the  integer  portion  of  the  floating-point 

argument  X.  The  description  ODD(X)  given  above  determines  the 
least-significant  bit  of  the  integer  portion  of  the  floating-point 
argument.  This  is  used  to  extract  the  entire  integer  portion  of 
the  argument  (bits  0  through  the  least  significant  bit). 


AD-A163  993  DEVELOPMENT  OF  A  RUN  TIME  HATH  LI8RARV  FOR  THE  1739A 
AIRBORNE  MICROCOMPUTERS)  AIR  FORCE  INST  OF  TECH 
MRIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENGINEERING 
UNCLASSIFIED  S  A  HOTCHKISS  DEC  83  8FIT/GCS/NA/85D-5  F/G  9/2 


DATE:  18  JULV  1885 

UERSION:  1.0 

NATE:  Cos 

rmJLE  NUneER:  1.0 
DESCRIPTION: 

This  function  is  Invoked  to  cosputa  tha  fixad-point  value  for  tha 
cosina  of  an  angle  that  has  baan  axprassad  in  pi -rod ions.  A 
pi-radian  can  ba  axprassad  as  a  value  bataaan  -1.0  and  1 .0,  and 
ehose  value  ahan  aultipliad  by  pi  is  equivalent  to  an  angla 
axprassad  in  radian  naasura.  Fixad-point  has  tha  advantage  of 
speed  of  coaputation  oust  floating-point  algorithms,  and  using 
pi-radians  further  axploits  this  diffaranca  by  simplifying  tha 
critical  stop  of  ranga  raduction.  Tha  ranga  raduction  effort  is 
parforaad  in  SlnCos,  and  is  usad  to  reduce  tha  passed  arguaent 
into  tha  interval  <-.5  pi,  .5  pi).  This  range  raduction  insures 
aaxiaua  accuracy  of  the  approximated  function.  Both  the  argument 
*Xx*  and  the  returned  value  ‘Cos'  are  expressed  in  1750A  double- 
precision  fixad-point  representation.  Both  have  one  sign  bit,  one 
integer  bit,  and  30  fractional  bits.  Tha  actual  value  returned  is 
computed  by  the  function  'SlnCos' .  Because  the  identity  function, 
Sin<x)  ■  cos(x  -  pi/2),  'SlnCos'  can  be  invoked  to  compute  both 
sine  and  cosina  values. 

PASSED  UAR I HBLES : Xx  -  The  pi-radian  measure  for  ehich  cosina  computed. 

The  argument  is  in  double  precision  fixed-point. 
RETURNS:  Cos  -  the  computed  value  in  double  precision  fixad-point 
MOOULES  CALLED:  SlnCos 
AUTHOR:  Copt.  Steven  A.  Hotchkiss 

HIST0RV:  This  project  aas  undertaken  as  a  thesis  project  for 

partial  fulfil  leant  of  requirements  for  an  MS  degree 
In  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  AS0 
Language  Control  Branch,  Wright  Patterson  RFB,0h. 


Cosine  Procedure 


3 

3 


■ 


DflTE:  18  JULV  1985 

UERSION:  1.0 

NRME:  Sin 

MODULE  NUMBER:  1.0 
DESCRIPTION: 

This  function  is  invoked  to  coapute  tha  fixad-point  value  for  tha 
sina  of  an  angla  that  has  baan  axprassad  in  pi-radians.  R 
pi-radian  can  ba  axprassad  os  a  value  bataaan  -1.0  and  1.0,  and 
ehose  value  ahan  auitipliad  by  pi  is  aqui volant  to  an  angla 
axprassad  in  radian  aoasura.  Fixad-point  has  tha  odvantoga  of 
spaed  of  coaputation  ova r  f looting-point  algorithas,  and  using 
pi-radians  furthar  axploits  this  diffaranca  by  siaplifying  tha 
critical  stop  of  ranga  raduction.  Tha  ranga  raduction  affort  is 
parforaad  in  SinCos,  and  is  usad  to  raduca  tha  possad  orguaant 
into  tha  intarual  <-.5  pi,  .5  pi).  This  ranga  raduction  insuras 
aaxiaua  accuracy  of  tha  approx i salad  function.  Both  tha  orguaant 
'Xx'  and  tha  ratumad  valua  ‘Sin*  ora  axprassad  in  1750R  double- 
pracision  fixad-point  raprasantotion.  Both  hava  ona  sign  bit,  ona 
intagar  bit,  and  30  fractional  bits.  Tha  actual  valua  ratumad  is 
coaputad  by  tha  function  'SinCos'.  Bacausa  tha  idantity  function, 
Sin(x)  ■  cos<x  -  pi/2),  'SinCos'  can  ba  invokad  to  coaputa  both 
sina  and  cosina  values. 

PRSSED  URRIR8LES:Xx  -  Tha  pi-radian  saasura  for  ahich  sina  is  coaputad. 

Tha  orguaant  is  in  doubla  precision  fixad-point. 
RETURNS:  Sin  -  tha  coaputad  valua  in  doubla  precision  fixad-point 
MOOUUES  CRLL.ED:  SinCos 
AUTHOR:  Copt.  Stavan  fl.  Hotchkiss 

HISTORY:  This  project  aos  undertaken  as  a  thesis  project  for 

partial  fulfil laent  of  raquiraaants  for  an  NS  degree 
in  Infaraatlon  Science  fro»  the  Rir  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  RSO 
Language  Control  Branch,  Wright  Patterson  RFB,0h. 


ITEM  Xx 

END 

R 

1,30; 

mm  sin 

Procedure 

PflOC  Sin  RENT 

BEGIN 

<Xx) 

fl  1,30; 

ITEM  Xx 

Sin  -  SinCos<Xx); 

R  1,30; 

ORTH:  19  JULY  1983 

VERSION:  1.0 

NRNE:  SinCo* 

MOOULE  NUMBER:  1.1 
DESCRIPTION: 

This  function  is  called  by  either  ‘Sin*  or  'Cos'.  The  algor i the 
used  is  designed  around  a  polynoeial  approximation  of  the  sine 
function.  The  coefficients  were  derived  through  a  Chebyshev 
economization  of  the  poeer  series.  Since  'Sin*  and  ’Cos’  use 
pi-radians  as  an  argument,  the  economized  coefficients  eere 
multiplied  by  an  appropriate  poser  of  pi.  The  first  step  of  the 
algorithm  is  to  reduce  the  argument  to  the  interval  for  shich 
the  polynomial  is  valid  <-.5  pi,  .5  pi >.  The  next  step  is  to 
determine  if  the  argument  mould  cause  cn  under  floe,  if  it  does 
the  value  of  the  argument  is  the  value  returned  <sin(x) — >x  for 
small  x).  In  fixed-point  arithmetic,  underflow  causes  doesn't 
cause  problems,  but  checking  for  undarfloe  prevents  unnecessary 
computations.  The  computed  result  of  the  polynomial  reflects  the 
use  of  Homers  rule.  This  function  and  its  argument  are  double 
precision  fixed-point  values. 

PR88ED  VARIABLES:  Xx  -  dochle  precision  fixed-point  representation  of 

an  angle  expressed  in  pi-radians  (radians 
divided  by  pi  equals  pi-radians>.  Legal  range 
_ of  values  include  -1.0  to  1.0 

RETURNS:  R  1730R  double  precison  fixed-point  representation 

of  'Sin'  or  'Cos' 

MODULES  CALLED:  None 

AUTHOR:  Capt.  Steven  A.  Hotchkiss 

HISTORY:  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  MS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  RSO 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


PROC  SinCos  RENT  (Xx)  A 

1,30; 

BEGIN 

CONSTANT 

ITEM 

Eps 

fl  1,30 

m 

0.0003348317; 

CONSTANT 

ITEM 

PiFixed 

R  2,29 

m 

3.1415926530; 

CONSTANT 

ITEM 

R1 

F  39 

m 

-1.644934061140; 

CONSTANT 

ITEM 

R2 

F  39 

m 

♦0.8117421707751 

CONSTANT 

ITEM 

R3 

F  39 

m 

-0. 1907476226769 

CONSTANT 

ITEM 

R4 

F  39 

m 

♦0.0261162053162 

CONSTANT 

ITEM 

A3 

F  39 

m 

-2.23522403740561 

CONSTANT 

ITEM 

PiFloat 

F  39 

m 

3. 141592653589; 

CONSTANT 

ITEM 

One 

F  39 

M 

+  1.0: 

a. 

3 


i 


DRTE:  18  JULY  1985 

UERSION:  1.0 

NRNE:  CosF 

NOCULE  NUHBER:  1.0 
DESCRIPTION: 

This  function  is  invoked  to  coaputa  the  floating  point  value  for 
the  cosine  of  an  angle  that  has  been  expressed  in  radian  eeasure. 
Both  the  a-gueent  'Xx‘  and  the  returned  value  'Cos'  are  in  1750A 
extended  precision  floating  point  representation.  The  actual 
value  returned  is  coeputed  by  the  function  'SlnCosF'.  Because  of 
the  identity  function,  sin<x)  ■  cos(x  -  pi/2),  'SinCosF*  can 
be  invoked  to  coepute  both  sine  and  cosine  values. 

PASSED  UARIABLES:Xx  -  The  radian  eeasure  for  ehich  sine  is  coeputed. 

The  argueent  is  in  extended  precision  float. 
RETURNS:  CosF  -  the  coeputed  value  in  extended  floating  point. 

NODULES  CALLED:  SinCosF 

AUTHOR:  Capt.  Steven  R.  Hotchkiss 

HISTORY:  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfil  leant  of  requireoents  for  an  MS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


PROC  SinCosF  RENT  <Xx>  F  3 
BEGIN 

I  TEN  Xx  F 

END 


Cosf  Procedure 


PROC  CosF 


RENT  <Xx) 


BEGIN 


CONSTANT  I TEN  Vmax 
CONSTANT  I TEN  PiDivByTeo 
CONSTANT  I  TEN  Zero 


+0. 232935020988 1E+7 
1.57079032679489; 
0.0; 


Yy  -  Xx; 

IF  <Yy  >  Yaax)  OR  <Yy  <  -Yeax); 

CosF  ■  Zero; 

ELSE 

IF  Vy  <  Zero; 

Yy  -  -  Yy; 


Yy  ■  PiOivByTeo  +  Yy; 


>  n  inim  i 

DOTE:  18  JULV  1985 

UERSION:  1.0 

NAME:  SinF 

nOOULE  NUMBER:  1.0 
DESCRIPTION: 

This  function  is  invoked  to  coeput*  th*  floating  point  value  for 
th*  sin*  of  an  angl*  that  has  been  expressed  in  radian  **asur*. 
Both  th*  arguMnt  'Xx'  and  th*  r*tum*d  value  'Sin1  or*  in  1750R 
extended  precision  floating  point  r*pr*s*ntation.  Th*  actual 
value  r*tum*d  is  co*put*d  by  th*  function  ‘SinCosF'.  B*caus*  of 
th*  Identity  function,  Sln<x>  -  cos(x  -  pi/2),  'SinCosF'  can 
b*  invoked  to  coaput*  both  sin*  and  cosin*  values. 

PASSED  UARIABLES:Xx  -  Th*  radian  **asur>  for  *hich  sin*  Is  co*put*d. 

Th*  arguaant  is  in  extended  precision  float. 
RETURNS:  SinF  -  th*  co*put*d  value  in  extended  floating-point. 

NODULES  CALLED:  SinCosF 

AUTHOR:  Capt.  St*v*n  R.  Hotchkiss 

HISTORY:  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfilleent  of  requirements  for  an  NS  degree 
in  Inforeation  Science  fro#  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASO 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


START 


REF  PROC  SinCosF 
BEGIN 

RENT  <Xx>  F 

39; 

ITEM  Xx 
END 

F 

39; 

"i|H|nMD| 

*********** 

Sinf  Procedure 

OEF  PROC  SinF 

RENT  «x> 

F  39; 

BEGIN 

ITEM  Xx 

F  39; 

SinF  »  SinCosF(Xx); 


DATE:  19  JULY  1985 

UERSION:  1.0 

NAME:  SinCosF 

MODULE  NUMBER:  1.1 
OESCRIPTIOH: 

This  function  is  called  by  either  'SinF'  or  ‘CosF'.  The  algorithm* 
used  is  a  polynomial  approximation  of  the  sine  function.  The 
coeffmcients  mere  determined  through  a  Chebyshev  Economization 
of  the  power  series,  and  mere  carried  out  to  I750R  machine 
precision.  For  efficient  computation,  the  polynomial  mas 
computed  using  Homer's  Rule.  This  function  returns  an  extended 
precision  floating-point  value. 

PASSED  UARIABLES:  Job  -  a  flag  indicating  whether  to  compute  either 

the  sine  or  cosine  of  an  angle  expressed  in 
radians 

AbsX  -  the  absolute  value  of  the  angle  under 
consideration 

Vy  -  originally  set  to  AbsX,  but  modified 

by  the  algorithm  for  range  reduction  of 
the  original  angle. 

RETURNS:  An  Extended  precision  floating-point  approximation  to 

a  user  requested  call  to  either  'SinF'  or  'CosF' 

MOOULES  CALLED:  Hone 

AUTHOR:  Copt.  Steven  A.  Hotchkiss 

HISTORY:  This  project  mas  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  US  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASO 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


PROC  SinCosF  RENT  <Xx>  F 

39; 

BEGIN 

CONSTANT 

ITEM 

R1 

F  39 

-  - 1 . 5449340668480 ; 

CONSTANT 

ITEM 

A2 

F  39 

■  +0.81 17424252778; 

CONSTANT 

ITEM 

R3 

F  39 

«  -0. 1907518239486; 

CONSTANT 

ITEM 

R4 

F  39 

■  +2 . 6 147845 1583 10E-2; 

CONSTANT 

ITEM 

R5 

F  39 

-  -2 . 3460587938600E-3 ; 

CONSTANT 

ITEM 

R6 

F  39 

■  +1.4832528223590E-4; 

CONSTANT 

ITEM 

R7 

F  39 

■  -6. 7236447557 180E-6; 

CONSTANT 

ITEM 

Eps 

F  39 

-  +0. 1348699 1523486  IE-5 

CONSTANT 

ITEM 

Pi 

F  39 

■  3.1415926535898; 

CONSTANT 

ITEM 

One 

F  39 

-  1.0; 

CONSTANT 

ITEM 

Zero 

F  39 

-  0.0; 

CONSTANT 

ITEM 

OneHalf 

F  39 

*  0.5; 

CONSTANT 

ITEM 

OneOivByPi 

F  39 

■  +0.31 8309886 18379E0; 

ITEM  Sign 

F 

39; 

ITEM  Gg 

F 

39; 

ITEM  Xx 

F 

39; 

ITEM  P  i  Pad  i  arts 

F 

39; 

ITEM  Rasul  t  F 

TRBLE  Overlays  (0)  U  6; 


BEGIN 

ITEM 

Ff 

F 

30 

POSC0,0) 

ITEM 

Fexp 

S 

7 

P0SC8, 1 ) 

ITEM 

Fbits 

B 

48 

POSCO, 0) 

ITEM 

Yy 

F 

30 

POSCO, 3> 

ITEM 

Y  «Xp 

S 

7 

P0S<8,4) 

ITEM 

Vbits 

B 

48 

POSCO, 3 > 

END 

Ff<0)  -  Xx  *  OneOivByPi ; 
Sign  *  One; 

IF  Ff<0)  <  Zero; 

BEGIN 

FfCO)  *  -FfCO); 

Sign  *  -  Sign; 

END 


Vy<0)  ■  Zero; 

IF  FexpCO)  >-  1; 

BEGIN 

IF  B I  TCFbitsCO),  FexpCO), 1)  ■  IB* 1  * ; 

Sign  ■  -Sign; 

VexpCO)  ■  FexpCO); 

BIT<Vbi  tsCO ),0, FexpCO >+1)  «  BITCFbl  tsC0),0,FexpC0H1 ); 

END 

PiRadians  »  FfCO)  -  VyCO); 

IF  PiRadians  >  OneHalf; 

PiRadians  ■  One  -  PiRadians; 

IF  PiRadians  <  Eps; 

Result  -  PiRadians  *  Pi; 

ELSE 

BEGIN 

Gg  ■  PiRadians  ♦  PiRadians; 

Result  ■  CCCCCCCC  R7  *  Gg  +  R6)  *  Gg  +  R5)  *  Gg  +  R4)  *  Gg  +  R3)  *  Gg 
+  R2)  *  Gg  +  R1)  *  Gg>  +  One)  *  PiRadians  *  Pi; 

END 

SinCosF  ■  Result  *  Sign; 

RETURN; 


*  DflTE:  6  August  1985  * 

*  UERSION:  1.0  * 

*  NAME:  Cot  * 

*  MODULE  NUMBER:  1.0  * 

*  DESCRIPTION:  * 

*  Cot  is  invoked  to  coaputa  the  cotangent  of  a  usar  passed  angle.  The* 

*  angle  is  expressed  in  pi-radians:  pi-radians  — >  radians/pi .  * 

*  The  range  of  legal  values  for  the  angle  are  <-1.0,1.0»  Both  the  * 

*  arguaent  and  result  are  expressed  as  fixed  point  results.  The  * 

*  value  returned  is  coaputed  by  the  function  TanCot.  * 

*  PASSED  VARIABLES:  Xx  -  the  angle  in  pi-radians  * 

*  RETURNS:  a  value  for  the  cotangent  of  Xx.  * 

*  plus  and  einus  infinity  are  considered  as  -4095  and  4095  * 

*  MODULES  CALLED:  TanCot  * 

*  AUTHOR:  Capt.  Steven  A.  Hotchkiss  * 

*  HISTORV:  This  project  was  undertaken  as  a  thesis  project  for  * 

*  partial  fulfillaent  of  requireaents  for  an  NS  degree  * 

*  in  Information  Science  froe  the  Air  Force  Institute  ♦ 

*  of  Technology.  Sponsoring  organization  is  the  ASD  * 

*  Language  Control  Branch,  Uright  Patterson  AFB,0h.  * 


PROC 

TanCot  RENT<Arg, 

Job) 

A  12,18; 

BEGIN 

ITEM 

Arg  A  1,30; 

ITEM 

Job  B; 

END 

V 

"X 

— 

* 

Cot 

Procedure 

PROC 

o 

o 

c+ 

a 

? 

V 

A  12,18; 

BEGIN 

DEFINE  Tangent  ' 

IB’O' 

1  M  » 

» 

DEFINE  Cotangent  * 

«•  r 

1  M  , 

* 

ITEM 

Xx 

A  1 

,30; 

ITEM 

Job 

B; 

slob  * 

Cotangent; 

Cot  ■ 

Tar>Cot<Xx,  Job); 

falij 


OATE:  6  August  1069 

VERSION:  1.0 

NRNE:  Ten 

NOODLE  NUMBER:  1.0 
DESCRIPTION: 

Tan  is  invoked  to  cosputa  ths  tangont  of  a  ussr  passsd  angle.  Tbs 
angle  is  expressed  in  pi-radians:  pi-radians  — >  radians/pi . 

The  range  of  legal  values  for  the  angle  are  <-1.0, 1.0>>  Both  the 
argusent  and  result  are  expressed  as  fixed  point  results.  The 
value  returned  is  coeputed  by  the  function  TanCot. 

PASSED  VARIABLES:  Xx  -  the  angle  In  pi -rad I  arts 
RETURNS:  a  value  for  the  tangent  of  Xx. 

plus  and  sinus  infinity  are  considered  as  -4009  and  4009 
MODULES  CALLED:  TanCot 
AUTHOR:  Capt.  Stevai  A.  Hotchkiss 

HISTORV:  This  project  sas  undertaken  as  a  thesis  project  for 

partial  fulfil  leant  of  requi resents  for  an  NS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  AS D 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


PAOC  TanCot  RENT<Arg,  Job)  A  12,18; 
BEGIN 

ITEH  Arg  A  1,30; 

ITEM  Job  B; 

END 


Tan  Procedure 


PAOC  Tan  HEHT<Xx)  A  12, 18; 
BEGIN 

DEFINE  Tangent  "1B'0'“; 
DEFINE  Cotangent  “IB'  I'"; 


ITEH  Xx 
ITEH  Job 


A  1,30; 

B; 


Job  *  Tangent; 

Tan  ■  TanCot<Xx,  Job); 


RETURN: 


I 


DATE:  6  August  1969 

VERSION:  1.0 

MANE:  TanCot 

MODULE  NUNBER:  1.1 
DESCRIPTION: 

This  function  is  invoked  by  either  ‘Tan*  or  'Cot*  to  coepute  the 
appropriate  value  for  a  given  angle.  The  angle  'Arg'  is  given  in 
pi-radians:  pi-radian  — >  radians/pi.  Since  tan<x>  -  1/cot(x> 

both  functions  can  call  this  routine  to  coepute  the  desire  values. 
Since  Tangent  and  Cotangent  approach  plus  or  sinus  infinity  at 
certain  angles,  a  check  Is  required  to  prevent  degradation  due  to 
underfloe  or  overfloe,  plus  and  sinus  infinity  for  this  procedure 
are  considered  to  be  4099  and  -4099  respectively.  The  coefficient 
used  in  the  teo  polt/toelals  eere  determined  through  the  PACE 
method.  The  original  polynomial  mas  determined  through  a 
Chebysheff  econimizat ion  of  the  poeer  series  for  tangent. 

PASSED  VARIABLES:  Arg  -  an  angle  expressed  in  pi-radians 

Job  -  tel  Is  ehether  to  compute  tangent  or  cotangent 
RETURNS:  Either  the  tangent  or  cotangent  of  a  given  angle. 

there  are  12  integer  bits  and  18  fraction  bits  for  this 
function. 

MOOULES  CALLED:  None 

AUTHOR:  Copt.  Steven  A.  Hotchkiss 

H I STORY :  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  NS  de9*ee 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


START 


OEF  PfiOC  TanCot  AENKArg,  Job)  A  12,18; 
BEGIN 


DEFINE  Tangent 
DEFINE  Cotangent 
DEFINE  Pg 
DEFINE  Qg 


"IB'O”; 

“IBT"; 

■Polynomial  <0>“; 
■Polynomial <1>"; 


ITEM  Job 
ITEM  Arg 
ITEM  Xx 
ITEM  Xfloat 
(TEN  Gg 
ITEM  Sign 


ITEM  One  STATIC  F  39  ■  1.0; 


CONSTANT  ITEM  PiFloat  F  39-3.141992693990 

CONSTANT  ITEM  PIFixed  A  2,29-3.1419926936; 

CONSTANT  ITEM  Eps  A  1,30  -  0.00007773138; 

CONSTANT  ITEM  UpperLimit  A  12,18  -  4099.0; 


ITEM  PO  STATIC  F  30; 

CONSTANT  ITEM  PI  F  30  — 1 .000003361855; 

CONSTANT  ITEM  P2  F  30  -  0. 104720845707; 

OUERLAV  PiFloat:  PO; 

ITEM  QO  STATIC  F  30; 

CONSTANT  ITEN  Q1  F  30  —4 . 38806 1470 150; 

CONSTANT  ITEN  Q2  F  30  -  1.555053508405; 

OUERLAV  On*:  QO; 


TABLE  frray  <0:1 >; 

BEGIN 

ITEN  Polynomial  F  30; 

END 

ITEN  Numerator  U  8; 

ITEN  Denominator  U  8; 

(TEN  Num  B  8; 

I  TEN!  Den  8  8; 

OUERLAV  Numerator  :  Num; 

OUERLAV  Denominator  :  Den; 

Sign  *  One; 

Xx  -  Arg; 

IF  Job  ■  Cotangent; 

BEGIN 

Numerator  ■  1; 

Denominator  *  0; 

END 
FI  SF 
BEGIN 

Numerator  ■  0; 

Denominator  *  1; 

END 

IF  Xx  <  0.0; 

BEGIN 

Xx  ■  -Xx; 

Sign  ■  -Sign; 

END 

IF  Xx> 1 .0; 

Xx  ■  Xx  -  1.0; 

IF  Xx  >  .5; 

BEGIN 

Xx  -  1.0  -  Xx; 

Sign  »  -Sign; 

END 

IF  Xx  >  .25; 

BEGIN 

Xx  *  .5  -  Xx; 

BIT<Num,7,1>'»  NOT  BIT<Num,7, 1>; 
BIT<0*n,7, 1)  -  NOT  BIT<D*n,7, 1>; 
END; 


Xfloat  ■  <*  F  39  *>  <Xx>; 

IF  Xx  <  Ep»; 

BEGIN 

Pg  -  Xfloat  *  PIFIoat; 

Qg  •  &*; 

END 

ELSE 

BEGIN 

Gg  »  Xfloat  *  Xfloat; 

Pg  ■  «P2  *  Gg  ♦  P1>  *  PIFIoat  ♦  P0>  *  Xfloat; 

Qg  -  <  02  *  Gg  +  Q1>  *  Gg  +  QO; 

END 

TanCot  -  <*  fl  12, 18  *)  (Sign  * 

PoIgnoaiaKNu— rotor  )/Polgnoaiol<Dono»ina  tor)  > 

RETURN; 


DATE:  19  JULV  1985 

UERSION:  1.0 

NAME:  CotF 

tIOOULE  NUMBER:  1.0 
DESCRIPTION: 

This  function  is  invoked  to  coapute  the  cotangent  of  an  angle 
expressed  in  radian  eeasure.  Because  of  the  identity  function 
Tan<X>  *  1/Cot<X>,  the  teo  functions  'Tan*  and  'Cot'  can  both 
invoke  the  function  'TanCot'  to  coopute  their  respective  values. 
Both  the  ar queen t  and  the  result  of  this  function  are  expressed 
in  17S0R  extended  precision  floating-point  representation. 
PASSED  URRIR6LES:  Xx  -  the  angle  of  interest,  expressed  in  extended 

floating-point  representation.  The  angle  oust 
lie  beteeen  (  -2329349.90332,  2329349.90332  > 
RETURNS:  Cotangent  of  the  angle  Xx  <-1.0  to  1.0>  in  extended 

floating-point  representation 
ntXXJLES  CALLED:  TanCotF 
AUTHOR:  Capt.  Steven  A.  Hotchkiss 

HISTORY:  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfil leent  of  requireeents  for  an  NS  degree 
in  Inforeation  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASO 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


PR0C  TanCotF  RENT  <Arg, 

Job)  F  39; 

BEGIN 

ITEM 

Arg 

F 

39; 

ITEM 

Job 

B; 

END 

Cotf  Procedure 

PA0C  CotF  RENT  <Xx)  F 

39; 

BEGIN 

DEFINE 

Tangent 

"  IB 1 0  * w ; 

DEFINE 

Cotangent 

"IB* I*"; 

ITEM 

Xx 

F 

39; 

ITEM 

Job 

B; 

Job  *  Cotangent; 

CotF  ■ 

TanCotF<Xx, 

Job); 

DOTE:  19  JULV  1989 

UERSION:  1.0 

NAME:  TanF 

MOOULE  NUMBER:  1.0 
DESCRIPTION: 

This  function  is  invoked  to  coepute  the  tangent  of  on  angle 
expressed  in  radian  Masura.  Because  of  the  identity  function 
Tan<X)  ■  1/Cot(X),  the  too  functions  'Tan*  and  'Cot'  can  both 
invoke  the  function  'TanCot'  to  coepute  their  respective  values. 
Both  the  ergueent  end  the  result  of  this  function  are  expressed 
in  17S0R  extended  precision  floating-point  representation. 
PASSED  UARIA6LES:  Xx  -  the  angle  of  interest,  expressed  in  extended 

floating-point  representation.  The  angle  oust 
lie  beteeen  (  -2329349.90332,  2329349.90332  > 
RETURNS:  Tangent  of  the  angle  Xx  (-1.0  to  1.0)  in  extended 

floating-point  representation 
NOOULES  CALLED:  TanCotF 
AUTHOR:  Copt.  Steven  A.  Hotchkiss 

H I STORY :  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfil  leant  of  requireeants  for  an  NS  degree 
in  Inforeation  Science  froe  the  Air  Farce  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Urirfit  Patterson  AFB,0h. 


START 

REF  PA0C  TanCotF  RENT  (flrg.  Job)  F  39; 
BEGIN 

ITEM  Arg  F  39; 

ITEM  Job  B; 

END 


Tanf  Procedure 


0EF  PA0C  TanF  RENT  (Xx)  F  39; 

BEGIN 

DEFINE  Tangwit  "IB'D"*; 

DEFINE  Cotangent  "IB' 1  * ", 

ITEM  Xx  F  39; 

ITEM  Job  B; 

Job  ■  Tangent; 

TanF  ■  TanCotF (Xx,  Job); 

RETURN; 

END 

TERM 


OflTE:  19  JULV  1989 

UERSION:  1.0 

NAME:  TanCotF 

HOOULE  NUMBER:  1.1 
DESCRIPTION: 

This  function  is  called  by  either  'Tan*  or  'Cot'  to  coepute  the 
tangent  or  cotangent  of  their  respective  angles.  The  result  is 
cooputed  by  using  a  Rational  approxieation  t  P<X)/Q<X>  I.  The 
coefficients  used  eere  found  by  a  Pade  approxieation  of  the 
Chebyshev  econoeizatian  of  the  poeer  series  for  tangent  functions. 
The  result  of  this  function  Is  In  1730R  extended  precision 
floating-point  representation. 

PRSSED  VARIABLES:  Xx  -  the  angle  under  consideration,  expressed  in 

radians.  Xx  eust  lie  beteeen 
(  -2329349.90332,  2329349.90332) 

Vy  -  absolute  value  of  Xx,  used  to  prevent  overflow 
I  flag  -  designates  whether  to  coepute  Tan  or  Cot 
RETURNS:  1790R  extended  precision  float  for  Tan  or  Cot 

MODULES  CflLLED:  None 
AUTHOR:  Capt.  Steven  R.  Hotchkiss 

HISTORV:  This  project  was  undertaken  as  a  thesis  project  for 

partial  fulfil  leant  of  requi regents  for  an  MS  degree 
in  Information  Science  from  the  Rir  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


39 


(TEH  Nuaarator  U  8; 

ITEM  Oanoainator  U  8; 

ITEM  Mum  8  8; 

ITEM  Dan  8  8; 

OUEBLAV  NuMrator  :  Nub; 
OUERLAV  Oanoainator  :  Dan; 


ITEM  Gg 

F 

39; 

ITEM  Sign 

F 

30; 

CONSTANT 

ITEM 

PO 

F 

30 

■+1.0; 

CONSTANT 

ITEM 

PI 

F 

30 

—  1.2661071041410; 

CONSTANT 

ITEM 

P2 

F 

30 

■+0. 27332 10453881; 

CONSTANT 

ITEM 

P3 

F 

30 

—7. 1046857833630E-3; 

CONSTANT 

ITEM 

00 

F 

30 

—►1.0; 

CONSTANT 

ITEM 

Q1 

F 

30 

—4 . 5550752378380; 

CONSTANT 

ITEM 

02 

F 

30 

-+2. 2740008037200; 

CONSTANT 

ITEM 

03 

F 

30 

—0.2003006071354; 

CONSTANT 

ITEM 

Vaax 

F 

30 

-+0 . 23203400033203E+7 

CONSTANT 

ITEM 

Zaro 

F 

30 

-  0.00; 

CONSTANT 

ITEM 

Ona 

F 

30 

»  1.00; 

CONSTANT 

ITEM 

OnaHalf 

F 

30 

-  0.50; 

CONSTANT 

ITEM 

OnaFourth 

F 

30 

■  0.25; 

CONSTANT 

ITEM 

TaoOiwByPi 

F 

30 

»+0. 6366 107723675; 

CONSTANT 

ITEM 

PIFIoat 

F 

30 

*+3. 1415026535808; 

CONSTANT 

ITEM 

Eps 

F 

30 

■+0. 1348600 152348E-5; 

Sign  ■  Ona; 

Yy<0)  »  Zaro; 

Xx<0)  •  Arg  *  TaoDivByPi; 

IF  Job  ■  Cotangant; 

BEDIM 

Nuaarator  »  1 ; 
Oanoainator  *  0; 

END 

ELSE 

BEGIN 

Nuaarator  *  0; 
Oanoainator  ■  1; 

END 

IF  Xx(Q)  <  Zaro; 

BEGIN 

Xx<0>  ■  -Xx<0); 

Sign  ■  -Sign; 

END 


IF  «oxp<0>  >  0)  AND  <BIT<Xbi  ts<0>,Xoxp<0>,  1  >  «  IB’ V  >; 

BEGIN 

Sign  •  -Sign; 

BIT<Nua,7, 1 )  »  NOT  BIT<Nua,7, 1); 

BIT<0an,7, 1 >  ■  NOT  BIT<0an,7, 1 ); 

END; 


IF  Xaxp<0>  >  1; 

BEGIN 

Ytxp(O)  ■  Xaxp(O); 

BIKYbi  t*<0),0,Xaxp<0>+1 )  ■  BITOfei  t*<0),0,Xaxp<0>+1  >; 

END; 

Xx<0>  -  Xx<0>  -  Yg<0); 

IF  Xx<0>  <>  Zaro; 

Xaxp(O)  »  Xaxp(O)  -(*87  *)(  1  >; 

IF  Xx(O)  >  OrwFourth; 

BEGIN 

Xx<0>  ■  OnaHal f-Xx<0); 

BIT<Nua,7, 1)  ■  NOT  BIT<Nua,7, 1); 

BlT<Dan,7, 1>  -  NOT  BIT<D«i,7, 1); 

END 

IF  Xx(O)  <  Ep»; 

BEGIN 

Pg  «  Xx<0>  *  PiFloat; 

Qg  «  Ona; 

END 

ELSE 

BEGIN 

Gg  ■  Xx(O)  *  Xx(O); 

Pg  -  «<P3  *  Gg  ♦  P2>  *  Gg  ♦  P1>  *  Gg  *  PiFloat  ♦  PO>  *  Xx<0> 
Qg  ■  «  03  *  Gg  +  02)  *  Gg  +  Q1>  *  Gg  +  QO; 

END 

TanCotF  »  Sign  *  Polgnoaial <Nuaarator)/Polynoaial (Danoainator); 
RETURN; 


DflTE: 

UERSION: 

Ill  L . 

MODULE  NUNBER: 
DESCRIPTION: 


21  August  1985 
1.0 
RCos 
1.0 


This  program  is  used  to  call  ASInCos.  FtSinCo*  is  the 
routine  that  actually  computes  the  inverse  cosine  for 
argument  Xx.  They  are  coded  this  eay  because  both  A8in 
RCos  can  use  the  same  routine  for  their  coeputations. 
This  is  because  of  the  identity  function 
ACos(X)  >  pi/2  -  RSin(x) 

PASSED  UARIABLES:  Xx  -  the  cosine  for  ehich  an  angle  is  to  be 

coeputed . 

RETURNS.  RCos  -  expressed  in  fixed-point  pi -radians 

NOOULES  CALLED:  RSinCos 

AUTHOR:  Capt.  Steven  A.  Hotchkiss 

HISTORV:  This  project  wot  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  NS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASO 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


PROC  ASinCas  RENTCArg,  Job)  A  1,30; 
BEGIN 

ITEN  Arg  A  1,30; 

I  TEN  Job  U  8; 

END 


ACos  Procedure 


PROC  ACos 
BEGIN 


RENT<Xx) 


DEFINE  ArcSine 
DEFINE  ArcCosine 


ITEN  Xx  A  1,30; 
ITEN  Job  U  8; 


Jab  ■  ArcCosine; 

ACos  ■  ASInCos<Xx,  Job); 


DOTE: 

VERSION: 

NRME: 

NOOULE  NUMBER: 
DESCRIPTION: 


21  August  1983 
1.0 
RSin 
1.0 


This  program  is  used  to  call  RSinCos.  RSinCos  is  the 
routine  that  actually  computes  the  inverse  sin*  for 
argumnt  Xx.  They  are  coded  this  way  because  both  RSin 
RCos  can  use  the  saee  routine  for  their  computations. 
This  is  because  of  the  identity  function 
ACos(X)  ■  pi/2  -  RSIn<x> 

PASSED  VARIABLES:  Xx  -  the  fixed-point  inverse  for  ehich  an  angle  is 

is  to  be  coeputed 

RETURNS:  RSin  -  the  angle  In  pi-radians. 

riOOULES  CALLED:  RSinCos 

AUTHOR:  Copt.  Steven  R.  Hotchkiss 

H I STORY :  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  US  degree 
in  Information  Science  from  the  Rir  Force  institute 
of  Technology.  Sponsoring  organization  is  the  RSD 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


START 


REF  PROC  RSinCos  RENTCflrg,  Job)  R  1,30; 
BEGIN 

ITEM  flrg  A  1,30; 

ITEM  Job  U  8; 

END 


RSin  Procedure 


DEF  PROC  RSin  RENT(Xx) 
BEGIN 

DEFINE  ArcS in*  "0"; 
DEFINE  RrcCos in*  "1"; 

ITEM  Xx  fl  1,30; 

ITEM  Job  U  8; 


Job  -  RrcS in*, 

RSin  ■  RSinCosCXx,  Job); 


RETURN; 


DATE:  21  August  1983 

UERSION:  1.0 

NATE:  ASinCos 

MOOULE  NUMBER:  1.1 
DESCRIPTION: 

This  program  is  eritten  to  return  the  correct  angle  for  a  user 
passed  argument.  The  argument  represents  the  Sine  or  Cosine  of 
the  angle  to  be  returned.  The  returned  angle  is  in  pi-radian 
eeasure  (pi-radian  ->  radians/pi).  This  routine  can  be  called 
by  either  ASin  or  ACos  because  of  the  identity  function 
ACos(x>  ■  pi/2  -  ASin<x>.  The  coefficients  sere  determined  by 
of  the  method  of  Chebysheu  expans  ion  as  described  in  'A  First 
Course  in  Numerical  Analysis'  by  Anthony  Ralston. 

PASSED  UARIABLES:  Arg  -  the  sine  or  cosine  of  the  angle  to  be 

determined.  Uariable  is  in  fixed-point 
Job  -  tells  ehether  to  compute  for  ACos  or  ASin 
RETURNS:  The  angle  representation  for  the  argument.  The  angle 

is  in  pi-radians.  Leagal  values  fall  in  the  range 

(-1.0,  1.0) 


PASSED  UARIABLES: 


RETURNS: 


NODULES  CALLED: 
AUTHOR: 

HISTORY: 


Capt.  Steven  A.  Hotchkiss 

This  project  eas  undertaken  as  a  thesis  project  for 
partial  fulfillment  of  requirements  for  an  NS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASO 
Language  Control  Brcrtch,  Uright  Patterson  RFB.Oh. 


PROC 

Sqrt 

RENT(Xx) 

F 

39; 

BEGIN 

ITEM 

Xx 

F  39; 

END 

nsmvos  rroceaure 

PROC 

ASinCos 

RENT (Arg,  Job) 

A 

1,30; 

BEGIN 

DEFINE  ArcSIne 

“0"; 

DEFINE  ArcCosine 

T; 

DEFINE  Positive 

"0"; 

DEFINE  Negative 

T; 

ITEM 

Job 

U  8; 

ITEM 

li 

U  8; 

ITEM 

Arg 

A  1,30; 

ITEM 

Result 

A  1,30; 

ITEM 

Pg 

F  39; 

ITEM 

Qg 

F  39; 

TABLE  (XMrloysCO)  U  6; 

BEGIN 

ITEM  Vy  F  39  POS<0,0>; 

ITEM  V«xp  S  7F>OS<8, 1); 
ITEM  Gg  F  39  P0S<0,3>; 

ITEM  G*xp  S  7  P0S<8,4>; 


CONSTANT  TABLE  Constants  (0:1)  U  4  ■  0.0,  1.0,  0.5,  0.5; 
BEGIN 

ITEM  Aa  A  1,30  F>OS<0,0>; 

ITEM  Bb  A  1,30  P0S(0,2>; 

END 


CONSTANT 

ITEM 

Eps 

F 

39 

a 

9 . 587379924290E-5; 

CONSTANT 

ITEM 

On*Ov*rPI 

F 

39 

St 

0.31830988018379; 

CONSTANT 

ITEM 

On* 

F 

39 

a 

1.0; 

CONSTANT 

ITEM 

OrwHalf 

F 

39 

a 

0.5; 

CONSTANT 

ITEM 

Orwlnt 

S 

7 

a 

i; 

CONSTANT 

ITEM 

Xnax 

A 

1,30 

a 

1.999999999; 

CONSTANT 

ITEM 

PI 

F 

39 

a 

-0 . 275 10555290590E+ 1 ; 

CONSTANT 

ITEM 

P2 

F 

39 

m 

+0 . 29058702374859E+ 1 ; 

CONSTANT 

ITEM 

P3 

F 

39 

a 

-0. 59450 144193240E+0; 

CONSTANT 

ITEM 

00 

F 

39 

m 

-0.  16509933202424E+2; 

CONSTANT 

ITEM 

01 

F 

39 

a 

♦0. 24804728909 104E+2; 

CONSTANT 

ITEM 

02 

F 

39 

a 

-0. 10333807072 113E+2; 

CONSTANT 

ITEM 

03 

F 

39 

a 

+0. 10000000000000E+1; 

Yy<0>  *  <*  F  39  *)<  ABS(Arg)  ); 

* 

IF  Vy<0>  <■  OrwHalf; 

BEGIN 
I  i  »  Job; 

IF  Vy<0>  <  Ep*; 

BEGIN 

A*sul t  ■  <*  A  1,30  *>  <Vy<0>  *  0n*0v*rPi  >; 

GOTO  LI; 

END 

Gg<0>  -  Vy<0)  *  Vy<0>; 

END 

ELSE 

BEGIN 

I  i  *  1  -  Job; 

IF  Vy<0)  >  On*; 

BEGIN 

ASinCos  ■  Xhx; 

ABORT; 

END 

Gg  (0)  ■  On*  -  Vy<0>; 

G*xp<0)  ■  G*xp<0)  -  On*lnt;  “  Gg  ■  Gg/2  or  Gc  •  Gg  *  2 
Vy  (0)  ■  -Sqrt<Gg<0)>; 

V«xp<0>  ■  V*xp<0>  ♦  On*lnt;  "  Vy  ■  Vy*2  or  Vy  ■  Vy  *  2 


Pg  ■  <<P3  *  Gg<0>  ♦  P2>  *  0g<0>  *  PI)  *  0g<0); 

Qg  ■  <<  Gg<0)  +  02)  *  Gg<0>  +  Q1>  *  Gg<0)  +  QO; 

R**uit  »  <*  fi  1,  »  *)  <<Vy<0)  ♦  Vy<0>  *  Pg  /  Qg)  *  OneOverPi  ) 

IF  Job  ■  RrcSine; 

BEGIN 

Result  ■  <*  A  1,30  *>  <  Result  *  Ratli)  ); 


IF  Rrg  <  0.0; 

Result  ■  -Result; 

END 

H SE  'ELSE  Job  ■  RrcCosine 

IF  Rrg  <  0.0; 

Result  ■  <*  A  1,30  *>  <  Bb<l i )  ♦  Result  ); 

ELSE 

Result  ■  <*  fl  1,30  *>  <  flat  I i >  -  Result  >; 


RSinCos  *  Result; 


RETURN: 


'  T7T  77  lJ[*  V  VKmV'*J  LrtMA1  h  r*  «  ■ 


*'7  ^ 


-  •*  k^  i-^W  '. 


DflTE: 

UERSION: 

rift: 

MODULE  NUMBER: 
DESCRIPTION: 


21  August 
1.0 
RCosf 
1.0 


1089 


This  program  is  used  to  call  ASinCosf.  ASinCosf  is  th* 
rout in*  that  actually  co*put*s  th*  invars*  cosin*  for 
argument  Xx.  Thay  ar*  codod  this  aay  bacaus*  both  ASinf 
RCosf  can  us*  th*  same  rout  in*  for  thair  computations. 
This  is  b*caus*  of  th*  identity  function 
RCos(X)  ■  pi/2  -  RSIn(x) 

PASSED  UARIABLES:  Xx  -  th*  cosin*  for  ehich  an  angle  is  to  b* 

co*put*d. 

RCosf  -  *xpr*ss*d  in  floating-point 
ASinCosf 

Capt.  St even  A.  Hotchkiss 

This  project  was  undertaken  as  a  thesis  project  for 
partial  fulfil  leant  of  requirements  for  an  NS  deg r** 
in  Information  Scianc*  from  th*  Air  Fore*  Institute 
of  Technology.  Sponsoring  organization  is  th*  ASD 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


RETURNS: 
N00ULES  CALLED: 
AUTHOR: 

HISTORV: 


REF  PROC  RSinCos 

RENT<Arg,  Job)  F  39; 

BEGIN 

ITEN  Arg 

F 

39; 

ITEM  Job 
END 

U 

8; 

nuoiT  rroewurt 


DEF  PROC  ACos  RENT(Xx)  F  39; 

BEGIN 

DEFINE  ArcSine  ”0“; 

OEFINE  AroCosine  T; 


ITEN  Xx  F  39; 
ITEM  Job  U  8; 


Job  *  ArcCosine; 

RCosf  «  ASinCosf<Xx,  Job); 

RETURN; 


I 


DflTE: 

UERSION: 

rVTC! 

nOOULE  NUMBER: 
DESCRIPTION: 


21  August  1985 
1.0 
RSinf 
1.0 


This  program  is  used  to  call  RSinCosf.  RSinCosf  is  th* 
routins  that  actually  computes  tha  invars*  sin*  for 
arguMnt  Xx.  Th*y  arm  cod*d  this  *ay  because  both  RSinf 
RCosf  can  us*  th*  sqm  routine  for  th*ir  computations. 

This  is  because  of  th*  identity  function 
ACos(X)  >  pi/2  -  ASin(x) 

PRSSED  URRIR6LES:  Xx  -  th*  floating-point  invars*  for  which  an  angle  is* 

is  to  b*  co*put*d 

RETURNS:  RSinf  -  th*  th*  angle  . 

I100ULES  CflLLED:  RSinCosf 

AUTHOR:  Capt.  Steven  R.  Hotchkiss 

HISTORV:  This  project  *as  undartakan  as  a  thesis  project  for 

partial  fulfil  leant  of  requirements  for  an  NS  d*gr** 
in  lnfor*ation  Sci*nc*  fro*  th*  Rir  Fore*  Institute 
of  Technology.  Sponsoring  organization  is  th*  RSD 
Language  Control  Branch.  Wright  Patterson  RFB,0h. 


PROC  RSinCosf 
BEGIN 

RENTCRrg,  Job)  F  39; 

ITEH  flrg  F 

39; 

ITEM  Job  U 

END 

8; 

.a**,******, 

*****  RSinf  Procedure 

PROC  RSin 

RENT<Xx>  F  39; 

BEGIN 

DEFINE  ArcS in* 

"0"; 

DEFINE  ArcCosin* 

T; 

ITEM  Xx  F 

39; 

1  TEN  Job  U 

8; 

Job  ■  RrcSin*; 

RSinf  «  ASinCosfCXx,  Job); 

3 


DflTE:  21  August  1965 

UERSION:  1.0 

NRME:  ASinCosf 

MODULE  NUMBER:  1.1 
DESCRIPTION: 

This  program  is  mrltten  to  return  ths  correct  angle  for  a  user 
passed  argues nt.  The  argueent  represents  the  Sine  or  Cosine  of 
the  angle  to  be  returned.  The  returned  angle  is  in  radian 
Measure.  This  routine  can  be  called 

by  either  RSinf  or  ACosf  because  of  the  identity  function 
ACos(x)  ■  pi/2  -  AS!n<x).  The  coefficients  sere  determined  by 
of  the  eethod  of  Chebyshev  expansion  as  described  in  'A  First 
Course  in  Nuesr ical  Analysis'  by  Anthony  Ralston. 

PASSED  UfflIABLES:  Arg  -  the  sine  or  cosine  of  the  angle  to  be 

determined.  Uariable  is  in  floating-point 
Job  -  tel  Is  another  to  compute  for  ACosf  or  RSinf 
RETURNS:  The  angle  representation  for  the  argument.  The  angle 

is  in  radians.  Leagal  values  fall  in  the  range 
NOOULES  CALLED:  None 
AUTHOR:  Capt.  Steven  A.  Hotchkiss 

HISTORV:  This  project  mas  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  NS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


RETURNS: 


PROC  Sqrt 
COIN 
TEM  Xx 
■NO 


RENT(Xx) 
F  30; 


F  39; 


AsinCos  Procedure 


PROC  ASinCos  RENT  (Arg,  Job)  F  39; 
COIN 


CFINE  ArcSine 
CFINE  ArcCosine 
CFINE  Positive 
CFINE  Negative 


Job 

u 

8 

li 

u 

8 

flrg 

F 

39 

Result 

F 

39 

P9 

F 

39 

Qg 

F 

39 

wjt  r_*  -jt  -ji  -yi 


TABLE  Ov«*lays<0>  U  6; 

BEGIN 

ITEM  Vy  F  39  POSCO, Q>; 

ITEM  Vexp  S  7P0SC8, 1>; 

ITEM  Gg  F  35  POSCO, 3>; 

ITEM  Gexp  S  7  P0SC8,4>; 

END 

CONSTANT  TABLE  Constants  (0:l)US> 

0.0,  1.570756320755,  0.7853561633574,  0.7853561633574; 
BEGIN 


ITEM 

Aa 

F 

35 

POSCO, 0); 

ITEM 

Bb 

F 

35 

POSCO,3); 

END 

CONSTANT 

ITEM 

Eps 

F 

35 

9 

5 . 587375624250E-5; 

CONSTANT 

ITEM 

OneOverPi 

F 

35 

m 

0.31830568618375; 

CONSTANT 

ITEM 

Ons 

F 

35 

m 

1.0; 

CONSTANT 

ITEM 

OnsHalf 

F 

35 

u 

0.5; 

CONSTANT 

ITEM 

Onslnt 

S 

7 

9 

i; 

CONSTANT 

ITEM 

Xmax 

F 

35 

9 

MAXFL0BTC35); 

CONSTANT 

ITEM 

PI 

F 

35 

m 

♦0. 85372 16436677E+1; 

CONSTANT 

ITEM 

P2 

F 

35 

m 

-0. 1 34287075 1 343E+2 ; 

CONSTANT 

ITEM 

P3 

F 

35 

9 

♦0. 55683 1576 1775E+1; 

CONSTANT 

ITEM 

P4 

F 

35 

n 

-0.6540406855Q34E+0; 

CONSTANT 

ITEM 

00 

F 

35 

9 

+0.512232586201 1E+2; 

CONSTANT 

ITEM 

Q1 

F 

35 

9 

-0. 1 0382273 18840E+3; 

CONSTANT 

ITEM 

02 

F 

35 

9 

40. 587 155576538 1E+2; 

CONSTANT 

ITEM 

03 

F 

35 

9 

-0. 1642555755750E42; 

CONSTANT 

ITEM 

04 

F 

35 

9 

40. 10000000000000E4-1; 

ir 

r-v1' 


U 


VyCO)  »  ABSCArg); 


IF  Vy<0)  <-  OnoHalf; 

BEGIN 
II  »  Job; 

IF  Vy(0>  <  Eps; 

BEGIN 

Result  -  Vy<0>; 

GOTO  LI; 

END 

GgC0>  «  VyCO)  *  VyCO); 

END 

ELSE 

BEGIN 

li  •  1  -  Job; 

IF  VyCO)  >  One; 

BEGIN 

ASinCosf  *  Xmax; 

GOTO  L2; 

END 

Gg  <0)  -  One  -  VyCO); 

GexpCO)  ■  Gsxp<0>  -  Onslnt;  “  Og  ■  Og/2  orOg«0g*2**-1 
Vy  <0>  ■  -Sqrt<0g<0)>; 


L- 

H 


115 


Vaxp(O)  «  Vaxp(O)  Onalnt;  “  Vy  -  Vy*2  or  Vy  ■  Vy  *  2  **  t  " 
END 


Pg  a  <<<P4  *  Qg(0>  ♦  P3>  *  Gg<0>  ♦  P2)  *  Gg<0>  +  PI)  *  Gg<0>; 

Qg  ■  «<  Gg<0)  +  Q3)  *  Gg<0>  +  02)  *  Gg<0>  +  01)  *  Gg<0)  ♦  QO; 

Rasul t  -  Vy<0)  +  Vy<0>  *  Pg  /  Qg; 

IF  Job  >  RrcSina; 

BEGIN 

Rasul t  a  Rasul t  ♦  Radi); 

IF  flrg  <  0.0; 

Rasul t  ■  -Rasul t; 

END 

ELSE  "ELSE  Job  -  RrcCosina" 

IF  ft~g  <  0.0; 

Rasul t  ■  Bbdi)  +  Rasult; 

ELSE 

Rasult  *  Radi)  -  Rasult; 


RSinCos  *  Rasult; 


RETURN; 


DATE:  19  JULV  1985 

VERSION:  1.0 

NAME:  MathLib 

MOOULE  NUMBER:  1.0 
DESCRIPTION: 

This  compool  Is  required  by  any  JOUIflL  program  that  needs  to 
reference  any  of  the  oath  functions  eritten  for  f I oa ting-point 
or  fixed-point  computations 
PASSED  UARIABLES:  N/fi 
RETURNS:  N/R 

MODULES  CALLED:  N/R 

AUTHOR:  Capt.  St even  A.  Hotchkiss  and 


HISTOAV: 


Capt.  Steven  A.  Hotchkiss  and 
Capt  Jam i for  Fried 

This  project  mas  indertakon  as  a  thesis  project  for 
partial  fulfil  leant  of  requirements  for  an  NS  dagrm 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


.  NathLib; 

PROC  Exp 
BEGIN 

Rent(Arg)  F 

ITEM  Arg  F 
END 

39; 

PROC  ALog 
BEGIN 

Rent  (flrg)  F 

ITEM  Arg  F 
END 

39; 

PROC  ALog 10  Rente Arg)  F 

BEGIN 

ITEM  Arg  F 
END 

39; 

PROC  Sqrt 

RENT (Arg)  F 

BEGIN 

ITEM  flrg  F 
END 

39; 

PROC  Sin 

RENT(Xx)  A 

BEGIN 

ITEM  Xx  A 
END 

1,30; 

PROC  Cos 

RENT(Xx)  A 

BEGIN 

ITEM  Xx  A 

1.30; 

3>  I 


PROC  Tan 

RENT<Xx> 

R 

12,18 

BEGIN 

ITEM  Xx  fl 
END 

1,30; 

PflOC  Cot 

RENT(Xx) 

R 

12,18 

BEGIN 

ITCH  Xx  fl 
END 

1,30; 

PROC  ftSin 

RENT(Xx) 

R 

1,30 

BEGIN 

ITEM  Xx  fl 
END 

1,30; 

PROC  RCo* 

RENT(Xx) 

R 

1,30 

BEGIN 

ITCH  Xx  fl 
END 

1,30; 

PROC  flTan 

RENT<Xx> 

R 

1,30 

BEGIN 

ITCH  Xx  A 
END 

1,30; 

PROC  Sinf 

RENT(Xx) 

F 

39 

BEGIN 

ITCH  Xx  F 
END 

39; 

PROC  Cosf 

RENT(Xx) 

F 

39 

BEGIN 

ITCH  Xx  F 
END 

39; 

PROC  Tanf 

RENT<Xx> 

F 

39 

BEGIN 

ITCH  Xx  F 
END 

39; 

'  PROC  Cotf 

REHT<Xx> 

F 

39 

BEGIN 

ITCH  Xx  F 
Bf) 

39; 

'  PROC  RSinf 

RENT(Xx) 

F 

39 

BEGIN 

ITCH  Xx  F 
END 

39; 

■  PROC  RCosf 

RENT<Xx> 

F 

39 

DATE:  29  August  1985 

VERSION:  1.0 

NAME:  loRefs 

MODULE  NUMBER:  1.0 

DESCRIPTION: 

This  Cospool  is  necessary  to  reference  routines  that  were 
necessary  for  testing  and  performance  evaluation  of  all  math 
functions  developed  for  the  1750. 

PASSED  UAAIABLES:  N/A 

RETURNS:  N/A 

MODULES  CALLED:  N/A 

AUTHOR:  Capt.  Steven  A.  Hotchkiss  and 

Capt.  Jennifer  Fried 

HISTORV:  This  project  eas  undertaken  as  a  thesis  project  far 

partial  fulfillment  of  requirements  for  an  MS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Wright  Patterson  AFB,Qh. 


START 


COMPOOL  loRefs; 


'  The  following  ITENs  are  required  to  print  a  carriage  return  and 

'  line  feed  on  a  terminal  connected  to  a  MIL-STD-1750  computer 
«• 


DEF  ITEM  Carriage  STRTIC  U  16  -  2573; 

DEF  ITEM  CRLF  STATIC  C  2; 

OUERLAV  Carriage:  CRLF; 


The  fol laming  referenced  subroutine  is  eritten  in  1750  Assembly  language 
and  is  used  to  print  character  strings  only.  Noncharacter  types  will 
have  to  be  converted  before  calling  this  routine.  The  following  DEFINE  is 
recommended  for  all  routines  calling  ObcSim: 

DEFINE  WRITE ’STRINGS)  ’  ,Printc<UOAOSIZE<!A>,LOC<!A»,  *; 

An  example  of  a  typical  call  follows: 

ITEM  Example  C  2; 


URITE'STRIMG(Example); 


REF  PROC  Printc  RENT (Length,  Message); 


ITEM  Langth  U  <BITSINU0RD-1 >; 
ITEM  Massaga  P; 

EMC 


Tha  fol loving  rmlmrmn cad  rout ina  is  nacasscry  for  routinas  aishing 
to  conuart  floating-point  uaiuas  to  a  charaetar  string 


’  PHOC  FltToChor  <flrg>  C  20; 
BEGIN 

ITEM  firg  F  39; 

END 


Tha  foi losing  rafarancad  rout ina  is  nacasscry  for  routinas  aishing 
to  conuart  fixad-point  valuas  to  a  charaetar  string.  Tha  war iabla 
IntOuarlay  oust  ba  ouarlayad  on  top  of  a  fixad-point  uarlabla  and 
BitsinFrac  is  an  intagar  uaiua  indicating  tha  nuabar  of  fractional 
bits  in  tha  fixad-point  valua. 


:  PflOC  FixToChcr  (IntOowlay,  BitsinFrac)  C  20; 
BEGIN 

ITEM  IntOuarlay  8  31; 

ITEM  BitsinFrac  U  8; 


i 


DATE:  29  August  1985 

UERSION:  1.0 

MATE:  FixToChar 

MODULE  NUMBER:  1.0 

DESCRIPTION: 

This  routine  is  used  to  convert  fixed-point  values  into 
character  representation.  This  routine  eas  necessary  for 
testing  and  performance  evaluation  of  math  routines  developed 
for  the  1750 

PASSED  VARIABLES:  IntOverlay  -  An  Integer  Variable  Over  I aged  on  top 

of  a  fixed-point  value 

BitsInFrac  -  the  number  of  fractional  bits  of 
the  fixed-point  argument 

RETURNS:  a  20  character  representation  of  the  argument 

NOOULES  CALLED:  FltToChar 

AUTHOR:  Capt.  Steven  A.  Hotchkiss  and 

Capt.  Jennifer  Fried 

HISTORV:  This  project  mas  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  MS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


PROC  FltToChar  <Arg)  C  20; 
BEGIN 

ITEM  Arg  F  39; 

END 


FixToChar  Procedure 


PROC  FixToChar  (IntOverlay,  BitsInFrac)  C  7 
BEGIN 

ITEM  IntOverlay  S  31; 

ITEM  BitsInFrac  U  8; 

TABLE  Overlays  <0)  U  3; 

BEGIN 

ITEM  Arg  F  39  POS<0,0>; 

ITEM  ArgExp  S  7P0S(8,1); 

END 

Arg<0)  -  <*  F  39  *)<  IntOverlay  ); 

ArgExp(O)  -  ArgExp(O)  -  (*  3  7  *>(  BitsInFrac  ); 

FixToChar  ■  FI tToChar<Arg<0>) 


OflTE:  29  August  1989 

UERSION:  1.0 

NAME:  FltToChor 

MOOULE  NUMBER:  1.0 

DESCRIPTION: 

This  routine  is  used  to  convert  floating-point  values  into 
character  representation.  This  routine  eas  necessary  for 
testing  and  performance  evaluation  of  eath  routines  developed 
for  the  1750 

PR8SED  UARIABLES:  flrg  -  the  value  to  be  converted 

RETURNS:  a  20  character  representation  of  the  argument 

NODULES  CALLED:  none 

AUTHOR:  Copt.  Steven  A.  Hotchkiss  and 


H I STORY : 


Capt.  Steven  A.  Hotchkiss  and 
Capt.  Jennifer  Fried 

This  project  eas  undertaken  as  a  thesis  project  for 
partial  fulfillment  of  requirements  for  an  US  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASO 
Language  Control  Branch,  Wright  Patterson  AFB,0h. 


PR0C  FltToChor  <Arg)  C  20; 

BEGIN 

DEFINE  Ves 

1BT"; 

DEFINE  No 

IB'O’ "; 

ITEM  Arg 

F 

39; 

ITEM  Fraction 

F 

39; 

ITEM  Temp 

F 

39; 

ITEM  Result 

C 

20; 

ITEM  lx 

U 

8; 

ITEM  ly 

u 

8; 

ITEM  ExpCnt 

u 

8; 

ITEM  NegExp 

B; 

ITEM  CharUal 

u 

8; 

ITEM  CharRep 

c 

i; 

0UERLAV  CharRep:  CharUal ; 

ITEM  ZeroRep 

STATIC  C 

1  - 

*  0  * ; 

ITEM  ZeroUal 

STATIC  U 

8; 

0UERLAV  ZeroRep:  ZeroUal; 

CONSTANT  ITEM 

Zero 

F 

39  * 

0.0 

CONSTANT  ITEM 

One 

F 

39  - 

1.0; 

CONSTANT  ITEM 

TenFloat  F 

39  ■ 

10.0 

CONSTANT  ITEM 

PtFIve  F 

39  ■ 

0.5 

CONSTANT  ITEM 

PtOne  F 

39  ■ 

0. 1 

•  *r*T*r*y * * . » v t ■  u ^ * W* T 7 ■  u ^ ■  y ■  1 pjj wy m p 1 pf mwjuj ■  ■■  1  ■  1  ■■fv; y  .^., 


'  VL  V  *JT 


W.VV. 


r  .t  .— 


i 


i 


K 


3 

I  U 

I 


1 


B 


V 


Rosult  •  •  0.0000000000000E+00‘, 

IF  flrg  <  2*ro; 

BEGIN 

Fraction  ■  -flrg; 

BVTE<B«wltJ0/ 1)  « 

END 

ELSE 

Fraction  ■  flrg; 

IF  Fraction  <  PtOna; 

NagExp  -  Vos; 

ELSE 

NagExp  ■  No; 


ExpCnt  ■  0; 

WHILE  (Fraction  >  Or*); 

BEGIN 

ExpCnt  •  ExpCnt  +  1; 

Fraction  *  Fraction  /  TanFloat; 

END 

IF  (NagExp  ■  Vos)  AND  (  Fraction  <>  Zaro); 

BEGIN 

BVTE(Rasul t, 17, 1 )  - 

WHILE  (Fraction  <  PtOna); 

BEGIN 

ExpCnt  ■  ExpCnt  +  I; 

Fraction  ■  Fraction  *  TanFloat; 

END 

END 

ly  «  0; 

WHILE  ((Fraction  <>  Z«ro>  AND  (ly  <  13)); 

BEGIN 

Top  »  Fraction  *  TanFloat; 

IF  ly  »  12; 

Taap  ■  Taup  +  PtFiwa; 

CharUal  -  (*  U  8  *)(  Taup  ); 

Fraction  ■  T«np  -  (*  F  39  *)(  CharUal  >; 

CharUal  »  CharUal  ♦  ZaroUal; 

BVTE(Rtsu) t, ly+3, 1 )  *  CharRap; 
ly  «  ly  ♦  1; 

END 

CharUal  >  (*  U  8  *>(E xpCnt  MOO  10)  *  ZaroUa I; 

BVTE(R«ult,  19, 1)  -  Charflap; 

CharUal  «  (*  U  8  *>(ExpCnt  /  10)  +  ZaroUal; 

BVTE(Ra*ul t, 18, 1 )  ■  Charfiap; 

FltToChar  ■  Result; 

RETURN; 

END 

TERN 


t 
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TITLE  HOL(PRINTC) 

MOOULE  PR1NTC 


DATE: 

VERSION: 

UAMC  • 

Iff  IL . 

NOOULE  NUMBER 
DESCRIPTION: 


4  Sapteaber  1985 

1.0 

Printc 

1.0 


This  aodula  is  called  to  print  a  character  string  onto 
a  consols  that  is  connected  to  a  Mi  l-Std-1750  coaputar 

PASSED  VARIABLES:  .  .  ..  _ 

LENGTH-3  -  this  variable  contains  a  count  of  the  number 
character*  to  print 

MESSAGE-3  -  this  is  a  location  pointer  for  the  string  to  be 
printed 


RETURNS: 

MOOULES  CALLED: 
AUTHOR: 

HISTORV: 


prints  eessoges  on  user  consols 

Copt .  Steoen  A.  Hotchkiss  and 
Copt.  Jennifer  Friad 

This  pro j act  acts  undartakan  as  a  thasis  project  for 
partial  fulfil  leant  of  raquiraaants  for  an  US  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 


*  $  4-SEP-85/16: 09: 29  $ 

PR INTOFF 

* 

*  START  OF  META  DEFINITIONS 


ORTAS 

LF<0> 


NORMA 


_J DRTAS 


SECTION 

ic<DL)N*<_) 


META 

EQU 

LOOP 

VOID 

GOTO 

LABEL 

DATA 

LOOPTEST 

MEMO 

META 

LOOP 

DATA 

LOOPTEST 

MEND 

LENGTH 

META 

LOOP 

CSECT 

LOOPTEST 

MEMO 


DO  MOT  LIST  HETAS 


HkFLHrED  PRESET  META 


2, 1,NUM<GF>-1 

GF<_,  1),  NORMA, -J3ATAS 

TEST 

GF<-> 


1, 1,GF<_,  1) 
GF<_> 


25, 9399 
0 

2,1,31 


CSECT  META 


lit 

LOOPTEST 

MEND 

♦  END  OF  META  DEFINITIONS 

* 

* 

*  BASE  RED  EQUATES 

ik 

B12 

EQU 

12 

B13 

EQU 

13 

B14 

EQU 

14 

B15 

EQU 

15 

* 

*  CONDITION  COOE  EQUATES 

* 

_NOP 

EQU 

0 

-LT 

EQU 

1 

_EQ 

EQU 

2 

JLE 

EQU 

3 

JOT 

EQU 

4 

_NE 

EQU 

5 

JOE 

EQU 

6 

JCY 

EQU 

8 

jclt 

EQU 

9 

JCEQ 

EQU 

10 

_CLE 

EQU 

11 

-COT 

EQU 

12 

JCNE 

EQU 

13 

JCGE 

EQU 

14 

jun 

EQU 

15 

* 

*  END  OF  EQUATES 

* 

REG 

SECTION 

PRINT 

DEFINE 

PRINTC 

PSSOATAI 

EQU 

3 

PSSCONSS 

EQU 

4 

PSSCOOEI 

EQU 

2 

*  NO  REF  OflTfl  DECLARATIONS 

*  MO  BYREF/TYPE/ABSOLUTE  DECLARATIONS 

*  LOCAL  AUTOMATIC  DATA  ***  SIZE  IN  UOROS  —  2  00 

*  LOCAL  flUTOMRUC  OflTfl  FOR  PflOC  PRIHTC 

*  STACK  FRflME  ***  SIZE  IN  UOflDS  —  2  DECIMAL  :  2 


BK-J003EF 

EQU 

HEX<0) 

LENGTH-3 

EQU 

HEX(O) 

MESSAGE-3 

EQU 

HEX<1> 

*  END  OF  LOCAL  AUTOMATIC  DECLARATIONS 

*  PSECT  $OATA  IS  EMPTY 


I 


I 

I 

1 

I 

I 

! 


I  *J 


* 

* 

*  R2  «  NUMBER  OF  CHARACTERS  IN  STRING 

*  R3  ■  LOCATION  OF  CHARACTER  STRING 

* 


* 


PSSCOOE 

ORIGIN 

HEX(O) 

PRINTC 

EQU 

$ 

41 

ORIGIN 

HEX (0002) 

AISP 

R2J 

.  ADJUST  CHARACTER  COUNT 

SAA 

R2J 

.  18T  TWO  COMMANDS  EQUIUALENT  TO 

* 

R2  -  R0UHD(R2/2> 

BLE 

LBL-D002 

.  BRANCH  OUT  IF  ILLEGAL  CHAR  COUNT 

OUTPUT 

EQU 

$ 

XIO 

R9,  ACS 

.  AERO  CONSOLE  STATUS 

TBR 

1,R9 

.  CHECK  STATUS  BIT  1 

BEZ 

OUTPUT 

.  IF  OFF,  LOOP  BACK  UNTIL  CONSOLE  READ 

L 

R9,0,R3 

.  GET  NEXT  TUO  CHARACTERS  OF  MESSAGE 

XIO 

R9,C0 

.  PRINT  BOTH  CHARACTERS 

AISP 

R3,1 

.  POINT  TO  NEXT  TUO  CHARACTERS 

£ 

SOJ 

R2,  OUTPUT 

.  DECREMENT  LOOP  COUNT,  GO  BACK  IF  HOF 

LB _ 0002 

EQU 

$ 

AISP 

R19,2 

on DM 

rUrn 

R2,R3 

URS 

R19 

ORIGIN 

HEX(OOOO) 

PSHM 

R2,R3 

SIS P 

Af9,2 

ORIGIN 

HEX <00 13 ) 

a 


OflTE:  10  October  1985 

VERSION:  1.0 

NRME:  Ref Mat 

MOOULE  NUMBER:  1 

DESCRIPTION1 

This  routine  is  used  to  convert  ITS  LINK  files  into 
a  format  that  can  be  loaded  into  the  SPERRV  1531 
computer  <1750R  arch i tec ture).  The  ITS  files  are 
'.SO'  files  and  oust  be  in  the  80  column  record  format 
described  in  the  EURO  ITS  Load  Nodule  ICO  <C0RL  *1005 
contract  *E33657-83-C-0244>.  Use  of  the  command  file 
LINK1750.C0N  to  link  all  compiled  modules  ei 1 1  insure 
that  these  records  are  of  the  right  format.  The  format 
of  the  SPERRV  loader  records  are  defined  in  Rppendix  B 
of  its  programmer  reference  manual .  The  bgtms  of  al I 
binary  data  fields  must  be  swapped  (i.e.  the  high 
order  bits  of  a  word  are  swapped  with  the  lorn  order  8) 
The  only  type  ITS  records  converted  are  binary  and 
end  record  types.  It  also  ignores  all  protection 
indicators,  and  can  not  handle  expanded  memory  jobs. 
Uhen  all  object  files  are  copied  into  a  single  object 
far  linking  by  the  ITS  LINKER,  the  main  procedure  must 
be  copied  into  Vie  file  first!!!!!!!  Otherwise,  this 
application  will  have  no  way  of  determining  the  point 
that  execution  is  to  begin.  The  'end'  record  created 
by  the  ITS  linker  contains  the  lowest  addrress  of  the 
load  module,  and  this  application  assumes  that  the 
routine  begins  at  that  point.  The  ITS  file  contains 
datafields  that  are  in  HEX  character  representation, 
and  the  SPERRV  1631  expects  binary  data  fields; 
therefore  the  ITS  data  must  also  be  converted  to 
binary 

PRSSED  VARIABLES:  N/fl 
RETURNS:  N/fl 

NODULES  CALLED:  GetHdr 

Readf 
Printf 
ClnUp 
IntFil 
UriteRcd 

AUTHOR:  Copt.  Steven  R.  Hotchkiss  and 

Copt.  Jennifer  Fried 

HISTORY:  This  project  was  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  NS  degree 
in  Information  Science  from  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  ASD 
Language  Control  Branch,  Uright  Patterson  RFB,0h. 


START 

IC0MP00L  <' loOata' ); 
IC0MP00L  (' loCal Is' ); 
IC0MP00L  < 'RfHtCpI 1 >; 
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PR0GRAT1  Raffia  t; 


BEGIN 

ChkSua  ■  48*0000'; 

FirstPoss  *  Trua; 

LdPt  -  -1; 

Eof  ■  Falsa; 

Buff  -  0; 

BufPtr<0>  -  1; 

BufPtr<1>  ■  1; 

■  Initial iza  10  Filas  * 

IntFil; 

“  Gat  Haador  Info  for  Loadar  FI  la  " 

GatHdr; 

WHILE  NOT  Eof; 

BEGIN 

"  Road  tha  first  80  coluan  racord  " 

Raadf < : I tsRcd,Eof ); 

“  Put  Loadar  Info  into  Contiguous  Naaory  Locations  " 
UdsInRcd  -  CntlCO)  -  RscilO; 

Adcfr'C(O)  ■  flddr<0>; 

Ud1C<0>  •  Wd1<0); 

Ud2C(0>  ■  Ud2<0>; 

Ud3C(0>  •  Ud3(0); 

Ud4C<0)  «  Ud4<0); 

Ud9C(0>  -  WdS<0);  . 

UdBC(O)  -  Ud0<0); 

Ud7C<0)  -  Ud7<0); 

■  Initial  iza  tha  Output  Buffars  11 
FOR  lx:  1  BV  1  WHILE  lx<33; 

CharToBinOx)  ■  0; 

FOR  lx:  0  BV  1  WHILE  lx<63; 

OutBuff(lx)  ■  0; 


“  Conwart  Char  To  Bin  and  Pack  it  " 

FOR  lx:  1  BV  1  WHILE  lx<«32; 

BEGIN 

IF  <RscliO<-CharToBln<lx»  RNO  <CharToBin<  IxX-flsci  I9>; 

CharToBinOx)  ■  CharToBln<lx)  -  RsciiO; 
n  SE 

IF  <Rsciifi<-CharToBin<lx>>  AND  <CharToBln< IxX-flscliF), 
CharToBinOx)  ■  CharToBinOx)  -  Rsciifl  ♦  10; 
HalfBytaOx)  •  NibblasOx); 

END 


IF  Typ<0)  ■  *  *; 

BEGIN  "This  is  a  binary  racord" 

IF  BufPtr<Buff)  +  UdsInRcd  <-01  AND  LdPT  -  Loddr<0>; 


BEGIN  ”  Old  Record  and  still  rooa  for  •ora  data  fields  " 

"  Flip  Flop  the  position  of  each  but*  of  a  1730A  word  " 

FOR  lx:  0  BV  1  UHILE  lx<UdslnRcd; 

BEGIN 

Buf8yt*0<BufPtr<Buff>+lx>  ■  FieldLC lx+1 >; 

BufBy  t*  1  (BufPtr  <Buf f H I x >  >  Fi*ldH<lx+1); 

END 

"  Point  to  share  info  froa  naxt  ITS  80  coluan  record  " 

"  is  to  be  placed  into  this  loader  record 
BufPtr(Buff)  •  BufPtKBuff)  ♦  UdsInRcd; 

*  Update  load  point  so  the  next  ITS  record  can  be  checked  to  “ 
"  see  If  it  belongs  in  this  loader  record 
LdPt  -  LdPt  ♦  UdsInRcd; 

IF  BufPtr<Buff>  -  61; 

BEGIN  *  Loader  Record  is  full  and  needs  to  be  eritten  “ 

UdsInBuffer  -  60; 

AcdTypI  ■  Rsci iB; 

UriteRcd; 

END 

ETC 

ELSE 

BEGIN  '  Old  record  and  not  enou^t  rooe  —  or  nee  record  " 

IF  LdPt  -  Laddr<0); 

BEGIN  "  Saaa  loader  record,  but  not  enough  rooe  for  all  " 

"  data  fields  in  ITS  record 

’  Soap  Bytes  of  eords  going  into  loader  record  * 

FOR  lx:  0  BV  1  UHILE  BufPtrtBuff Hlx  <  61; 

BEGIN 

BufByt*0<8ufPtr<Buff Hlx>  -  FieldL<  lx+1  >; 
BufByte1<BufPtr<Buff Hlx>  -  FieldH< lx+1 ); 

LdPt  »  LdPt  ♦  I; 

ETC 

9  srite  the  full  record  out  9 
UdsInBuffer  ■  60; 

RcdTypI  ■  Rsci iB; 

Uri teRcd; 

9  Sat  the  load  point  for  this  nee  loader  record  9 
LdRd(O)  ■  LdPt; 

9  Seap  bytes  of  the  other  ITS  data  fields  and  place  than  into 
9  record.  If  the  next  ITS  record  doesn't  have  the  load  point 
9  computed  here,  it  should  be  the  first  entries  for  another 
I oader  r  ecord 

FOR  ly:  lx  BV  1  UHILE  ly  <  UdsInRcd; 

BEGIN 

8ufByteO(BufPtr(Buff >+ly)  ■  Fi*ldL(ly+1); 
BufByt*1<BufPtr<Buff >+ly>  -  Fi*ldH<ly+1>; 

LdPt  ■  LdPt  ♦  1; 


END  "Soee  record  not  enogh  rooa 


ELSE 

BEGIN  "  this  is  the  start  of  a  nee  loader  record  " 

IF  NOT  FirstPass; 

BEGIN 

IF  BufPtr<Buff)  <>  1; 

BEGIN  "the  last  record  didn't  get  filled  up,  so  it 
"  hasn't  been  eritten  yet.  The  routine 
"  UriteRcd  sets  BufPtr  to  1  before  exit  " 

RcdTypI  ■  flsci IB; 

UdsinBuffer  >  BufPtr<Buff>  -1; 

UriteRcd; 

END 

END  "  end  not  first  pass  " 

FirstPass  *  Faise; 

"  Set  the  load  point  for  this  loader  record  " 

LdRdCO)  -  Laddr<0>; 

"  Seap  bytes  of  ITS  data  fields  going  Into  loader  record 
FOR  lx:  0  BY  1  WHILE  lx  <  UdsInRcd; 

BEGIN 

BufBy teO< I x+ 1 >  -  FieldL< lx+1 >; 

BufByteK  lx+1 )  ■  FieldHC lx+1 ); 

END 

BufPtrCBuff)  -  UdsInRcd  +  1; 

LdPt  -  LoddrCO)  +  UdsInRcd; 

END  "  end  nee  record  " 

END  "end  of  old  record  not  enough  rooe  —  or  nee  record  " 
END  "end  of  this  is  a  binary  record" 


BEGIN  "this  is  an  execution  address  record  " 
IF  Typ<0>  -  'E'; 

BEGIN 

RcdTypI  ■  8261;  ‘blank  E" 

OutBuff(O)  -  LoddrCO); 

OutBuffC 1 )  >30;  "  ascii  record  seperator  " 

UriteRcd; 

END 


END  "end  execution  address  record  " 
D  "end  ehi le  loop  " 

Urite  end  of  file  loader  record  " 


DOTE:  10  October  1083 

UERSION:  1.0 

NAME:  Uritefied 

MOOULE  NUMBER.  7 
DESCRIPTION: 

This  routine  is  called  by  RefMat  to  do  10  stuff  that 
needs  to  be  done  throughout  the  ia i n  procedure .  Three 
types  of  SPERRY  1631  loader  records  are  eritten: 

Binary  , Execution,  and  End  of  file.  If  the  record  type 
is  a  binary  record,  this  routine  coeputes  a  checksum 
for  it  and  teks  it  on  to  the  end  of  the  record.  Then 
the  record  type  is  eritten  out  folloeed  by  the  binary 
record,  if  the  record  type  is  an  execution  record  or 
an  end  of  file  record,  the  record  type  is  eritten  out 
folloeed  by  the  record.  The  variable  'Buff'  is  a 
global  variable  that  points  to  the  record  to  be 
eri tten. 

PflSSED  VARIABLES:  None 
RETURNS:  Nothing 

MODULES  CflLLED:  Printf  —  a  FORTRAN  10  routine 

AUTHOR:  Copt.  Steven  R.  Hotchkiss  and 

Capt.  Jennifer  Fried 

HISTORY:  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfil  leant  of  requirements  for  an  MS  degree 
in  Information  Science  froe  the  Rir  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  R8D 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 

eeeeeeeeeeiMeeenreiii  i  MitiMteeeeeeeemment  1 1  \i*+*v*+**+++m*+—**+**+**********> 


!CQMP0QL< 'RfMtCpI ' ); 

1C0MP00L  <* loOata*  >; 

REF  PROC  Printf (RcdTyp, Buffer); 
ILIMKRGE  FORTRAN; 

BEGIN 

ITEM  RcdTyp  S  15; 

ITEM  Buffer  C  126; 

END 

OEF  PROC  UriteRcd; 

BEGIN 

LoopCnt  -  UdsInBuffer; 

If  RcdTyp I  >  flsci iB; 

BEGIN 

ChkSue  -  46*0000'; 

0utBuff<0)  -  LdAd<0); 

ChkSue  -  ChkSue  XOR  OutBuffB(O); 


OutBuff(l)  -  UdsInBuffer; 

ChkSue  -  ChkSue  XOR  OutBuf fB< 1 ); 


FOR  lx:  1  BV  1  WHILE  lx  <-  LoopCnt; 
BEGIN 

0UTBUFF< I x+ 1  >  >  BufWddx); 

ChkSua  ■  ChkSua  XOR  0utBuffBdx+1> 
END 

OutBuffBCIx+l)  -  ChkSua; 

Pp i ntf (RcdTyp I , OutF I d >; 


ELSE 

Printf<RcdTgpl  ,OutFld>; 


BufPlr<Buff )  -  1; 
Buff  «  RBS< 1-Buff >; 


RETURN; 

END 


OEF  ITEM  RcdTyp  C  2; 

OUERLRV  RcdTyp I :  RcdTyp; 

OUERLRV  Infi I,  OutFil,  Filrwa,  Hu* fcr,  lt*Rcd,  (XitFId,  Eof,  RcdTyp 


*  DflTE: 

*  UERSIOM: 

*  NAME: 

*  MODULE  NUMBER: 

*  DESCRIPTION: 


PRSSED  URRIRBLES 
RETURNS: 

MOOULES  CALLED: 
AUTHOR: 

HISTORY: 


10  October  1985  * 

1.0  * 

RfMtCpI  ♦ 

10  * 

* 

This  coapool  contains  all  the  variables  and  tables  * 

that  are  used  to  unpack  ITS  I  inker  records,  packs  thee  * 
and  converts  the  HEX  characters  to  binary  data  fields,  * 
and  then  places  thee  into  a  SPERRY  1631  loader  record  * 

format  * 

N/fl  * 

N/A  * 

N/A  * 

Copt.  Steven  A.  Hotchkiss  and  * 

Copt.  Jennifer  Fried  * 

This  project  eas  undertaken  as  a  thesis  project  for  * 

partial  fulfil  leant  of  requirements  for  an  MS  degree  * 

in  Information  Science  from  the  Air  Force  Institute  * 

of  Technology.  Sponsoring  organization  is  the  RSD  * 

Language  Control  Branch,  Wright  Patterson  AFB,0h.  * 


COMPOCL  RfMtCpI; 


ITEM 

ChkSue 

B 

16, 

ITEM 

FirstPass 

B 

1, 

ITEM 

LdPt 

S 

15 

ITEM 

Buff 

U 

8 

ITEM 

lx 

U 

8 

ITEM 

<y 

U 

8, 

ITEM 

UdsInRcd 

S 

13, 

ITEM 

LoopCnt 

S 

15 

ITEM 

UdsInBuffer  S 

15, 

OEF  ITEM  Zero  STATIC  C  1  »  *0’ 

OEF  ITEM  AsciiO  STATIC  S  7; 

OUERLAV  Zero:  Asci  iO; 

OEF  ITEM  Mine  STATIC  C  1  ■  ’9’ 

OEF  ITEM  Asci  (9  STATIC  S  7; 

OUERLAY  Nine:  Ascii9; 

OEF  ITEM  AA  STATIC  C  1  -  ‘A’ 

OEF  ITEM  Asci iA  STATIC  S  7; 

OUERLAV  AR:  Asci iA; 

OEF  ITEM  FF  STATIC  C  1  ■  T’ 

OEF  ITEM  Asci iF  STATIC  S  7; 

OUERLAY  FF:  AsciiF; 

OEF  ITEM  B8  STATIC  C  2  -  *  E 

OEF  ITEM  Asci iB  STATIC  S  15; 

OUERLAY  B8:  AsciiB; 


3EF  TRBLE  LoodPoint  CO); 
BEGIN 

I TEN  Ldftd  S  IS; 

END 


DEF  TABLE  BufStuf  (0:1); 
BEGIN 

I TEH  BufPtr  S  7; 

END 


DEF  TABLE  PockedRcd 
BEGIN 

I TEN  AddrC  C  4 

I TEH  UdIC  C  4 

I  TEH  Ud2C  C  4 

I  TEH  Ud3C  C  4 

ITEM  Ud4C  C  4 

ITEM  Ud5C  C  4 

ITEM  Ud6C  C  4 

ITEM  UdTC  C  4 

END 


<0>  U  8; 

POSCO,  0), 
POSCO, 1) 
POSCO, 2) 
POSCO, 3) 
P0SC0,4) 
POSCO,  5) 
POSCO, 6) 
POSCO,?) 


DEF  TABLE  CharConvert  < 1:32)  T  8  U; 
BEGIN 

ITEM  CharToBIrt  S  7  POSCO, 0); 

ITEM  Nibbles  S  3  POSC4,0>; 

END 

OUERLAV  PackedRcd :  CharConvert; 


DEF  TABLE  HexBuf  <1:32>  T  4  U; 

BEGIN 

ITEM  HalfByte  8  3  POSCO, 0); 

END 

DEF  TABLE  Paklts  <0:7>  T  16  H; 

BEGIN 

ITEM  Laddr  S  15  POSCO, 0); 

END 

OEF  TABLE  BfnFields  <0:7>  T  (6  H; 
BEGIN 

ITEM  Field  S  15  P0S(0,0>; 

ITEM  FieldH  S  7  POSCO, 0); 

ITEM  FieldL  S  7  P0SC8,0); 

END 

OUEALAY  HexBuf, lx:  Paklts:  BinFields; 

OEF  TABLE  DatFields  CO)  U  1; 

BEGIN 

ITEM  BufByteO  S  7  POSCO, 0); 

ITEM  BufBytel  S  7  P0SC8,0); 

ITEM  BufUd  S  15  POSCO, 0); 


ooooooooooooooooooooooooo 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC^^ 

10  October  1989 
1.0 
IntFil 

6 

This  routine  is  cal  lad  by  tha  JOUIflL  routina  cal  lad 
RefNat.  its  purposa  is  to  proapt  tha  usar  for  tha 
noaa  of  a  f i la  that  oas  craatod  by  an  ITS  link, 
proapt  tha  usar  for  tha  naaa  of  a  f i la  that  tha 
reforaottad  ITS  fila  is  to  ba  erittan  to,  and  than 
opans  both  filas.  Tha  input  f i L ^  aust  ba  a  *.S0"  fila 
and  tha  output  fila  is  a  “.OAT"  fila. 

PASSED  VARIABLES:  Nona 
RETURNS:  Nothing 

GLOBAL  VARIABLES:  All  uarioblas  usad  ora  global,  and  haua  baan 
dafinad  in  tha  coaaon  (CQMPOOL)  cal  lad  loOata 
NODULES  CALLED:  Nona 

AUTHOR:  Copt.  Stavan  A.  Hotchkiss  and 

Capt.  Jormi far  Friad 

HISTORV:  This  pro j act  aas  undartakan  as  a  thasis  pro j act  for 

partial  fuifiiiaant  of  raquiraaants  for  an  NS  dagraa 
in  Inforaation  Sclanca  from  tha  Air  Forca  Institute 
of  Technology.  Sponsoring  organization  is  tha  ASO 
Language  Control  Branch,  Uright  Patterson  AFB,0h. 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC^^ 


DATE: 

VERSION: 

UAUC  . 

URL. 

NODULE  NUNBER: 
DESCRIPTION: 


Subroutine  IntFil 

INPLICIT  INTEGER  <A-Z) 

CHARACTERS  Filnaa 
CHARACTER* 10  Infil,  Outfit 

URITE<*,*)'  Enter  Fila  Nona  <Hox  6  Characters)  1 
REA0(*, 10)Fi Inan 
10  FORMAT  <R6 ) 

I  «  INDEX(Fi Inaa, ' . ' )  -  1 
IF  CI.LE.O)  THEN 

I  »  INDEX(Fi Inaa, '  * >  -  1 
IF  (I.LE.O)  THEN 
I  •  6 
END  IF 
END  IF 

Infil  ■  Filnoad:  !>//'.  SO' 

Outfil  -  Fi lnaa(1: I )//* .OAT' 

URITE<*,*)’ Input  Fila  »  ’, Inf i I , ’Output  fila  »  ’, Outfil 

OPENCUNIT  -  2,  NAME  -  Infil,  TVPE  -  ’OLD’,  FOAM  -  ’FORMATTED’ > 
OPEMUNIT  »  3,  NAME  ■  OutFII,  TVPE  ■  ’NEU’, 

1F0AM  ■  ’UNFORMATTED’) 

END 


141 


ooooooooooooooooooooooooo 


oooooooooooooooooooo 


DATE:  10  October  1989 

UERSION:  1.0 

NAME:  QetHdr 

NOOULE  NUMBER:  2 

DESCRIPTION: 

This  routine  per fores  10  for  a  JOUIRL  routine  called 
RefMat.  It  requests  a  user  to  input  a  one  line 
header  that  eill  be  placed  in  a  loader  file. 

PRSSED  URRIRBLES:  None 

RETURNS:  Nothing 

OLOBRL  URRIRBLES:  Rll  variables  used  are  global,  and  are  defined  in 
the  coeeon  (C0NP00L)  called  loOata 
NOOULES  CRLLED:  None 

AUTHOR :  Copt.  Steven  R.  Hotchkiss  and 

_  Capt.  Jennifer  Fried 

HISTORY:  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfil  leant  of  roquiroeonts  for  an  NS  degree 
in  Information  Science  from  the  Rir  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  RSD 
Language  Control  Branch,  Uright  Patterson  flFB.Oh. 


oooooooooooooooooooo 


nooooooooooooooooooooo 


oooooooooooooooooooooo 


ccccccccccccccccccccccccccccccccccccccccc^^ 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


OfiTE: 

VERSION: 

NAME: 

MOOULE  NUMBER: 
DESCRIPTION: 


10  October 
1.0 
Read  f 
3 


1989 


This  rout Ins  is  call'd  by  the  JOUIflL  routin' 
RefMat.  Its  is  used  to  rood  80  coluan  record 
by  the  ITS  linker.  Th«  no"  of  the  fi I  being 
stor'd  in  the  global  variable  Infil  ehich 


PRSSED  VARIABLES: 
RETURNS: 

GLOBRL  URRIRBLES: 


cal  I'd 
:  creat'd 
read  is 

set  in 

the  routin'  call'd  IntFil.  This  f i i'  is  a  *.80“  f i l« 


MOOULES 

AUTHOR: 

HISTORV: 


CALLED: 


Nothing 

All  variables  used  are  global,  and  are  defined  in 
the  co— on  (COMPOQL)  call'd  loOata 


for 


Copt.  Steven  A.  Hotchkiss  and 
Capt.  Jennifer  Fried 

This  project  «as  undertaken  as  a  thesis  project 
partial  fulfil  leant  of  require— nts  for  an  MS 
in  Inforeation  Science  froe  the  Air  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  A8D 
Language  Control  Branch,  klright  Patterson  AFB,0h. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ccccccccccccccccccccccccccccccccccccccccccccccccc^^ 


II 


10 


20 

30 


Subroutine  Readf < I tsRcd,  Eof) 

IMPLICIT  INTEGER  <A-Z> 

CHARRCTER*80  ItsAcd 
LOG  I  CALM  Eof 

Eof  ■  .FALSE. 

READ<2, 10,  END  »  20)  ItsAcd 
FORMAT <A80) 

URITE<*,*)I tsRcd 
GOTO  30 
Eof  -  .TRUE. 

CONTINUE 

END 


ooooooooooooooooooo 


DfiTE:  10  October  1985 

VERSION:  1.0 

NRME:  ClnUp 

nOOULE  NUMBER:  5 

DESCRIPTION: 

This  routine  is  called  bg  the  routine  RefMat  to  close 
the  files  it  used  for  10 
PRSSED  VARIABLES:  None 

RETURNS;  Nothing 

NOOULES  CALLED:  None 

GLOBAL  UfiRIRBLES:  Allvariables  used  are  global,  and  are  defined  In 
the  eoeeon  <C0NP00L)  called  foOata 
AUTHOR:  Copt.  Steven  R.  Hotchkiss  and 

_  Capt.  Jennifer  Fried 

HISTORV:  This  project  eas  undertaken  as  a  thesis  project  for 

partial  fulfillment  of  requirements  for  an  NS  degree 
in  Information  Science  from  the  Rir  Force  Institute 
of  Technology.  Sponsoring  organization  is  the  RSO 
Language  Control  Branch,  Wright  Patterson  RFB,0h. 


ooooooooooooooooooo 


■odul* 


!  A8H1790  —  Until  a  1790  sourc* 

I 

I  M8H1790  f i  I* 

I 

I  fi la  ■  input  source  naa*  of  *odu I*  fila.81 

I 

!  Croat*  fliinti I t  input  fli*  Ui  that  daslgnatas  mi-Std-l7S0A  as  tha  targat 
I  rathsr  than  th*  altamat*  179QA  targat 
I 

CREATE  ‘PI'.UI 
ASSEMBLE  TRAGET*N1730A 

tv  ‘pr.ui 

i 

ASSIGN  ‘PI'.UI  UI  I  UPDATE  INPUT  C0HHIC6  FILE  (INPUT) 

ASSIGN  *P1  * .81  81  I  1750A  AB8BBLY  SOURCE  FILE  (INPUT) 

ASSIGN  'PI*. OBJ  00  I  OBJECT  OUTPUT 

AB8IGN  'PI* .80  80  I  SYMBOLIC  OUTPUT 

ASSIGN  ‘P1‘.LO  LO  I  LIST  INS  OUTPUT 

ASSIGN  LIBUOUIAL.1790A  01  I  LIBRARY  INPUT 

I 

SET  VERIFY 

H1790A 

I 

0EA88IGN  81 
0EA88IGN  UI 
0EA88IGN  00 
0EA88IGN  80 
0EA88IGN  LO 
0EA88IGN  01 
I 

OELETE  ‘Pr.UI;* 

SET  NOUERIFV 


I  LINK 1730  —  Link  one  or  sort  17303  target  object  aodulas. 

I 

f  CL INK 1790  fl la 
I 

I  flla  ■  objact  file  (containing  one  or  aora  object  nodules) 

I  create  abject  file  by  first  deleting  all  .abj  files  for 

!  CONPOOLs  that  don't  oontain  aiy  OEFs.  Then  use  the 

i  fol losing  coeeands  to  create  the  object  file 

I 

I  COPY  *.0BJ  file.O 
I  RENRHE  file.O  flle.OBJ 
I 

I  OBJ  files  created  by  the  coepiler  and  the  assoabler  can  be  copied  to  the 
I  sane  OBJ  file,  but  the  UAX  sill  glue  an  incoapatible  files  earning.  Ignore 
!  the  earning,  the  copy  is  eade  anyeay 
I 

I  Create  Li r*er  irput  file  Ul" 

I 

SET  UERIFY 
I 

CREATE  'PI'.UI 
LINK  ORTH, LI ST, DEBUG, INPUTS 
ALLOCATE  LOCUTION- 1000  NODULES  . 

INKEND 

I 


ASSIGN  ‘PI'.UI 

Ul 

1  LINKER  CONTROL  (INPUT) 

ASSIGN  'PI' .OBJ 

00 

1  OBJECT  HOOULE(S)  (INPUT) 

ASSIGN  'PI'. 80 

SO 

1  LORO  NODULE  (OUTPUT) 

ASSIGN  'PI' .LO 

LO 

1  LINKER  LIST  FILE  (OUTPUT) 

A88IGN  LIB-J0UIAL-1730A 

01 

1  LIBRARY  OBJECT  FILE  (INPUT) 

I 

ITSLINK  I  RUN  1730R  Linker. .  .reads  logic  deulce  Ul 

I  Output  on  80  and  LO 

0ER8SIGN  Ul 
0ER88IGN  00 
0EAB8IGN  80 
0ER88IGN  LO 
0ER88IGH  01 
I 

DELETE  'PI'.UI;* 

SET  NOUERIFV 


I  LOGIN. COM  This  coooand  procedure  is  Inwokod  mi th  soch  login. 

I  and  Bay  bo  changed  to  tai  lor  yoir  any  Iron— nt. 

I 

I  Sot  standard  aliasos.  No to  that  sovaral  UNIX-1 iko  allasos  arm  sot  up. 
I 

SET  NOUERIFV 

SET  PROTECT  1 0tMSV8TEH:R,0UHER:  A®),  GROUP  :RU,  WORLD:  RUE) /DEFAULT 
8VHB0LS: 

BQ  SHOW  QUEUE /BATCH 

CD  SET  DEFAULT 

D8  OIAECTORY  /8I2E 

E  EDIT 

hfiir  qct  ncrai  t  rwytflnn  •  ronra  unrmjfioi 

LO  CL000UT.C0H 

LS  OIAECTOAV 

PQ  :—  SHOW  QUEUE  SYSlPAINT 

PS  :—  SHOU  PROCESS  I  Llko  UNIX  ps  coooand 

PIO  SHOU  DEFAULT  I  Llko  UNIX  pod  coooand 

R  RUN 

SO  SHOU  OEUICES 

SO  SHOU  SVNBOLS  /GLOBAL  /ALL 

ST  :—  SHOU  TEAHINAL 

UHO  SHOU  USERS  I  Llko  UNIX  oho  coooand 

SHQ  SHOU  QUEUE  SLAH$QUEUE/ALL 

S80  :—  SET  TEAHINAL/UIDTH-60 

S132  :—  SET  TEAHINAL/UIDTH-132 

JOG 1790  :—  SJ0U173O 
LINK  1730:—  SLINK  1730 
8 IN 1730  SS IN  1730 

ASH 1790  SASH  1730 

UNLOCK  SUNPAOTECT 

I 

I  End  uoor  do f I nod  koyins. 

I 

I  DEFINE  JOUIAL  LIBAAAV  FDA  AUTOMATIC  SEARCH  I  NO  FDR  UAX  TARGET 
I 

I  ASSIGN  JOULI0U: JQULIBU.OLB  LNMLIBAARV 
I 

I  Tho  fol  losing  dofinos  tho  1730A  support  tools  pioudo  coooands: 

I 

LINK90A  LINKITS 
AAIDX  STOOLS: RAID 

I 

I  END  LOGIN. CON 
I 

IFINISH: 


proc«kr«  TCHEBVSHEF_ECQMOniZRTIOH  Is 


This  procadura  is  tha  Min  dr  Ivor  for  ths  Tchabyshaf  aconoaizotion 
of  o  polytoalol. 


;ifli  3t  :  □ 


procoduro  DI8PUWJJECT0R  (PRINT  JJECTOR:  in  UECTOR)  is 
— Ths  sols  purposs  of  this  routine  is  to  display  an  intsgor  use tor 

package  INT-IO  is  nos  INTEGER- 10  ( intsgsr); 
use  INT.I0; 

begin  — Display  Doctor. 

for  I  In  0. .OEGflEELOFJ’OLVNOfllflL  loop 
put<l>; 
put  <"  -); 

put  (PRINT JUECTOR  (I  >>; 
nosi_i  ins; 
and  loop; 

end  0I8PUW JJECTOR: 


procedure  DI8PUfi'_FL0fiT_UECT0R  (PRIHTJJECTOR:  in  FLORTJJECTOR)  is 
— The  sols  purposa  of  this  rout ins  is  to  display  a  floating  point  woe tor 

package  INT-IO  is  noo  INTEGER- 10  (integer); 
use  INT-IO; 

package  FLT-IO  is  nos  FLOflT-IO  (float); 
use  FLT-IO; 

begin  — Oieplay  Float  Uector. 
for  I  in  0. . DEGREE jQF_pq.VtOimL  loop 
put  (I); 
put  ("  “); 

put  (PRINT JUECTOR  (I)); 
neo-line; 
end  loop; 


end  OI8PLRV-FLORTJUECTOR; 


procadr*  0I8PUWJWTRIX  (PfllNT-NRTRIX:  In  fIBTRIX)  is 
— Tht  sols  purposo  of  this  rout I no  is  to  display  an  intagor  oatrlx 

packogo  INT-IO  is  nos  INTEGER- 10  (Intagor); 
uso  I  NT-10; 

packogo  FLT-10  is  nao  FLQRT.IO  (float); 
uso  FLT-10; 

bogin  — Display  Natrix. 

far  I  in  0.  .OEGfiEEJOF-PGLYNOhlflL  loop 
put  (I); 
put  (“  •>; 

far  J  In  0. . OEOREE JOF-POLVHOH I flL  loop 
put  (FRItfTJWTRIX  (l,J»; 
put  C  "); 
and  loop; 
noo-lino; 
and  loop; 

and  DISPLRV-ttflTRIX: 


procedure  D»8PLflV_FL0flTJtflTBIX  (PRINT-MATRIX:  In  FLOflT-MBTAlX)  is 
—The  sole  purpose  of  this  routine  is  to  display  a  floating  point  oatrix 

package  IHT-10  is  nee  INTEQER-10  (integer); 
use  IMT-10; 

package  FLT-IO  is  nee  FLOAT- 10  (float); 
use  FLT-10; 

begin  —Display  Float  Matrix, 
for  I  in  0.  .DEOREEJF-POLVMOMiflL  loop 
put  (I); 
put  ("  ">; 

for  J  in  0.  .DEDflEEJF-FOLVNOHIflL  loop 
put  (PRINT-MATRIX  (l,J»; 
put  ("  ">; 
end  loop; 
nee-l  ine; 
end  loop; 

end  OISPLflV-FliMTJfrrRIX: 


igin  — 'Tchdbgahaf  Econoaization 


INPUT-COEFFICIENTS; 

put  <“ Input  Coafficionts'); 

n— _l  ina; 

OISPLflVJTJJflTJKECTOR  (COEFFICIENTS); 

CONPUTE-TCHEBV8HEF-POLVNONIRL; 
n— -I  lr»; 

put  ("Tchabgihaf  Polgnoalal"); 
naaulina; 

DISftflVJmTOIX  (TCHDVMFJPOLVNOHIALS); 

CONPUIEJNMER8JOFJTCHEBV8HEF; 

moulino; 

put  (“Posars  of  Toh«bgihofa  >; 
noa-Jino; 

OISFUW-FUMT-mTRIX  (POUER8-OF-TCHEBV8HEF >; 

— tianarota  tho  oork  aatrix  und  in  tha  final  calculations  of  tho  oconooizod 
—  polgnoalal .  Again  tho  aatrix  is  looar  triangular, 
for  I  In  0.  .OEGAEEJ3F_POLVNOniRL  loop 
for  J  in  0. . I  loop 

UOfKJVrmiX  <|,J>  :■  float (HULTIPLIER  <l»  •  P0UER8_0F_TCHEBVSHEF  <i,j> 
*  COEFFICIENTS  <l>; 
and  loop; 
and  loop; 

■ftecuaulata  tho  sua  of  tho  oork  aotrix  co loans 
for  I  In  0.  .OEBREEJOF-FOLYNOHIflL  loop 
for  J  in  0.  .OEQfEEJOF-POLYNOniAL  loop 
8UH(J)  StffKJ)  ♦  (UOAKjnmiX<l,U»; 
and  loop; 
and  loop; 

— Porfora  tha  final  additions  and  aul tipi  loot ions  to  fora  tho  rasult. 
for  I  in  0. .  (OEGAEEJDF-POLWOniflL  -  1)  loop 
for  J  in  0. . I  loop 

ECONOtllZED-POLVNOfllflL  (J)  ECONOHI2BLPOLYNOHIflL  (J>  ♦ 
flooterOEBVSNEF-FOLVNOHIALS  <I,J»  *  SUN  (I); 
and  loop; 
and  loop; 

nao.1  ina; 

put  ("Econoaizod  Polgnoaiol*); 
naoul  Ina; 

OI8PLAV-FLOATJJECTOR  (ECONOMIZED-POLVHONIflL); 


I  if  i  if  ii 


STRINOuTQ_INT 
IWUT-C0EFFICIENT8 
COJWTE_TCHEBV8HEF_POLVMOh  I  AL 
COnFUTE-POUER8dOF-.TCHEBVSHEF 
DI8PLAVJJECT0R 

Diflpuwjnmix 

DESCRIPTIONS:  Provided  aith  Mch  routine. 

PR88ED  URRIABLE8:  The  input  to  this  systaa  Is  ths  daseription  of 
tha  polyoaial  to  ba  aconoai zad . 

RETURNS:  Tha  result  of  procassing  is  tha  coafficiants  of  tha 
aconoai  zad  polynoalol. 

ORLLINO  MODULES:  TCHEBV8HEF-EC0N0NIZATI0N 
AUTHOR:  Copt  Jam)  far  Friad  and 
Copt  Steuen  Hotchkiss 
HISTORV:  Original  version,  Dae  1,  1088 


This  packoga  receives  tha  coafficiants  of  a  polyoaial  that  is  to  bo 
aconoai  zad,  coaputas  its  Tchabyahag  polynoalol,  and  tho  poaors  of 
Tchobyshaf  aatrlx. 


— Uncons  trai  nod  typo  dac  I  orations 

MATRIX  is  array  <intagar  ronga  <>,  intagar  ranga  <>>  of  intagar; 
latrlx  of  Intagar  ualuas,  usod  to  contain  tho  Tchobyshaf  polynoaials 
FLOAT-MATRIX  is  crray  < intagar  ronga  <>,  intagar  ranga  <>)  of  float; 
Matrix  of  floating  point  ualuas,  used  to  contain  tho  possrs  of 
Tchobyshaf 

VECTOR  is  array  (intagar  ronga  o)  of  intagar; 

Jactor  of  intagar  ualuas,  usod  to  contain  tho  nut  tipi  lor  of  tho  aatrlx 
FLOAT JJECTOR  is  cavay  (intagar  ranga  «>  of  float; 

—Vector  of  floating  point  ualuas,  usod  to  contain  tho  coafficiants  of 

—  tha  polynoalol 

— Variable  declarations 
MAX-DIGIT:  in toga'  :-  19; 

—Tho  ooxlaua  nuabar  of  digits  para  it  tad  in  a  nuabor  is  nine. 

—  This  ualue  roprosants  tha  aaxiaua  input  string  length  for  too  nuabora 

—  end  a  slash,  V". 

MAX-DEGREE:  intoga'  :-  9; 

—Tho  aaxiaua  ualuo  of  tho  largest  exponent  of  tho  polynoalol 
DEDREE-flF-PCLVNOniflL:  intagar  :•  0; 

— Tha  actual  ualuo  of  the  largest  exponent  as  input  by  tha  usor 
COEFFICIENTS:  FLOAT  JJECTOR  (0.  .MAX-DEGREE)  :■  (0.  .MAX-DEGREE  »  0.0); 

— Contains  a  coefficient  for  each  day  so  of  the  pofyoaioi  that  sas 

—  specified  by  tho  user 

MULTIPLIER:  VECTOR  (0.  .MAX-DEGREE)  :■  (0.  .MAX-DEGREE  »  0); 

—This  uactor  contains  the  reciprocal  of  tha  ualuas  contained  on 

—  tho  diagonal  of  the  Tchobyshaf  polytoalal  natrix. 

—  Usod  in  generating  tho  aconoai zad  polynoalol. 


i  ration 


■The  Matrix  obtained  ahon  using  tha  Tchabgshaf  foraula. 
R8U0F.TCMEBV8NEF:  FLORTJffTRIX  <0.  .HRXJEOflEE,  0.  .URDU 
(0.  .IVDLDEGREE  »  <0.  .mX-OEGREE  ■>  0.0»; 

■Tha  aatrix  for— d  ahon  applying  tha  sacond  stop  of  tha 
■  algor i tha 


function  STRIN8-TQ_INT  <8:  string)  return  intagar; 

— This  function  is  used  to  convart  tha  input  coofficiant  string  into  on 
—  intagar  valua  that  aquatas  to  tha  nuaarotor  and  tho  donooinator. 


Tha—  procaduras  perform  tho  functions  specif  lad  by  this  package 
procedure  IMPUT_OXFFICIENT8; 

—Get  tho  input  coefficients  for  tho  polynomial 
procedure  GOnPUTEJTCHEBVSHEF-PQLVNQniflL; 

— Ganarato  tho  Tchabyshef  polynomial  matrix 


function  8TRINGLT0LINT  <8:  string)  return  integer  is 
— String  to  integar  equivalent  conversion. 

CHRH  :  character;  —‘individual  nuabar  in  each 

—  placeholder  of  the  Input  string. 
DIGIT  :  Integar;  —Individual  nuabar  in  each  placeholder 

— »  of  thi  outxxjt  intiCMr . 

MULTIPLIER  :  Integer  :■  1;  -Tee  value  of  the  Integv* 

—  pointer. 

FINRLJE9ULT  :  Integar  :■  0;  —Output  integer  being  generated. 

P08ITI0N  :  integar  :■  S' last;  —Pointer  Into  input  string 

—  < coves  rltftt  to  left). 

begin  — String  to  integer  conversion. 

—Starting  froe  the  end  of  the  Input  string,  process  each 
—  successive  character  unti  I  al  I  characters  have  been  converted, 
chile  POSITION  >-  S' first  loop 

—Get  one  character  digit  froe  the  input  string, 
am  :■  S(PQSITION); 

— If  this  is  a  valid  character  digit  representation,  convert  the 
—  character  Into  Its  nuMric  representation,  and  eultlply  It  by 
—  Its  tons  value. 

If  am  in  'O' . . '0*  than 

OIGIT  :■  character 'pos<CHRR)  -  charactar'pasC '0' ); 

DIGIT  :»  DIGIT  •  MULTIPLIER; 

—If  the  final  value  el  1 1  be  the  cost  negative  nunbor, 

—  desipnte  It  as  the  aost  negative  nuabar  and  stop 

—  processing.  The  reason  this  is  done  Is  to  adjust  for  the 

—  problae  that  the  absolute  value  of  the  aost  negative  nuabar 

—  is  1  digit  larger  than  the  aost  positive  nunbor  and  sill 

—  result  in  an  out-of-bound  condition. 

if  integer* last  ■  (FINRL-RESULT  -  1)  ♦  DIGIT  thw 
FINRLJESULT  :■  In  tegs' '  first; 

POSITION  S'  first; 
else 

— Others  I  se,  this  Is  not  the  aost  negative  nunbor.  Thus, 

—  add  the  current  digit  to  the  rest  of  those  fowd,  and 
—  Incresant  the  tens  value  to  the  next  larger  nuabar. 
FINRL-RESULT  :•  FINRLJESULT  ♦  DIGIT; 

MULTIPLIER  MULTIPLIER  •  10; 
end  If; 

—If  the  original  Input  aas  negative,  than  negate  the  results, 
elsif  CHRR  -  '-'  then 
FINRL-RESULT  :■  -FINRL-RESULT; 
and  if; 

—Adjust  the  pointer  into  the  input  string  to  point  to  the  next 


rocodure  INPUT-COEFFICIENTS  is 

— This  procedure  obtains  the  Information  about  tha  Input  polynoaial  and 
—  converts  tha  coaff iclants  into  floating  point  foraot 


package  INT-IO  is  nos  INTEQER_IO(  intagr); 
usa  IMT-IO; 


POWERS:  Intagar  :■  3; 

—Indicates  ehethr  all  poors,  only  tha  ouan,  or  only  tNi  odd  poors 

—  ora  prasant  in  tha  input  polynoaial .  Original  ly  sot  \‘a  out  of 

—  bounds  condition  to  urify  propr  Input. 

STEPS:  intogr  :■  2; 

— Incranant  ualua  for  antring  tha  coaff iei ants  of  tha  polynoaial 
INITIAL:  Intogr  :■  0; 

-•Starting  ualua  for  tha  ualua  of  tha  exponent 
COUNTER:  intagr; 

—Loop  counter  through  tha  input  string 
NUNERRTOR:  intagr; 

Nuorator  of  the  coefficient 
OENOHINRTOR:  intagr; 

— Oanoainator  of  tha  coaff icient 
CONUERT-STRING:  string  (1.  .HROLOIOIT); 

— String  representation  of  tha  coefficient 
LAST-DIGIT:  Intagr; 

—Actual  length  of  tha  Input  string 

begin  —Input  Coefficients. 

— Obtain  the  ualua  of  tha  largest  exponent  of  the  polymetal . 

—  It  oust  bo  be  tear  2  and  9. 

ehi  la  OEOAEEJOF-POLVNOniAL  <  2  or  OEDAEEjOF-POLWOHIAL  >  MAX-DEGREE  loo 
put  ("Enter  tha  degree  of  polynoaial  desired.  (Niniaun  Is  2>:  ">; 
gat  (OEGAEEJOF-PQLVNOn I AL ); 
nosullna; 
end  loop; 

—Obtain  an  indicator  for  tha  typo  of  tha  polynoaial 's  exponents 
ehi  la  POWERS  <  0  or  PQWER8  >  2  loop 
put  ("Enter  0  for  coefficients  for  ALL  poors  of  X">; 
naaul  ina; 

put  <"Ehtr  I  for  coefficients  for  000  poors  of  X">; 
nosullna; 

put  ("Enter  2  for  coaff  Iclants  for  EUEN  poors  of  X">; 


gat  (P0WER8); 
end  loop; 

— Sat  the  Initial  and  Incraaontal  ual 
—  coefficients.  Saues  tiao. 
if  P0WER8  >  0  than 
STEPS  :■  I; 
alsif  POURS  ■  1  then 
INITIAL  :•  1: 


for  obtaining  the  polynoaial 
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—Obtain  tha  ooaff iciants  for  aoch  alaaant  of  tha  polynomial 
put  ("Entar  tha  ooaff  iciants  of  tha  sari  as  being  expanded  by"); 
naaJIna; 

put  ("  entering  a  fraction,  i.a.  -2/3  or  +2/3  or  2/3“ >; 
n>iu|  |  • 

put  (“Coefficient  for  X**  “>; 
naaL.1  ina; 

— Loop  through  all  alaaants 
■hi  la  INITIAL  <■  OEGREEjOF-POLVNOHIAL  loop 
put  (INITIAL); 
put  <"  •  “); 

gat. lino  (CONUEATJBTRINB,LAST-DIGIT); 
nea-l  Ina; 

COUNTER  :■  1; 

—Stop  through  tha  input  string  looking  for  tha  “/“  ahlch  soparatas 

—  tha  nuaarator  froa  tha  danoalnator.  If  ona  doas  not  exist,  or  It 

—  appears  in  eithar  tha  first  or  tha  lost  position  In  tha  string, 

—  than  tha  coefficient  oust  be  reentered. 

■hi  la  COUNTER  <■  LAST-DIGIT  loop 

if  (CONUEHT-STRINO  (COUNTER)  ■  */'  >  and 
(COUNTER  /-  CONUERTJSTRING' first  and 
COUNTER  /-  LAST-DIGIT)  than 

declare 

—Once  the  “/“  has  boon  located  and  Is  In  a  proper  location 

—  obtain  tha  nuaarator  string  and  tha  danoalnator  string. 
NUNERRTQR-STRINB:  string  rananas 

C0HUERTJ3TRINQ  (CONUEAT.JSTRINB* first.  .(COUNTER  -  1»; 
DENOHINRTOR-STRINO:  string  rewMt 
CONUERT-STRINO  ((COUNTER  +  1).  .LAST-DIGIT); 

begin  — Block 

— Convert  tha  too  strings  into  integers 
NUMERATOR  :■  STRINOlTOlINT  (NUMERATOR-STRING); 

DENOMINATOR  :»  STRINOLTOLINT  (DENOMINAfTOR-STRINO); 

— If  the  danoalnator  is  a  valid  value,  than  ganarata  tha  floating 

—  point  value  for  the  coefficient 
If  DENOMINATOR  /■  0  t hen 

COEFFICIENTS  (INITIAL)  :■  float(NUNERATOR)  /  float(OENOMINATOR); 
— Incr— nt  to  the  next  alaaant  In  tha  polgnoaial. 

INITIAL  :■  INITIAL  ♦  STEPS; 

«nd  if; 

—Indicate  that  this  coefficient  has  bean  found  and  converted 
COUNTER  :■  LAST-DIGIT; 
and;  —Block 

end  if; 

—Point  to  tha  next  character  in  tha  input  string 
COUNTER  :■  COUNTER  ♦  1; 
and  loop; 
and  loop; 

and  INPUT-COEFFICIENTS: 


— Gonaroto  tho  Mtrix  of  tha  Tchabyohof  polgnooial .  Tho  procoduro  i 
— -  valuos  of  tho  Mtrix  olooonts  that  hauo  alroady  boon  found. 

—  Tho  algor!  tho  is  roeuraivo  in  that  roapoet. 

bog  In  — Coaputo  Tchobgshof  Polgnooial . 

— Tho  first  too  olononts  oust  bo  initial izod  to  alios  tho  foi losing 

—  paaaos  to  uao  thoo. 

TCNEBYSHEFJ’OLVNOHIflLS  <Q,0>  :■  1; 

TCHEBV8HEFJ«tJLVH0hlflLS  <t,t>  1; 

—Loop  ttrou^i  tho  looor  triangular  portion  of  tho  Mtrix 

—  and  oaloulato  tho  Tchobgrtwf  polgnooial  ualuM. 

for  I  in  2.  .HRXJEQflEE  loop 
for  J  in  0. . I  -  2  loop 
T€HEBV8HEFJP0LVNQI1lflL8  <I,J> 

TCHEBV8HEF_P0LVM0MlflLS  <I,J>  -  TCHEBV8HEF_P0LVM0h I flLS  (I  -  I 
and  loop; 

f or  J  In  0. . I  -  1  loop 
TCHEBVSHEF-POLVNOnmLS  <I,J  +  1> 

TCHEBV8HEF_PQLVN0niflLS  <I,J  ♦  1>  ♦ 

<2  «  TCHEBVSHEF-POLVNOflIRLS  (I  -  1,J»; 
and  loop; 
ond  loop; 

and  CtXVUTE-TCHEBVSMEF-PGLVrOllflL: 


procadura  C0WUTEJP0MER8-0F 
— Coaputa  tha  aatrix  for  tha 

COEFFICIENT-LIST:  FLOflTJUE 


INDEX: 

STEP: 


POINTER:  intagar; 

bag  in  — Coaputa  Poaars  of  Tchabyahaf . 
ahlla  INDEX  >■  0  loop 

MULTIPLIER  (INDEX)  1  /  TO€BV8*F_POLYNONIflLS  (INDEX, INDEX); 

STEP  :-  INDEX; 
ah  I  la  STEP  >*  0  loop 

COEFFICIENT-LIST  (STEP)  f lootOCHEBVSWF-POLVNONIflLS  (l«EX,8TEP» 

STEP  :■  STEP  -  1; 
and  loop; 

P0UER8L0F-.TCHEBV9HEF  (INDEX, INDEX)  :■  1.0; 

STEP  INDEX  -  2; 

■hi  la  STEP  >■  0  loop 
POUERSjOF-TCHEBVSHEF  (INDEX, STEP)  :- 
-  (COEFFICIENT-LIST  (STEP)) 

/  f loaKTCHEBVSHEF-PQLVNOMIflLS  (STEP, STEP)); 

POINTER  :-  STEP; 

■hi  la  POINTER  >-  0  loop 
COEFFICIENT-LIST  (POINTER)  :- 

COEFFICIENT-LIST  (POINTER)  ♦  P0UER8-0F-TCHEBVSHEF  (INDEX, STEP) 
•  f I ooKTCHEBVSHEF-POLVNOM I RLS  (STEP, POINTER)); 

POINTER  :•  POINTER  -  2; 
and  loop; 

STEP  :■  STEP  -  2; 
and  loop; 

irflEX  :•  IfCEX  -  1; 
and  loop; 
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28  Nooaahar  1989 

1.0 

f^jprcxJJHvar 

1.0 

This  routina  loops  until  a 
shichouar  function  ha  dasii 


Is  dona  approx iaa ting 


Choica 

MENU 

Copt.  8tauan  A.  Hotchkiss  and 

Capt.  Jam  I  far  Friad 

Oaualopad  as  a  thasis  otd  ADA  pro j act 


sith  GLOBAL-DATABASE;  usa  GLOBAL DRTnonoC; 

si th  APPROXIMATORS;  usa  APPROX I HAT0R8; 

sith  TEXT.IO;  usa  TECT-IO; 

procadura  APPRQXJ3AIUER  is 

HUH:  Intagsr  :■  0; 

DEN:  intagsr  :■  1; 

CHOICE,  KEY  :  choroctar; 

QUIT  :  choroctar  :-  •?*; 


package  INT-IO  Is  nas  INTEGEAJQ<  INTEGER); 
usa  INT-IO; 


pockaga  FLT-IO  Is  naa  FLOAT-I 0<LQNGlFL0RT >; 
usa  FLT-IO; 


bag  in 

sat_pagaulangth<24>; 

—  initialize  data  points 

camE-TQOvacv; 


—  let  tha  user  approxiaata  as  aany  functions  as  noodad 
■hi  la  (CHOICE  /■  QUIT>  loop 

—  select  function  to  approxiaata 

—  by  gluing  users  a  aanu  of  options 
NENU(CHOICE>; 

■—  use  the  built  functions  to  sake  a  sore  accurate  approxiaation 
CONPUTE-PADE-APPROXINATION8; 

CCNPUTEJCNK; 

If  CHOICE  /■  QUIT  than 
for  I  In  0. .N  loop 

If  C<NUH, I )  /■  0.0  or  C(Den,l >  fr  0.0  than 
put<*0***>; 
put(l); 
putC  '); 
put(C(NUN.  I »: 
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Data: 

Uarsion: 


Nom: 

Modulo  Nuabor: 
□ascription: 
Poassad  Uariablas: 
Ra turns: 

Globols  Usad: 
Modules  Cal  lad: 
Author: 

History: 


28  Novaabar  1 

1.0  _ 

QL0GRL_DHTHGft8C 

2.0 

Contains  all  global  uoriablas 

N/R 

H/R 

ALL 

H/R 

Copt.  Stauan  R.  Hotchkiss  and 

Copt.  Jennifer  Frlad 

Coop I a ted  for  Thasis  and  ADR  pro j act 


0L0BRLi)RTRBR8E  is 


LOHOLFLORT  is  digits  8; 

HECTOR  is  arroydntagar  rartga  0..2S>  of  LONGlFLORT; 

MATRIX  is  orraydntagar  ranga  0.  .25,  intagar  ronga  0.  .25)  of 
LONQ-FLORT; 

PAOE-URTRIX  is  array  (intagar  ronga  0..23,  intagar  ranga  0..1, 
intagar  ronga  0.  25>  of  LONB-FLORT; 


T:  HRTRIX; 

R:  PAOEJtITRIX; 


0:  HECTOR; 
h:  intagar; 

K:  intagar; 

H:  intagar; 
MACLRURIH:  HECTOR; 


COEFFICIENT:  stringCI.  .33); 

EPS:  LONO-FLORT; 

C:  HATRIX; 


—  Matrix  containing  tho  coafficionts  of 

—  diffarant  poaars  of  Tchabyshav  pojywnials 

—  Usad  to  contain  tha  sari  as  of  PACE  approx 

—  R(S,H  or  D,C> 

—  S  is  tha  sorias  nua bar 

—  H  or  D  M  -  0  for  tha  nuaarotor 

—  D  -  1  for  tha  danoainator 

—  C  -  coaff  iciant  for  a  pooar  of  X  for  tha 

—  particular  sorias'  nuaarotor  or 

—  danoainator _ 

—  Error  values  of  PACE  approx i nations 

—  Paaar  of  tha  nuaarotor  polynoaial 

—  Pooar  of  tha  danoainator  polynoaial 

—  Pooar  of  tha  initial  pooar  sorias 

—  Contains  tha  coafficionts  for  tha 

—  diffarant  poo**  of  "X"  for  tho  pooar 

—  sari os  axpcnslon  of  a  function 

—  Usad  to  contain  user  antarad  coafficionts 

—  a  pooar  sorias  expansion 

—  Convargant  apsilon 

—  Final  rational  approxlaation 


DATE: 

Version: 


28  Novestof  1885 

1.0 

ncM.  umiums 

nodule  Nuetoer:  3.0 

Description:  This  pockoge  contains  procedures  that  are  Invoked 

throughout  the  systee 
Passed  Variables:  N/A 

Returns:  M/fl 

Olabals  Used:  H/R 

Rodules  Cal  led:  N/fi 

Author:  Copt.  Steven  A.  Hotchkiss  and 

Copt.  Jennifer  Fried 

History:  Developed  as  a  thesis  and  flOfi  project 


th  OL08RL- DflTflBRSE;  use  QLOBRL-DRTRBRSE; 

skage  COIIIOILPROCS  is 

procedure  POWER-PROMPT <NUM,  DEN:  out  integer;  Epsi  Ian:  out  LONOlFLORT) 
procedure  OETjCOOTICIE?fT8<STRUtrrunE:  In  character;  POUER:  in  Integer); 
function  PROOUCT(FROM,  TO,  BV:  integer)  return  LONOLFLORT; 
function  FRCTORIRL  (NUMBER:  integer)  return  LONOLFLORT; 


package  body  COWOLPHOC8  is 


package  INT.IQ  is  noo  IHTEQER-KX  intogor); 
un  1MT-I0; 

package  FLT.I0  is  nss  FLQRT-I 0<L0H0_FL0flT ); 
us«  FLT-IO; 

proosdur*  POUER-PBOfPKMUM,  DEM:  out  Integer;  Epsi  Ion:  out  LOKLFLQRT)  is 
begin 

s*t_pog*_l ong th<24 >; 
loop 


put<"Entor  ths  pooor  of  ths  nu— ratoKoust  bo  intogor)  "); 
got  (Mill); 
noo-l  ino; 

put< ‘Enter  tho  pooor  of  tho  donooina  tor  (oust  bo  Intogor)  *>; 

get(OEN); 

neoul  ino; 

put<  "Enter  tho  opsilon  of  convergence  ■*>; 

putCThis  oust  bo  a  roal  fraction  and  ontorod  as  0.x"  ); 

put<"Uhare  x  is  any  string  of  digits  up  to  0  In  length  ">; 

got<EP8ILQN); 

nooul  Ino; 

* 

oxit; 
ond  loop; 
oxeoptlon 

ehen  datauorror  ■> 

put-1  ine<”  Invalid  Entry.  Roan  tor  data"); 


procedure  QET-CCEFF 1C IENT8 (STRUCTURE:  in  character;  POUER:  in  integer)  is 


COEFF  :  LONQlFLORT  O.O; 

FRON,  TO,  BV  :  intagar; 

procadura  GET-POUER<NUNBER:  in  Intagar;  COEFF:  out  LONQlFLORT)  Is 

LRBTJBMRP  :  intagar  COEFFICIENT'  last-1; 

MIERRTOR  :  boo  I  son  TRUE; 

OUT-COEFF  :  LOHOlFLORT; 

CHRR-PTR  :  intagar; 

Mil,  DEN  :  stringCI. .  19); 

INPUT-ERROR  :  exertion; 


procadura  COHPUTEJERLJXEFF<NUN,  DEN:  in  string; 

COEFF:  out  LONG-FLORT)  is 

CHRR-PTR:  intagar  :-  NUN' first; 

MfCRRTOR:  UONGlFLORT  :•  0.0; 

OENOniNRTOR,  SIGN:  LONQlFLORT  :-  1.0; 


bagin  —  CONPUTEJRERLJCOEFF 

if  <NUN<CHRR_PTR)  -  •♦*>  than 
QWLPTR  :■  CHRR-PTR  ♦  1; 
alsif  <NUH(CHHR_PTR>  -  •-*>  than 
SIGN  -SIGN; 

CHflR_PTR  :■  CHRR-PTR  +  1; 
and  if; 

■hi  la  ((MIKCHHRLPTR)  b  '  ’  >  and  (CHRR-PTR  <•  MIT  last))  loop 
NUNERRTQR  NUNERRTQR  *  10.0  ♦ 

LONG-FLORT (chcractor 1  pos<NUN<CHRR_PTR  >  )  - 
charaetor‘pos<*0* )); 
and  loop; 

CHRR-PTR  :•  OEN'first; 
if  (OEMCHRR-PTR)  -  >  than 

CHRR-PTR  :■  CHRR-PTR  ♦  1; 
alsif  COEN(CMVLPTR)  •  •-*)  than 
SIGN  :■  -SIGN; 

CHRR-PTR  CHRR-PTR  ♦  I; 
and  if; 

•hi  la  (COEN (CHRR-PTR)  >  and  (CHRR-PTR  <■  OEN'Iast))  loop 

NUNERHTOR  :■  HUNERRTOR  *  10.0  + 

LONO-FLORT(charactar'pos(OEN(CHRR-PTR))  - 
character 'pos( '0* )>; 
and  loop; 


COEFF  NUNERRTQR /DENON I NRT0R*8 1 GN; 


begin  —  QET-POEB 


loop 

—  proapt  the  user 

put (“Enter  the  coefficients  for  x**“>; 

put(NUNBER); 

putt*  — >  ■); 

gat  (COEFFICIENT); 


—  pock  and  asperate 

for  I  in  COEFFICIENT 'range  loop 

if  <'0‘  <■  COEFFICIENT < I  )  end  COEFFICIENTS  >  <-  ‘9‘>  or 
COEFFICIENTS)  ■  or  COEFFICIENTS)  ■  or 
COEFFICIENTS)  -  '/'  or  COEFFICIENTS)  •  '  *  than 

if  COEFFICIENTS)  -  */*  than 
CHHLFTR  :■  OEN' first; 

HliWMITflft  * ■  CflIOC • 

elslf  COEFFICIENTS)  -  or  COEFFICIENTS)  -  *-•  or 

CO*  <.  COEFFICIENTS)  and  COEFFICIENTS)  <-  ‘9*) 

then 

If  NUHERffTOR  then 

NUN<CHHR_PTR)  COEFFICIENTS); 

else 

OEH(CHflR_PTB)  :■  COEFFICIENTS); 
and  if; 

CHRRLFTR  CHflR-PTR  ♦  1; 
end  If; 

else 

raise  INPUT-ERROR; 
and  If; 
end  loop; 

exit; 
end  loop; 

COrPUTEJERLJCOEFFCNUn,  OEN,  OUTJCOEFF); 

COEFF  OUTJCOEFF; 
put  (OUTJCOEFF); 
nee-line; 

exception 

ahan  INPUT -ERROR  ■>  put-1  IneS  Input  Error.  Reenter  value.  ">; 
nesJ  ine; 


gin  —  GET-COEFFICIENTS 
set-pageJongth<24>; 


put-1  ine<‘Eniar  the  coefficients  for  eoch" ); 
put— I  ina(*pcNMr  of  tha  ‘X’  in  fractional  fcre.“>; 
put— I ina("lf  a  sign  is  entarad,  it  oust  bo  tha  ">; 
put-1  ine<" first  character.  No  blanks  are  al  loaad.  * > 
put-1  lna<*7ha  aax  allaaabla  size  is  9  digits  par* ); 
put-1  ina<"nuobar.">; 
wmLmI  ins* 

put— I  ine(a8onple  tries:  1/2  ,  +1/2  ,  or  -1/2"> 
nae-l  ine; 


TO  :•  POWER; 
case  STRUCTURE  is 


ahan  T  »  FRON 

■  0; 

BV 

■  1; 

ahan  *2*  ->  FROM 

■  0; 

BV 

-  2; 

ahan  *3‘  ■>  FRON 

•  1; 

BV 

•  2; 

ahan  othars»  FRON 

■  0; 

BV 

-  i; 

end  case; 


ehi  la  (FROM  <■  TO)  loop 
GETJWBKFROH,  COEFF); 
NRCLRURIN<FRON)  COEFF; 
FRON  :■  FROU  ♦  BV; 
end  loop; 


end  0ETJXEFFICIENT8: 


V.IT^  V^-  »;-  •».-  »■.•  -.-  tf’  o  -l-  H-  <■ 


u 


Data: 

28  Nouaobor  1885 

Uerslan: 

1.0 

Naoo: 

FUNCTION-PACKAGE 

Modulo  Nuobor: 

4.0 

Ooocription: 

This  pockags  contains  aoduias  that  era  cal  lad 
cooputo  a  predefined  poser  series  expansion  of 
or  alloe  a  users  to  enter  their  am 

Passed  War  icfcles: 

N/R 

Ro  turns: 

N/fl 

Olobals  (tood: 

ousRL-orrmonoc 
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a  function 


package  FUNCTION-PACKAGE  is 
procedure  SIN-SERIES; 
procedure  TRN-SERIES; 
procedure  AS IN-SERIES; 
procodiro  ATRN-SERIE8; 
procedure  EXP-SERIES; 
procodiro  BUILD-SERIES; 


«nd  FUNCTION-PACKAGE: 


T 


t 


package  body  RJNCTIOtLPRCKRQE  is 


procadura  SIIUSERIES  Is 
N:  intogar; 
bogin 

—  gat  tha  poaara  of  tha  nuaarator  and  donooinator  polynomials 

—  Also  proapt  tha  usar  for  a  conMorgont  apsl Ion. 

pa€R_pnofrr<n,K,EP8); 


—  Cooputa  tha  poaar  of  tha  noclaw'ln  sari  as.  It  Is  tha 

—  of  tha  nuaarator,  danoainotor,  and  tha  uolua  too 
H  :■  H  ♦  K  ♦  2; 

—  Cooputa  tha  initial  approximating  polynoaial 
for  I  In  0. .29  loop 

hRCLflURIIKI  >  0.0; 

and  loop; 

for  I  in  1..«hf1)/2>  loop 

mCLflURIH<  1*2-1 )  -1.0**l /FHCT0RlflL<2*l-1>; 

and  loop; 


and  3111-SERIES: 


procedure  7RNUSERIE8  is 


N:  integer; 
begin 

—  gat  tha  poeari  of  the  numerator  and  danoainator  polynoeial*. 

—  Also  praapt  tha  uaar  for  a  convergent  apsi  Ion. 
P0UER_PA()VT(n,K,EP8>; 

—  Coaputa  tha  poaar  of  the  hoc  I  car  in  series.  It  is  tha  sue  of  tha 

—  of  tha  numerator,  danoainator,  end  the  value  too 
N  :■  H  ♦  K  +  2; 

—  Coaputa  tha  Initial  approx  looting  polynoalal 
far  I  In  0. .29  loop 

mamiRimt)  o.o; 

end  loop; 

mCLflURINd)  :■  1.0; 

far  I  in  1..«N»1>/2>  loop 

HACLflUAIM<2*l+1)  :■  PA00UCT<2,  2*1,  2)  / 

FACT0RlflL<2*l+1>; 

and  loop; 
and  TAN-SERIES: 


procedure  R8INJSERIES  is 


N:  integer; 
begin 

—  gat  the  posers  of  the  numerator  and  denominator  polynomials. 

—  Also  proept  the  user  for  a  convergent  epsilon. 
POBLPWflPT<H,  K,  EPS  ); 

—  Coepute  the  poser  of  the  floe  lour  In  series.  It  is  the  sub  of 

—  of  the  nueerator,  denoeinator,  and  the  value  tso 
N  :■  H  ♦  K  ♦  2; 

—  Coepute  the  initial  approx i eating  polynomial 
for  I  In  0. .29  loop 

fflCUURIIKn  :•  0.0; 
end  loop; 

for  I  in  t..«NM>/2>  loop 

HACLflURIN<  1*2-1)  :•  PROOUCT<1,  «l-2)*2+1>,2>  / 

PBOOUCT<2,<  1*2-2), 2>  • 

LONILFLQHK  1*2-1  >; 

end  loop; 
end  R8 IN-SERIES: 


procadura  ATRKjSERIES  is 
N:  intagar; 
bag  in 

—  gat  tha  poaars  of  tha  nuaarator  and  danooinator  polynoaiols. 

—  Also  proapt  tha  usar  for  a  convargant  apsi Ion. 

PCNER-PROHPT  <H,  K,  EPS  ); 

—  Coaputa  tha  poaar  of  tha  hoc I our in  sari as.  It  is  tha  of  tha  pooar 

—  of  tha  nuMrator,  danooinator,  and  tha  volua  too 
M  :■  tl  ♦  K  ♦  2; 

—  Coaputa  tha  initial  approx  (noting  poignoaiol 
for  I  In  0. .25  loop 

mCLflURIIKI)  :•  0.0; 
and  loop; 

for  I  in  1..«N*1)/2>  loop 

f«CUURIN<2*l-1>  H.O*<M>/FfiCTORIflL<2*l-t>; 
and  loop; 


and  RTRNLJ8ERIE8; 


procadum  BUILDJSERIE8  is 

N  :  intagar; 

STRICTURE  :  choroctar; 

bag  in 

sat_pagn_langth<24); 

—  gat  tha  poaars  of  tho  nuarator  and  danoainaior  polyioaiais. 

—  Also  proapt  tha  usar  for  a  oonMorgant  apsl  Ian. 
P0tER_PAQflPT<H,K,EP8); 

—  Coaputa  tha  pooar  of  tha  hoc I our in  sari as.  It  is  tha  sub  of  tha 

—  of  tha  nunarator,  danoainotor,  and  tha  oalua  too 
h  :■  n  ♦  K  ♦  2; 

—  Coaputa  tha  initial  appraad noting  polynoaial 
for  I  in  0. .29  loop 

mCLAURIH<l>  0.0; 
and  loop; 

—  Proapt  tha  usar  for  tha  structura  of  tha  polynoaial 


LI: 

loop 

put_l  inaCEntar  1  If  all  poaars  of  X">; 
put-1  Inot’Entar  2  if  only  auan  poaars  of  X“>; 
put<*Entar  3  if  only  odd  poaars  of  X  — >  ">; 
gat  (STRICTURE); 
naa-llno; 

if  *  1'  >  STRUCTURE  or  STRUCTURE  >'3‘  than 
put_l  inaCBod  Entry.  Try  again.  ">; 

alsa 

OET_COEFFICIEMTS<  STRUCTURE^); 
and  If; 

axit  LI; 

and  loop  LI; 


and  BUILD-SERIES; 


i  ^TATOTT. 


MICROCOPY  RESOLUTION  TEST  CHART 

NAtmNAI  RtIBFAU  OF  STANDARDS  196‘  6 


EXPJBERIE8  is 


N:  in 


bag  in 

—  gat  tha  poaars  of  tha  nuarator  and  daraninator  polytoaials 

—  Also  prcapt  tha  usar  for  a  conwargant  apsi  Ion. 
POSLPM)nPT<n,K,EP8>; 


—  Caaputa  tha  posar  of  tha  hoc I our in  sari  as.  It  is  tha  sua  of 
—  of  tha  m— "ator,  danoainotor,  and  tha  uolua  too 
N  :■  ft  ♦  K  ♦  2; 

—  Caaputa  tha  initial  appro* I acting  polgnoaiol 
for  i  in  0.  .29  loop 
NKLNfflllK!)  :•  0.0; 
and  loop; 

for  I  in  0. .N  loop 

mCLflURirKI)  1.0  /  FnCTQRIflLd  >; 
and  loop; 

and  EXPJBBIIES; 
and  RJNCTION-PfCXAOE: 


/*/ 


—Data: 

—  Uersion: 

—  Noee: 

—  Nodule  Hueber: 

—  Description: 

—  Pomd  Uariables: 

—  Returns: 

—  Olobals  Used: 

—  Hodul—  Col  lod: 

—  Author: 

—  History: 


28  How— hor  1989 

ro _ 

APPAQXINRT0A8 

9.0 

This  package  Includ—  prooodur—  that  coapute 
approx  I —t  ions  to  uoor  —I— tod  functions 
N/A 

N/A  _ 

QLQBAL-DATABABE 

N/A 

Copt.  8 town  A.  Hotchkl—  and 

Copt  Jennifer  Priad 

Dowolopad  —  a  thesis  and  ISA  project 


package  APPROXIMATORS  is 

procedure  C0NPUTE.TCHEBV9HEU; 
procedure  fOtK  CHOICE:  out  character); 
procedure  C0fWTE_PA0E_APPA0XINATI0N8; 
procedure  COTPUTEJOK; 


end  APPROX  I HATORB: 


■Ith  GLOBAL-DKTAMSE;  un  GLOBAL-DATABASE; 

•ith  FUCT  ION-PACKAGE;  uat  FUNCTION-PACKAGE; 

•Ith  TECT-IO;  un  TECT-IO; 

poduaga  body  APPAOXIIWTOAS  is 

package  FLT.IO  is  noo  FUWT_IO<U3hOLFLtWT);  un  FLT-IO; 

procdtrt  COWTi_TCHEBW«U  Is 
bogln 

—  bui  Id  tho  global  table  "T*  containing  the  coofficiants  for 

—  each  of  a  sari  os  of  Tchobyohov  polipnoiais 

T<0,0>  :■  1.0; 

T< 1, 1 >  1.0; 

T<2,0>  :■  -1.0; 

T<2,2>  2.0; 

for  I  In  3.  .25  loop 

for  J  in  0.  .25  loop 

T<I,J)  :■  T<I,J>  -  T<I-2,J>; 
end  loop; 

for  J  In  0. .24  loop 

T<I,JM>  T<I,JM>  ♦  2.0  •  T<t-l,J>; 
and  loop; 

and  loop; 


and  COHPUTE-TOOVBHN; 


■VWJWWT*'. 


MENU  (CHOICE:  out  cho-acter)  is 


OUTJCHOICE:  chorocUr; 
BRDUCHOICE:  SMoaption; 

bog  in 

sat-poge-l  ongth(24  >; 

—  cl  oar  screen  and  print 


put-1  I no< "Choose  function  to  bo 
nooJino; 


put-1  InoOEntar  1  far  sin”); 
put— 1 1  not*  Enter  2  far  ton'); 


looted"); 


put-i  ino<  "Enter  3  far  arcs  in"); 

put-1  ino<  "Enter  4  far  arc  list"); 

put-1  ina< "Enter  S  far  asp"); 

put-1  (na<  "Enter  0  for  uaor  defined  function") 

put-1  Inot-Ehter  ?  to  quit"); 


far  I  in  HRCUUR  IN’ range 
HRCLRURINd  )  0.0; 

and  loop; 


noai.1  Ino; 
put<*— >  "); 
gat<0UT-£H0ICE); 

CHOICE  :■  OUTJCHOICE; 
caso  OUTJCHOICE  is 

ahan  T  ■>  8IN-SERIE8; 

ohan  ‘2*  ->  TRNL8ERIE8; 

ahan  ‘S'  •>  RBIILSERIE8; 

ahan  '4‘  ■>  RTRNJSERIES; 

ahan  ‘5‘  ->  EXP-SERIES; 

ahan  ‘O'  •>  BUILOLSERIES; 

ahan  *?'  ■>  null; 

ahan  othars  ■>  raise  BROlCHOICE; 


oxl  t; 


exception 

ahan  MLCHOICE  » 

put-1  lno("lnwal id  entry.  Try  again"); 


COWUTtFROEJFFROXINnTOHB  Is 


NJffK:  integer; 

Nil  :  integer  :■  0; 

OEM  :  Intag**  :•  1; 

TEMP  :  LOHOlFLOHT; 

UOFK  :  HRTRIX; 

B  :  VECTOR; 

bagln 

—  this  prooadm  converts  tha  initial  appratclaating  polymetal 

—  (tha  Haelaurin  peaar  sari  as)  Into  a  rational  approx I eat I on 

—  cl  oar  out  this  PROE  approx I eat I  on's  nuasrator 

—  and  dvMMlnator  polynoaials 
for  SERIES  In  0.  .25  loop 

far  NUHLOEN  In  0. .  1  loop 

for  COEFFICIENT  in  0..25  loop 

R<SERIES,NUHJEN, COEFFICIENT)  :■  0.0; 


for  I  In  0.  .H  l< 


—  loop  for  al  I 


far  J  In  0.  .K  loop  —  loop  far  al  I 


of  tha 


of  tha  donoainotor 


If  <1  >-  J>  than 

—  build  a  work  aatrlx  to  solve  siaul taneous  aquations 
B<0)  :•  1.0; 

NJWX  :■  I  ♦  J; 

far  3  In  0..<tUMX  -  I  -  1>  loop 
far  HI  in  0. . J  loop 

U0RK(8+1,N1)  HKlRURIN<abs<NJ«X  -  8  -  H1»; 
If  (HI  -  0)  than 

B<8+1)  :•  -HRCLRUR I N(abs<NJtfK  -  S  -  N1»; 
and  if; 
and  loop; 


—  Solve  siaul  taneous  aquations  for  donoainotor  coefficients 
for  HI  in  1..J  loop 

If  OOK(H1,N1>  ■  0.0)  than 
SETUP: 

for  N2  in  1..J  loop 

If  <U0n«N2,N1V*  0.0)  than 
TEIP  :■  B(N2>; 

B(N2>  B(N1>; 

B(N1)  :■  TEMP; 
for  HS  in  1..J  loop 
TEMP  U0RK(N2,N3); 

U0HKN2,N3)  M0fK(N1,N3>; 

U0NKN1,N3)  :■  TEMP; 
and  loop; 
exit  8ETUP; 
and  if; 

and  loop  8ETUP; 


TEMP  :■  U0NKH1,  Ml); 

If  TSP  b  0.0  than 
B<M1>  :■  B<N D/TEMP; 

also 

B<M1>  :■  0.0; 
and  if; 

for  M2  in  1.  .J  loop 
if  mr  b  0.0  than 

tOK(N1,N2)  iO«C(N1,N2>/TBV; 

•Im 

MOV((N1,N2)  :■  0.0; 
and  if; 
and  loop; 


for  M2  in  1.  ,J  loop 
If  <N1  b  M2)  than 
TEMP  -U0fK(M2,N1); 
far  N3  in  1..J  loop 

U0fVKN2,M3)  UGRK(N2,M3)  ♦  U0RK<M1,N3)  *  TEMP; 
and  loop; 

B<N2)  :■  B(N2)  ♦  B<N1)  •  TEMP; 
and  if; 
and  loop; 

and  loop; 

—  un  danoainotor  ooofftcionta  to  oooputo  tha  nu— rotor 

—  coofficiants,  and  build  tha  sar  las  of  PME  approx  loot  Ians 
for  Ml  In  0. . I  loop 

far  N2  in  0..Nt  loop 

R<l«d,NUH,N1)  :•  R<l+J,MUn,M1>  ♦  B(N2)  •  HRCLRURIH<N1-M2)/ 
B<0); 

and  loop; 
and  loop; 

for  N1  in  rovarsa  0..J  loop 
R<I+J,0EN,M1)  B<M1)/B<0); 

B<Nt)  B<N1)  /  B<0); 
ond  loop; 

—  Coaputo  tha  0's  that  arm  uood  to  oooputo  C(a,k) 

—  0<l+J*t>  -  8UMIL-0  to  J  <Hoelaurln<l+OH-0*8<L)l 

IKI+Jtl)  0.0; 

for  L  In  0. .  J  loop 

0<I«J»1)  0<I+Jf1>  ♦  MRCLRURIMO+vHI-L)  •  B<L); 

and  loop; 

and  If; 
and  loop; 
and  loop; 


and  CGMPUTE_PflGE_flFPRQXIMRTIOMB: 


adurt  CCWUTEjCHC  is 

A:  intsgar  :■  0; 

B:  intagar  :■  1; 
LftDA:  VECTOR; 

bag  in 


—  Cosputa  tha  1  (aipha^l) 

UTOACO)  -CDCIWC+1)  •  TCJWC,0»  /  C2.0**CWHO>; 

far  J  In  O..Cft*K-1>  loop 
if  DCJ»1>  f  0.0  than 

LAHDACOH)  :■  CDCmttH)  •  TC!WCM,JM»  /  «2.0  **CIWC»  •  0<J+1» 

alas 

UBOACJH)  :•  0.0; 
and  if; 
and  loop; 


—  Load  Pa<X)  and  <*CX>  oith  thoir  A  and  B  ooofficiants  raapactivaly 
for  I  in  0.  .11  loop 
CCA, i >  :■  ACn*K,A,l); 
and  loop; 


for  I  in  0. .K  loop 
CCB,  I )  :•  ACIHK,B,  I ); 
and  loop; 

—  Caaputa  ooaff iclonts  "A"  of  nuaarotor  and  ”B"  of  danaainator 
far  J  In  0..CIWC-I)  loop 
for  K  In  0. .29  loan 

ACJ,A,K>  ACJ,A,K>  *  LfWACJH); 

CCA, K>  :•  CCA,K>  ♦  ACJ,A,K>; 

RCJ,B,K)  :»  RCJ,B,K)  •  LHCACJM); 

CCB, K)  CCB,K>  ♦  RCJ,B,K>; 

and  loop; 

and  loop; 

CCA,0)  CCA,0)  ♦  LflADACO); 

for  I  in  rowans  0..29  loop; 

CCA,  I >  CCA, I )/CCB,0>; 

CCB,  I )  CCB, I >/CCB,0); 

and  loop; 
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This  preset  produced  a  run-time  math  library  for  the  MIL-STD-1750A 
embedded  computer  architectures.  The  math  library  consists  of  the 
circular  trigonometric  functions  and  their  inverses.  In  addition,  the  steps 
required  for  the  performance  analysis  of  the  math  library  have  been 
outlined. 

Several  approximation  methods  were  investigated,  but  the  Chebyshef 
Economization  of  the  Maclaurin  series  polynomials,  and  rational 
approximations  derived  from  the  second  algorithm  of  Rernes  were 
determined  to  be  the  best  available.  Each  functions  implementation  was 
designed  to  take  advantage  of  features  of  MIL-3TD-1750A  architectures. 
The  recommended  test  procedures  will  provide  measures  of  the  average 

and  worst  case  generated  errors  within  each  approximation. 
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